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Over  the  course  of  this  contract  a  variety  of  techniques  have  been 
employed  to  study  the  properties  of  small  atomic  and  molecular  clusters 
formed  in  the  gas  phase  via  homogeneous  nucleation.  The  clustering 
occurs  either  in  an  adiabatic  expansion  of  a  condensable  species  (e.g. 
argon,  krypton,  xenon,  or  sulfur  hexafluoride)  in  an  inert  carrier  gas 
(e.g.  helium),  or  as  a  mixing  process  using  a  hot  condensable  (e.g.  lead, 
silver,  copper.  Indium  or  bismuth)  and_a  eojd  cajrrier  gaj _ _ 
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helium,  carbon  dioxide  or  sulfur  hexafluoride).  A  continuous  development  ot 
several  types  of  cluster  sources  has  been  carried  out  and  includes  free 
jets,  very  small  hypersonic  laval  nozzles,  and  a  series  of  metal  ovens 
with  carrier  jjas  mixing.  Any  one,  of  these  sources,  constitutes  the  first 
stage  of  a  differentially  pumped,  molecular  beam  system  which  then  produces 
a  continuous  beam  of  clusters.  The  density  in  the  beam  is  so  low  that  it 
is  collisionless  and  thus  the  clusters  do  not  interact  with  each  other 
or  with  any  other  foreign  molecule  or  surface.  The  study  of  these  isolated 
clusters  is  carried  out  primarily  using  high  energy  electron  beams  (40  to 
75  KeV) .  The  resulting  diffraction  patterns  are  obtained  either  on  film 
or  through  use  of  a  single  channel,  scintillation,  pulse  counting  system 
employing  synchronous  detection.,-- The  diffraction  patterns  after 
interpretation,  provide  information  on  the  structure  of  the  clusters,  as 
well  as  an  estimate  of  their  size  and  temperature.  The  size  range  studied 
extends  from  50,000  down  to  50  atoms  per  cluster.  The  goal  is  always 
to  push  toward  the  smallest  size  possible  and  to  ascertain  whether  changes 
in  structure  occur.  It  there  is  a  change  in  structure  their -the  nature 
of  the  change  and  the  size  range  at  which  it  occurs  is  determined.  It 
is  found  that  as  cluster  size  is  reduced  there  were  indeed  structural 
changes  such  as,  the  transition  from  face  centered  tetragonal  to  face- 
centered  cubic  (fee)  in  indium,  the  appearance  of  icosahedral  structure 
in  noble  gas  clusters  and  the  appearance  ol  a  new  *amorphous^  phase  in 
metal  clusters  such  as  lead  and  silver.  Computer  programs  were  written 
to  determine  the  extent  and  magnitude  of  changes  to  the  diffraction 
patterns  due  to  the  possibility  of  electron  multiple  scattering  within 
the  cluster.  These  effects  are  significant  for  clusters  of  high  atomic 
weight  species  and,  of  course,  for  very  large  sizes.  Finally,  an 
investigation  was  undertaken  (but  not  yet  completed)  to  examine  the 
similarities  and  differences  in  the  structural  changes  that  occur  in 
the  noble  gas  clusters  Ar,  Kr  and  Xe ,  and  in  the  metals  Pb  and  Ag;  all 
of  which  have  f.c.c.  bulk  structure. 
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ATOMS  AND  MOLECULES  IN  SMALL  AGGREGATES 
THE  FIFTH  STATE  OF  MATTER 


Atoms  dnd  molecules  in 
small  aggregates 

The  fifth  state  of  matter* 


I  n  our  earliest  encounters  with  the  study  of  physics  we  were  told  that 
■  n  liter  could  he  subdivided  into  three  phases  solid,  liquid,  and  gaseous,  with 
physical  leatmes  that  made  them  distinct.  As  time  went  on  we  encountered  a 
loutlh  s', ile  ol  matter  that  of  plasmas  existing  naturally  in  stellar  and  inlet 
stellai  media  ami  artificially  in  numerous  research  laboratories.  Plasmas  consist 
of  a  mixtuic  of  positively  and  negatively  charged  particles  which  indeed  exhibit 
piopeihcs  very  distinct  from  those  of  the  more  familial  three  states  of  matter 
Now  tlnne  is  a  tilth  possibility.  Those  of  us  who  have  been  pursuing  the  study 
ol  atoms  and  molecules  itt  small  isolated  aggregates  or  clusters  have  made  physr- 
i.il  measrnemenls  on  these  systems  revealing  properties  significantly  different 
liom  those  of  the  bulk  phase,  so  much  so  that  one  might  deem  it  appropriate  to 
designate  matter  in  this  domain  as  the  fifth  slate  of  matter.  1  veil  when  clusters 
oi  panicles  are  not  in  isolation  their  properties  may  still  vary  appreciably  from 
bulk  as  in  lire  case  of  colloids  of  sufficiently  small  si/e. 

It  will  become  apparent  as  we  pursue  the  subject  at  hand  that  the  physics 
ol  the  boundary  between  the  three  traditional  phases  ti  e.,  interfacial  phenom¬ 
ena)  will  become  an  omnipresent  feature  in  the  study  of  small  clusters.  When 

•  me  considers  clusters  small  enough  so  that  they  are  essentially  all  surface,  one 
should  not  be  surprised  to  encounter  physical  phenomena  which  arc  at  least 
atypical  ■>!  the  bulk  phase  and  hard  to  calculate. 

Small  clusters  arise  naturally  as  an  important  entity  in  the  treatment  ot 
the  kinds  s  ot  phase  change  and  m  nearly  all  cases  elusleis  are  foum-d  .is  the 
tcsull  ol  .1  phase  change,  finis  an  interesting  interdepends nu-  arises  in  which 
the  physual  pioperdes  til  small  clusters  arc  nnporl.  nl  m  ihcoiclical  predictions 
loi  the  kinetics  ol  phase  change  and  conveisely  an  un  '.-islanding  ol  the  kinetics 
ol  phase  ,  lunge  is  Important  for  the  production  of  clusters  of  a  given  mean  si/e 
loi  use  m  the  measurement  of  then  properties  I  or  example  the  predicted  rale 
ol  loimahon  ol  water  droplets  in  a  condensing  vapor  li  e.,  ilx  nuclealion  rale) 
may  vary  by  as  much  as  a  factor  of  ten  billion  tl()IH  I  for  a  30'!  variation  in 
dioplel  suit. ice  tension.  Understanding  how  the  kinetics  of  water  nuclealion 
lake s  place  piovules  the  insight  necessary  to  produce  high  concentrations  ot 
.  litstc’t  s  ol  spec  it  K  mean  si/e  and  as  narrow  a  si/e  clistnbutnm  as  possible.  I  his 
is  impoilanl  loi  any  experiment  designed  to  study  cluster  properties 

Neatly  every  technique  for  the  formation  of  small  clusters  necessitates 
the  c  n  ation  ol  an  environment  where  the  species  to  be  clustered  is  highly  super 
salinated  llns  is  the  single  most  important  parainelet  tor  the  product  inn  ot 
high  c  one  c  ni rat  ions  of  clusters  of  very  small  si/e.  Small  duster  samples  whethei 
in  the  gaseous.  Iiquicl.  or  solid  stale,  or  on  substrates  (supporting  surlaees)  have 
been  cvammrd  with  photons  as  a  piobe.  such  as  light  sc  altering  and  absoiption 
\  lew  ol  these  I  cell  n  upics  include'  such  things  as  absoiption  spectroscopy  .  Iiglil 
■calle  imp.  both  clastic  (K.iylcigli  sc.illcimgl  and  inelastic  iK.inian  seal  Icrmg  I 
\  1 1  I  In-  modi-in  Ice  h  niqucs  ol  tunable  continuous  wave  KWi  and  pulsed  illu 
mmuf  foil  *  u  rltf  laser  pulse-  tunes  as  short  as  nanosec  ond  to  p  n  osec  curd  that  i- 
It)'*  to  In1'  seel  are  luricnlly  bemg  humph!  to  trial  in  tin  study  ot  these 
small  c  Instils  in  addition,  photons  are  being  used  to  ioni/c  (In  c  luster  spec  res 
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l.o  siihsci|iicnt  mass  analysis  using  mass  spectrometers. 
(  lusters  ate  also  hemp  examined  with  electrons  as  a  probe 
using  huih  loss  and  high-eticrgv  eleclton  diffraction  and  a 
s.incts  ot  electron  microscopic  techniques.  Vets  shorl- 
waselcnplh  photons  |  x-ray  si  have  also  been  used  to  study 
tin  slimline  ot  small  clusters  Magnetic  and  electrte  fields 
has c  also  been  useil  lot  examination  with  methods  such  as 
nu.  leal  mapnelic  tesonance  ansi  electron  spin  resonance 
I  malls  neutrons  have  been  ttsesl  as  ptohe  particles  to  study 
the  siiucture  ol  small  clusters. 

Small  cluster  formation 

\i  the  outset  one  might  ask  how  small  sloes  the 
[>.iii is  le  have  to  be  before  it  is  considered  to  tie  a  cluster 
lias  mp  propci  ties  distinct  from  the  hulk  material  ’  And  why 
ate  small  clusters  important  for  the  understanding  of  the 
■Is  n. nines  of  phase  change?  Let  us  begin  by  recalling  a  less 
is  alnics  from  classical  or  bulk  macroscopic  tliermodynum- 
i>  s  Wc  recall  that  a  single  component  thermodynamic 
'S  stein  could  he  completely  described  by  the  specification 
ni  t It t»v  properties.  Lor  example,  the  volume  of  ilu-  system 
s  1 111  LI  he  considered  as  a  function  ol  pressure,  temperature, 
and  mass,  I  ,  =  m  i  )  Now  it  this  system  is  cut  in 

hall,  each  half  may  be  described  by  I  .  -  r2t?’.  /,  m.h 
wheie  in  *  ==  m,/?.  and  we  see  that  the  system  is  still 
spec  tiled  with  three  variables,  However,  if  one  continues 
this  process  of  cutting  the  system  in  hall,  we  will  see  that 
at  some  point  the  system  begins  to  display  some  addi¬ 
tional  1'iopeilies  which  were  nol  noticed  previously  Lot 
,  sample'  when  the  system  si/e  begins  to  pet  small  wo 
noli,  c  1  h .1 1  the  components  begin  to  slick  ‘.ogether  or 
no.  to  solid  surfaces  or  nearby  walls  (  le.nly  some  olliet 
ihvim.idvnainic  properly  is  coming  into  play  ll  one  were 
to  lontmuc  this  halving  process  even  tin  I  lie  ■  Ihctc  would 
•  ”iii  a  |>omt  .it  which  (be  small  paiticlcs  would  begin  to 
li 'date  ni  llo.it  .lb,. ni  m  the  medium,  be  it  lli|Uld  ot  gave 
oils  displacing  Htiiwnian  motion  Again,  wc  are  beginning 
i"  obse'rve  another  physical  phenomenon  that  evas  not 
ipp.it cut  in  the  large  syslerns  In  mdei  In  di'linc  uiot|iicly 
the  .  lassie  at  Iflci  mu  uly  n. mill  slate  nl  the  system  ,n  hlilion.il 
p ropy; t ms  m ti si  he  specified,  lor  example,  st/c.  m  suil.ne 
hue.  .a  suil.ne  free  eiieigv  Ivpi.allv  win'll  paiticlcs 


approach  the  It)  to  100  jnn  I  I  jam  -  It)  h  ml  si/c  i.nigc 
they  will  begin  to  exhibit  the  above  mentioned  icatims 

L.el  us  now  constde'i  another  process  Heated  n 
introductory  physics  or  thetmody  namics  lutagmc  a  e.is 
enclosed  in  a  c>  lunler  wheie  a  piston  is  tree  to  move  an, I  i, 
subject  to  constant  pressure  As  heat  is  removed  lion,  the 
gas  the  piston  slowly  descends  until  at  some  point  phase 
change  begins  and  liquid  appears  in  the  bottom  ol  tlu 
cylinder,  litil  if  one  actually  portotms  this  cxpcutm'm  n 
is  seen  thal  the  ln|ind  phase  may  not  api'cat  at  the  v.ip"i 
lic|utd  equilibrium  lemperatuie  and  the  gas  will  iti'ic.nl 
supersaturate  Several  implicit  conditions  arc  icquinc  ni 
order  to  roali/e  I  he  idealized  phase  change 

I  In  order  I'm  the  phase  change  to  on  an  and  th. 
liquid  to  form,  a  surlace  onto  which  the  liquid  will  .. die.  I 
is  required  If  the  walls  of  the  cylinder  arc  suitable  tliev  .  m 
provide  the  needed  suit  ace 

?  If  the  surlace  is  nol  available,  that  is.  not  cite,  'u. 
as  an  active  site  lot  the  fiuination  ol  the  new  phase.  ;n,r 
the  new  phase  will  have  to  begin  homogeneously.  lh..i  ;s 
ll  will  have  to  form  upon  itself  Classical  tlieiniocly  naro. 
tells  us  thal  ill  otdei  to  have  a  dioplet  ol  linite  i.idm-  m 
equilihimm  with  the  gas.  the  gas  must  he  supeis.iiui.iU. I 
In  I  act.  the  smallest  possible  udiux  of  swell  .1  sjWiiifx.fi 
drop  is  invvisel>  proportional  to  the  natui.il  lop  ol  tie 
saturation  ratio  Once  the  dioplet  exceeds  this  Md.ns  :r 
will  glow  without  hunt  and  ptodue v  sut  lie  lent  x  oieUms.-d 
phase  to  establish  a  stable  thctmodyieinue  ewuilihi uim 

V  I  he  reason  the  eas  supei salinates  is  because  it  0 
cooled  at  a  finite  tooling  tale  as  opposed  to  ,m  inhndx  K 
slow  eoolmp  rate  lopmcd  lot  a  ijuaxistatn  ffex  1  \\  >/< 
classic  a  I  thei  mods  nainiv  x  I  ven  when  a  pas  is  stipt  1x.1i  11  j  a 
ted  the  new  phase  in. in  not  loim  at  ait  appreciable  oilr 
beeause  Ihete  is  .1  positive  amount  ot  x-ncic\  mpmed  \" 
torm  the  droplet  ol  me  minimum  radius  t»  I  xxeuhalK  the 
reason  foi  tins  is  due  to  the  enetes  ot  tomuiini)  ot  tlu 

siulaee  4nr*’<r  where  u  is  (he  mu  hue  tension  I  Inis  tlu 

higher  the  xupeixadnafiou  of  (he  pt(s  the  Mji.illrr  r.  1  » 

the  Milallei  is  tin-  eueipy  t'amei  I  * »  tin  t  ■  «i  1n  .1l  e  *n  ot  lit 

diop'el  <d  xiilual  si/i  .  and  Ilu  nmn  r.ipnj  t.  tin  1  1 1 
bn  illation  ot  (tie  ilmplets  of  x  lush  is  I  Inis  am  mu  lied  1 1 1 
fuithei  out  '*t  equtiibr iiim  the  thmnody  iuiuk  si  •  ; 
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TABLE  1 

Cluster  si/e  for  bulk  density 

No.  ot  interatomic 

2  a 

do  q„ 

(2n*  l)  ' 

Distances,  n 

a 

3 

6 

343 

4 

8 

729 

b 

10 

1331 

tin'  greater  is  the  rale  process  which  si  rives  the  system  hack 
low  .ml  ei|u  ilihr mm 

I  or  v  nr  u  a  1 1  y  every  type  ol  actual  phase  change, 
whether  lormed  homogeneously  as  pure  droplets  or  hetero¬ 
geneously  as  a  droplet  growing  on  a  surface  or  a  foreign 
site,  whether  a  phase  change  between  gas  and  liquid,  gas 
and  solid,  or  solid  and  liquid,  the  starting  or  critical  size 
involves  a  cluster  or  aggregate  whose  si/.e  g  is  less  than 
1000  molecules  (atoms  if  monatomic  species)  per  cluster 
Beyond  die  critical  si/e  the  new  phase  can  grow 
without  limit.  How  docs  this  si/.e  range  compare  with  the 
si/e  limit  at  which  hulk  properties  are  valid?  In  order  to 
estimate  this  si/e  let  us  consider  pr opertie.s  at  an  interlace 
between  two  phases.  For  example,  consider  the  density  as 
a  function  of  position  across  an  interface  as  shown  in  big. 

I  typically  the  density  change  from  gas  to  condensed 
phase  is  ol  order  10001.  Mathematically  we  might  con¬ 
sider  that  the  interface  occurs  at  a  plane  of  zero  thickness. 
Miciuscopically .  when  one  considers  that  a  phase  is  made 
up  of  discrete  atoms  or  molecules,  the  density  change  will 
occur  across  some  number  n  of  mean  intermoleeular 
( mleratom ic I  distances  a  Ihus  the  characteristic  length 
over  which  the  density  goes  from  hulk  liquid  to  bulk  gas  is 
\  na  where  n  would  be  no  smaller  than  .1  ami  perhaps  as 
high  as  10  or  more.  Now  consider  a  slab  of  condensed 
phase  with  gas  phase  on  two  sides  as  shown  in  Fig  IB. 

I  he  minimum  characteristic  length  at  which  bulk  properties 
would  obtain  at  the  center  of  the  condensed  phase  would 
be  roughly  2  A  II  we  consider  a  three-dimensional  con¬ 
densed  phase  tie.,  a  droplet-like  entity)  we  can  estimate 
the  numhci  of  atoms  or  molecules  required  m  order  to  have 
bulk  properties  even  though  only  al  Ihe  center  ol  the 
droplet  Ihe  results  of  this  estimate  are  shown  in  Table  I 
where  Ihe  number  of  interatomic  or  intermoleeular 
distances  required  for  the  density  change  from  that  of  gas 
to  condensed  phase  is  varied  from  3  to  5.  We  see  that  a 
cluster  must  have  more  than  343  (1,331)  molecules  or 
atoms  in  the  cluster,  if  n  =  3(5).  before  there  is  enough 
material  in  the  droplet  so  that  there  is  bulk  property  any¬ 
where  in  Ihe  cluster.  Note  that  if  n  =  5  then  virtually  all 
the  critical  si/e  clusters  for  the  formation  or  nucleation  of 
a  new  phase  tall  within  the  region  wherein  hulk  physical 
properties  will  not  be  present  anywhere  in  the  droplet 

t  he  above  estimate  was  made  with  density  as  the 
property  that  was  varied  there  are  many  physical  proper¬ 
ties  other  than  density  that  may  be  considered  and  numer¬ 
ous  estimates  have  hcen  made  in  the  literature  of  the  cluster 
si/e  at  which  these  other  properties  begin  to  deviate  hum 


those  of  the  bulk  these  estimates  aie  based  on  appioxi 
male  models  ami  are  llicieloic  subject  to  .  halh-ngc  I  >m 
example,  it  lias  been  eslimaled  dial  only  '  lo  1  al**nis  aie 
neeessaiy  lo  appioximalc  Ihe  elec  limn.  slimluo-  ol  .. 
m el. 1 1  e lux! cl  whereas  otlici  estimates  fui  piopcili.  vi.  !j  i 
Ihe  latent  heal  oi  surtax  e  tension  requite  sevei.il  I li.uixui.d 
alums  or  molecules  in  a  clustcl  b.-toie  bulk  value  -  ao. 
appioached 

Clusters  are  everywhere 

From  the  foregoing  discussion  il  is  evident  Ih.u 
clusters  are  important  wherever  there  is  a  physix.il  phase 
change  involved  l  lns  one  phenomenon  alone  xan  involve  a 
tremendous  sx'ope  ol  applications,  for  exjiiipb-  tin  nmlea 
turn  and  growth  of  small  metal  particles  on  subxliatcs  i 
very  important  tor  thill  tilm  and  solid-state  devices  Manx 
atmospheric  phase  changes  such  as  ihe  formation  ot  solid 
aerosols,  water  condensation,  and  hail  lotnialion  involve 

small  nuclei  to  initiate  the  processes  Hiysix.il  propciiu . . 

metals  are  determined  ill  pari  by  file  xluslei  .  mu  cut  nil  ion 
and  therefore  the  metal  gram  si/e  Many  chemical  engincci 
ing  processes  require  the  condensation  ol  cry  stables  o' 
formation  of  liquid  droplets  that  miliate  Ihe  loimalion  oi 
a  new  ptiase.  the  products  ot  eombustixm  such  as  walci 
or  carbon  dioxide  ami  steam  hi  many  types  ol  engines  am! 
turbines  can  condense  into  dioplets  and  cause  subsequent 
damage  to  paits  downstream  ol  the  condensation  pi-wcss 
Performance  of  boileis  aiixl  comlenseis  in  large  eip.iipmeu' 
depend  on  um  lealion  phenomena 

(loud  chambers  and  bubble  chambers  lo:  use  m 
high-energy  physics  are  examples  ot  ladiatioiiumlm  e  ! 
nuelealion  Numerous  geologic  ami  cosmologic  plus-, 
changes  are  import  ml  within  the  earlli  ami  bey  oml  Ihe 
hsl  can  go  oil  and  on.  One  very  important  appln.ition  toi 
small  clusters  is  in  (lie  field  ol  catalysis  wlieie  1  ■ »:  xoim 
chemical  reactions  it  has  been  found  that  small  metal  ami 
metal-oxide  clusters  are  much  more  vltcclive  lh.ni  bulk 
material  in  promoting  chemical  reactions  Oilier  divers, 
areas,  such  as  the  study  ol  the  theory  of  liquids,  involve  a 
model  which  is  essentially  that  ol  a  cluster  Ihe  use  ui 
heterogeneous  cloud  seeding  agents  m  the  atmosphere  toi 
Ihe  control  of  rainfall,  increasing  snow  pack  m  luonntaiu- 
n n cl  the  suppression  ol  hail  damage  to  ciops  ale  -ilu 
examples  ot  (lie  use  of  small  clusters 

Clusters  are  systems  which  are  small  m  thicc  dimer, 
sions  but  there  are  examples  ol  sy  stems  Dial  aie  small  m 
one  and  two  dimensions  which  also  exhibit  behavioj  mu 
typical  of  bulk  material.  For  example,  thin  hints  aie  small 
in  only  one  dimension  and  once  they  become  tlunnci  than 
about  lOOOAtlOO  nml  their  properties  begin  to  sinus 
deviation.  I  xamples  ot  systems  which  are  small  m  two 
dimensions  are  those  ol  libers  or  whiskers  their  piopeilies 
also  are  different  from  bulk  In  particular  lot  crystalline 
materials,  due  to  the  small  volume  involved.  Ihe  piobabilip, 
for  imperfections  or  dislocations  is  reduced  thereby  texuli 
mg  in  materials  with  unusually  high  strength  Is  it  possible 
that  some  materials  in  small  cluster  lornt  might  be  ideal  loi 
superconducting  mateiials'’  1  Ins  prosper  t  would  have  gic.it 
economic  implic  alums 

Examples  of  properly  changes  of  small  aggregates 

Ihe  problem  ol  coinpansoii  between  a  eluslei  and  a 
bulk  phase  of  the  same  material  may  be  considered  m 
terms  ol  a  two  step  change  I  tie  lust  step  would  involve 
identilyuig  a  ceitain  numbei  (gl  ol  uiolec  ules  m  the  bid) 


'iiis 


T  Ht  PI  I YS  ICS  I  I  Anil  u 


NOVI  MHt  K  1 '» /■» 


l-ig.  2.  I hc-Mclical  stay  diffraction  patterns  arc  shown  lor  a  b.Hiy-contsTes!  culm  latuv  .n 
a  I  unci  ion  of  lire  number  of  unit  cells  in  the  cluster.  I  Ire  y  court!  inure  is  a  relative  scale  l  a 
peak  height.  while  for  \.  1  is  a  function  of  the  angle  of  scattering  tor  the  \-ra%  I  hits  r It-, 
chances  in  tl  iff  ruction  pattern  tiecur  even  lor  perfect  unit  cells  as  the  cluster  si/c  acts  small,  i 
with  lirurnatrc  chances  below  ’0  cells  per  cluster,  d  rum  t'.  M->r-  /unri  anil  II  l  Km 


lig.  I .  t  he  minimum  energy  configuration  for  van  tier  Wauls  type 
aloimc  clnslers  with  g  <55  atoms  is  that  til  the  icosahetlral  packing. 
A  -ornplcte  I  1-atom  icosahedial  is  shown  here. 


phase  and  imagine  them  removed  of  isolated  Itotn  the  hulk 
with  exactly  the  same  structure  tie  .  atomic  positions! 
The  I'initeness  of  the  new  system  alone  will  cause  a  change 
in  some  physical  measurements,  tor  example  diffraction 
patterns.  The  second  stage  of  tile  process  then  would  he  to 
allow  the  cluster  to  relax  into  its  equilibrium  position  as  an 
isolated  entity,  lilts  then  would  resuli  m  further  tJiulion- 
of  some  physical  properties. 

Consider  the  ease  of  \-rj>  or  electron  Jtllr.it  //..■- 
from  a  crystal  which  is  perfect  but  of  Unite  extern  (  <>m 
plete  constructive  interference  at  the  Bragg  angles  and 
destructive  interference  at  the  angles  m-hetween.  resuli  mg 
in  infinitely  narrow  diflraetron  peaks,  require  a  crystal  ol 
infinite  extent.  As  the  crystal  becomes  tmite  m  si/c  the 
peaks  are  broadened  and  peak  heights  reduced,  but  with 
little  additional  changes  as  the  st/e  goo*  Horn  inlimly  down 
to  micron  diameter  lire  broadening  oecuis  lor  linile 
crystals  because  at  angles  slightly  » . 1 1  the  Biagg  angle 
(say  A(f )  there  are  not  enough  scalterers  in  the  uysi.il 
lattice  for  complete  destructive  interlerenee  I  he  lew-  i 
atoms  (molecules)  in  the  lattice  the  poorer  is  the  dcstuK 
live  interference  for  a  given  AO  Hence  the  peaks  are  in¬ 
creasingly  broadened  as  the  si/e  decreases  1 1  his  phenom 
enon  is  analogous  to  the  resolving  power  of  a  diffraction 
grating  like  those  used  in  spectrographs  and  monoehrom- 
eters.  Hie  greater  the  extent  of  the  grating  flu-  belief  is  its 
resolution  I  Below  this  si/e  I  Ip  in  diameter)  the  peaks  begin 
lo  ..Iter  more  tapully  1  and.  as  the  mimbei  ol  unit  cells 
deetease  below  about  1011.  the  patients  take  on  substantial 
changes  as  seen  in  1  ig  J  I  bus  even  wdh  pellet  l  mill  sells 
(body-centered  cubic-IH  t  ,  in  this  easel  which  have  n- 
interatomic  variations  m  distance,  no  thcimal  vibiulion  -•! 
the  atoms  about  cspnlibr  iiiin  positions  and  mi  sty -.id 
delects,  the  dittraction  pattern  is  slrastically  dillt-icnl  liom 
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ill, ii  nl  ihe  hulk  for  the  small  clusters.  In  order  to  model 
successively  more  realistic  clusters  tor  comparison  to  ex¬ 
periment.  one  must  add  to  the  perfect  model  the  above 
mentioned  el  fee  Is. 

Argon  clusters 

I  el  us  consider  the  structure  and  properties  of  a 
"simple"  monatomic  species  which  exhibits  little  or  no 
.licmicnl  bonding  but  only  the  much  weaker  van  der  Wauls 
alli.ielion  Argon  is  the  most  common  of  this  family  arid 
us  bulk  crystalline  structure  is  face-centered  cubic  (I  CC. 
i  e  .  a  cube  with  atoms  at  the  corners  and  centered  on  the 
sis  laces!  I  he  interaction  potential  between  two  argon 
atoms,  dial  is.  the  attractive  or  repulsive  force  which  is  a 
I  unci  ion  ol  ihe  distance  between  them,  operates  between 
any  two  argon  atoms  in  a  simple  additive  manner  losing  a 
simple  potential  function,  which  fitscxperimcnt.il  measure¬ 
ments.  minimum  energy  configurations  for  clusters  as  a 
t tint,  t ion  ol  si/e  have  been  computed  these  so-called  static 
.  a Ic ul.i I  ions  have  been  made  by  numerous  researchers  ll 
was  found'  that  for  clusters  up  through  size  55  the  equilth 
i  min  coul  iguialion  is  not  ICC  bul  lends  toward  the  ico- 
s.dicdral  packing  shown  in  lig  i.  With  single  atoms  being 
added  sequentially  the  icosabedral  structure  cannot  propa¬ 
gate  dselt  and  at  some  point  the  minimum  energy  con- 
I puliation  lor  (he  growing  cluster  will  revert  to  face- 
.  entered  ,  ubic 

I  hcio  is  experimental  evidence  tor  the  postulated 
iios.ihcdi.il  siiucluie  lor  ciy stalling  argon  obtained  using  a 
nozzle  mold  ular  beam  technique  hy  a  research  group  at 
the  I  'liivcisity  id  Paris  (Orsjy ). '  ■4  A  typical  molecular 
beam  .  onl  iguialion  is  shown  in  big.  4  and  has  been  used 
over  die  Iasi  ’0  years  for  a  great  variety  ol  research  invest  i- 
g.ilnms  ‘  ll  eonsisls  ol  a  gas  source  winch  is  cillici  a  small 


nozzle  oi  an  oriluc  tree  jet  li  e  .  a  pm  hole  m  a  linn  metal 
disci  and  is  operated  at  pressures  high  enough  so  dial  flic- 
expansion  through  the  hole  forms  a  tree  jet  or  plume  in  the 
evacuated  first  stage  of  the  beam  As  the  gas  expands  d 
cools  adiabatieally.  becomes  supersaturated,  anil  clusters 
begin  to  nucleate  rite  pressure  in  the  first  stage  of  pumping 
is  low  enough  so  that  the  density  in  the  free  jet  becomes  so 
low  that  the  expansion  becomes  essentially  collision  less 
Somewhere  near  this  point  in  the  expansion  an  aperture  is 
used  to  define  a  beam  ll  is  formed  as  a  cone  willi  a  hole  in 
it  to  minimize  aerodynamic  mterfeiencc  with  the  beam  that 
passes  through  ll  and  is  called  a  skimmer  I  icquenlly  i 
second  stage  ol  pumping  is  used  in  oi.lci  to  key  ih. 
pressure  m  the  scattering  chamhei  as  low  as  possible  I  > 
a  given  gas.  it  the  gas  mlcl  piessuic  po  is  low.  clusleimg  i- 
progressively  increased,  more  and  mole  clusleis  will  lorn: 
resulting  m  a  cluslet  beam  hot  siruclurc  analy  sis  ol  lie 
cluster  beam,  a  high-energy  electron  gun  is  placed  al  ng hi 
angles  to  the  cluster  beam  just  downstream  ol  the  .ollim 
ater  in  this  scattcung  chamber  Ditfiactron  paltcins  .tic 
taken,  usually  on  film .  with  a  beam  stop  to  .auh  the 
primary  undittractcd  electron  beam  I  lie  cludci  Ivan 
intensity  is  monitored  using  an  ionization  gauge  oi  i  ma  s 
spectrometer 

I  Ice tron  dilliaction  patterns  trom  clustered  argon 
beams  have  hcen  obtained  and  compared  lo  theoiclic  al 
siruclurc*  calculations  I  he  patterns  lot  huge  .Itisicis 
(g  '  5001  compare  well  to  lit  bulk  slim  lure  While  the 
size  ranges  liom  g  '  50  the  data  is  iepi  .uluced  well  using  a 

model  which  is  a  combination  ol  icosuhcdralclustei'S  and  be 
gas.  4  l  ire  tiactmn  which  is  lcosalicdial  goes  liom  Ml  to 
O'!  as  g  goes  liom  50  lo  ’0  I  ven  Ibougli  one  may  esp.  i  1  a 
small  range  t  *  I  sizes  m  Ihe  «  In  .lei  beam  I  lie  .  omp  n  t  .<  m  i  . 
nevertheless  very  good  arid  thus  pomades  Ihe  eypeittu,  nl  il 
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1  ig.  5. 1  lu  oiiMn  .il  molecular  dynamics  calculations  tor  argon  dusters 
t.m  be  presented  on  a  temperature  versus  energy  per  atom  plot  as  a 
turn  non  .*t  cluster  size.  I'he  so-called  van  der  Waals  loop  is  shown 
with  I  the  melting  point  temperature  and  Um  the  energy  of 
sapoi izatixii  I'he  results  shots  evidence  ot  a  van  der  Waals  loop  and 
thus  a  phuv  change  even  tor  the  size  4  oi  tetramcr 


Individual  snapshot  pictures  trout  1 3-atom  and  4  Tatum 
dusters  were  obtained  from  Monte  Carlo  simulation  at  various 
icmpcr.untcs  Note  the  leosahcdial  packing  at  very  low  temperature 


verification  for  the  appearance  of  these  minimum  o»org> 
icosahedral  argon  structures  under  proper  inohv wlar  beam 
operating  conditions. 

In  addition  to  static  structure  configuiat ions.  d>n.im 
ic  calculations  have  been  earned  out  by  a  tumihei  ot 
researchers  in  which  there  is  kmcln  atul  potential  eneigv 
interaction  within  the  cluster  Basically  two  iheoieinal 
techniques  have  been  employed  I'he  hrsl  is  vailed  mole,  >< 
iar  dynamics  in  which  each  particle  is  given  an  initial 
position  and  velocity  and  the  evolution  ot  the  dustei  is 
followed  with  ordinary  Newtonian  mechanics  I  he  second 
technique  is  called  Monte  Carlo  with  the  ejlcul.it i« ok 
initiated  using  an  ensemble  of  g  atoms  located  within  some 
arbitrary’  spherical  volume  Using  a  given  iiiteratorm* 
potential,  the  potential  eneigy  of  the  initial  eonfigutati.  o 
is  calculated.  I' hen  one  of  the  atoms  is  chosen  at  random 
and  given  a  random  displacement  A  new  potential  enetg, 
is  then  calculated  followed  by  a  second  random  ilispl.it  . 
men!  After  a  large  number  o!  these  steps  the  potenliii 
energy  is  ensemble  averaged  and  a  cluster  equilibrium  'date 
determined  All  of  the  tbermod  v  njimc  properties  i  an  thcii 
he  calculated  from  these  or  subsequent  configurations 

Results  obtained  using  molecular  dynamics'  show  a 
triple -valued  van  der  Waals-typc  loop  m  the  eneigy  |h  ’ 
atom  as  a  lunction  ot  temperature  as  the  cluster  size 
increased.  The  structuie  of  the  cluster  passes  from  solid  t- 
liquid  in  the  loop  region  as  temperature  is  uk  teased 
analogous  to  a  van  der  Waals  loop  lor  an  isotbeim  m  .i 
pressure-volume  diagram  for  a  bulk  phase  change  im. 
l  ig.  M 

Several  mterestmg  features  ate  piedicted  ' 

1  the  melting  point  of  a  cluster  is  dependent  <m  in 
size,  that  is.  the  smaller  the  v  luster  fhe  lower  the  me/fu.i* 


and  ev idem  e  of  a  transition  to 
a  liquid  like  or  amorphous 
Miueture  in  the  temperature 
lange  ot  4u  btt  K  in  agreement 
with  the  molecular  dynamics 
.  ah  (Rations  °  1  Reprinted  w  ith 
pemiisslon  Horn  I  I  Abra¬ 
ham  I 
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irmju't.iiiuv  I- von  lor  iho  100-atom  olusioi.  a  inciting  point 
tomperutute  of  about  40  K  is  mote  than  a  factor  of  1  lower 
Ilian  iho  S4  k  moiling  point  temperature  of  hulk  solid 


.11  gun 

'  I  lio  I  *  atom  cluster,  whic  h  is  a  c  omplete  loosafio- 
dian  has  a  "heal  of  fusion”  {ll|tl  ill  lig  51  which  is  larger 
than  the  dusters  above  ami  below  it.  I  Ins  again  points  to 
t lie  minimum  m  potential  energy  for  this  configuration  or 
an  optimum  m  van  dor  Wauls  bonding,  thus  tequinng  more 
eueipy  pel  atom  tor  Us  phase  change 

'  I  or  the  larger  si/e  elusleis.  as  then  lempeiature 
increases  toward  the  melting  point,  other  calculations 
pied ic t  that  the  duster  is  partially  liquid  and  partially  solid. 

1 1  has  bee  n  postulated  elsewhere  that  elusteis  m  this  regime 
would  be  liquid  like  near  the  snilaee  with  the  cote  ic  main-  T 
me  solid  Q 

-\n  example  ol  the  results  that  can  be  obtained  using 
Monte  t'ailo  calculations  art*  shown  m  I  ig  t>  Here  single 
\  onheiii at  ion  snapshots  are  shown  lot  a  I  s  atom  and  a  4  4 
atom  duster  at  temperatures  between  0.0 1  K  and  I  0(1  k 
We  see  that  as  the  temperature  nears  0  k  the  structure  is 
vety  dose  to  that  of  the  icosahodial  minimum  energy 
wuihgmation  It  is  also  noted  that  as  the  lempeiature 
uu  reuses  to  the  range  of  4()*M)  K  the  structure  is  detmitely 
mmiiii:  away  liom  icosahedral  m  qualitative  agreement 
with  the  molecular  dy  mimics calcul.it ions. 


X  /  D  «  39.2 


P  »  3  bar 
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Wilier  elusleis 

Hit-  above  see  lion  dealt  with  simple  monatomic 
■>\  'Iv  ins  possessed  ol  weak  interaction  vine  lo  the  van  .lei 
Wauls  .nn.ielive  force  Here  we  will  consider  .1  e.ise  ol 
hvdiogen  bonding  wlueli  is  imieli  sponger  than  van  .let 
W  .nils  al  ii.it  Inin  by  a  factor  of  about  '() 

\n  investigation  ol  l he  structure  ol  water  clusters  in 
Hie  pas  phase  using  high  energy  electron  dillraction  was 
un. lei  taken  m  0111  laboratory5  a  lew  years  ago  using  lire 
niolevul.ii  beam  apparatus  shown  in  big.  4.  I  lie  structure  of 
these  1  Insieis  was  mil  the  usual  hexagonal  lorm  ol  ice  but 
111sic.nl  was  .hatnonJ  t-uhic.  1  lie  diffraction  patients  are 
shown  in  I  ig  ’  with  stronger  paMerns  occurring  al  higher 


1  J _ 1 _ L 

12  3  4  5 

S  =  (4 II  /X)  sin  (6/  2) 


I  ig.  7.  I  lev  l  roll  diffraction  put  terns  for  water  with  the  background  subtracted  out 
uie  shown  in  A  as  obtained  Ironi  the  molecular  beam  apparatus  shown  in  1  ig  4 
lot  three  beam  starting  pressures.  The  average  si/e  for  these  clusters  is  in  the  ranee 
of  y  'nOO  to  .130  and  is  characteristic  of  the  bulk  diamond  cubic  structure  ot 
ice  I  he  diamond  cubic  structure  is  shown 
111  It  with  the  water  molecules  shown  only 
.is  v  mles  to  simplify  the  presentation  which 
is  basically  that  of  l(('  with  4  molecules 
micnoi  to  the  unit  cell  located  as  shown 
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llu-  utlililion  of  gas-phase  molecules  condensing.  Just  as 
with  argon,  as  the  size  gets  larger  the  minimum  energy 
configuration  reverts  to  one  of  the  more  usual  hulk  forms 

Metal  clusters 

l  or  a  number  of  years  in  our  laboratory  we  have 
earned  out  a  systematic  study  of  the  structure  of  small 
clusters  formed  in  the  vapor  phase.  Due  to  the  relatively 
low  vapor  pressure  of  metals,  plus  their  relatively  strong 
bonding  (a  factor  of  about  20  greater  than  the  hydrogen 
bonding  m  water),  and  high  surface  tension,  the  design 
criteria  for  a  cluster  molecular  beam  is  much  different  from 
Ilia!  for  gaseous  species.1  1  The  most  important  design 
wilerinn  for  a  metal  cluster  beam  source  which  will  pro¬ 
duce  small  si/e  clusters  with  sufficient  cluster  concentration 
to  conduct  experiments,  is  that  the  supersaturation  be  very 
high  Ibis  produces  an  enormous  cluster  nucleation  or 
lornijlion  rate  so  that,  on  the  time  scale  of  the  experi¬ 
ment.  cluster  concentrations  of  It)10  to  It)1 ''/cm'  are 
produced  With  that  many  clusters  competing  for  the  re¬ 
in. lining  metal  vapor,  the  average  size  remains  quite  small. 

\  metal  cluster  oven  design  used  in  our  laboratory  is 
shown  in  I  ig  *•  I  he  metal  is  vaporized  in  an  inert  atmo- 
xphore  pressure  at  about  I  I'orr  (1/760  of  an  atmosphere) 
and  the  gas  flowing  slowly  over  either  a  vapor  heater 
lilanicnt  oi  crucible  toward  the  first  orifice  N,  I’art  of  the 
sample  passing  through  N|  is  pumped  by  mechanical  pumps 
and  the  lemamder  passes  through  a  supersonic  orifice  Ns 
and  •becomes  the  cluster  beam  to  be  examined  using  high 
energy  electron  diffraction  Ibis  source  replaces  the  gas 
molecular  beam  in  I'ig.  4  and  is  connected  directly  into  the 
xcalteung  chamber  A  number  of  metals  have  been  studied 
using  this  soutcc.  resulting  in  some  surprising  and  exciting 
leslllts 

Consider  the  case  of  indium.  Its  hulk  structure  is 
lace  centered  tetragonal  (that  is,  face  centered  but  with  a 
height  ,  which  is  greater  than  its  square  base  a),  f  igure  10 
indicates  that  when  the  indium  clusters  have  5  000  atoms 
pel  i  Insti  l  oi  larger  they  exhibit  the  hulk  structures  just 
mentioned  However,  as  they  are  made  progressively 
smaller,  we  carl  see  that  the  height  dimension  c  decreases 
while  the  size  ot  the  base  a  increases  with  the  result  that  at 
about  g  2000.  c  =  a  and  the  structure  becomes  face- 
centered  cubic  (ICO  Ihe  structure  appears  to  be  more 
compact  (Are  the  surface  tension  forces  which  create 
higher  pressures  in  small  droplets  (lie  motivating  force 
here'1)  However  an  examination  of  the  density  of  the  small 
size  l  td  structure  shows  that  it  is  within  I  .  of  the  bulk 
density  Nevertheless  we  feel  that  the  structural  rearrange¬ 
ment  is  sign  il  lean  I  and  occurs  at  a  cluster  size  in  a  range  of 
’  000  to  4  000  atoms  per  cluster  which  can  be  considered 
as  iclai ively  large  ti  e  containing  500  to  I  000  unit  cells). 

As  a  second  example  consider  the  case  of  lead, 
figure  IIIK  shows  the  variation  of  the  unit-cell  parameter 
twlinli  is  ICC)  with  size  Note  the  very  small  variation  in 
unit  i  ell  dimension  with  size  (obseive  lliat  the  y  coordinate 
is  pie. UK  expanded  ovei  that  ol  pail  A)  so  lli.il  on  lust 
examination  the  data  (oi  lead  does  not  appear  interesting 
(  aielul  exam  illation  ol  the  data  shows  that  it  is  consist 
cut  with  a  model  ol  a  cluster  which  lias  a  core  which  is 
cry slallized.  plus  a  surface  layer  which  becomes  progres¬ 
sively  more  dixonlered  or  liquid-life  as  the  cluster  gets 
smaller  1 


fig.  10  Hie  unit  coll  lattice  parameters  lie 42 hi  c.  base1  a.  ami  then 
ratio,  for  indium  microclusters  are  plotted  as  a  function  ot  the 
number  atoms  in  the  cluster.  g.  and  the  average  cluster  diameter.  1> 
Broken  lines  indicate  the  corresponding  values  in  bulk  crystalline 
materials  at  various  temperatures.  Typical  experimental  eir»>r  hao 
are  shown  by  vertical  linos.  Deviation  from  bulk  structure  hct-mN.o 
the  diameter  approaches  5000  atoms  per  cluster  average  si/o  In  b 
Ihe  unit  cell  parameter  for  lead  is  plotted  as  a  function  id  g  Vm 
the  expanded  scale  for  unit  cell  parameter  with  very  little  chance  in 
the  lead  unit  cell  with  size.  I  he  indium,  on  the  other  hand,  lias 
undergone  a  structural  change  from  face-centered  tetragonal  to  lace 
centered  cubic  as  the  si/o  is  decreased  to  the  2  nOO  atom  pci  dustoi 
range. 


Because  the  atomic  weight  of  lead  is  so  large.  Hie 
possibility  of  multiple  electron  scattering  within  the  cry  dal 
lite  is  of  some  concern  Iheiefore  another  I  VC  monatomic 
structure,  that  ot  argon  mentioned  above  (g  '  MUM.  Ii.h 
been  examined  ami  the  same  cited  seen  I  bus  with  do 
creasing  si/e  a  larger  fraction  of  the  cluslet  becomes  stir  hue 
with  Ihe  crystalline  center  becoming  progressively  .1  sm.ilb  1 
fraction  of  the  atoms  in  the  cluster  Piesumahly .  then,  the 
structure  of  the  clusters  in  the  100  to  I  000  atom  pi  t 
cluster  range  would  be  completely  amorphous  01  1 1« |u h I 
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like  .ts  suggested  by  the  argon  data  reteried  to  earlier  for 
g  50  Hut  then  again  we  may  he  surprised  by  the  small 
elustei  structure  as  with  the  previous  examples. 

Cluster  studies  with  the  laser 

With  the  advent  of  tunable,  highly  monochromatic, 
high  powei,  and  very  short-pulsed  lasers  many  traditional 
and  new  reseaich  areas  are  blossoming  forth  l'he  study  of 
small  clusters  is  one  ol  these,  and  recent  advances  have 
occulted  at  a  very  rapid  rate  Some  beautiful  spectroscopic 
studies  have  been  carried  out  by  scientists  at  the  University 
ol  Chicago  1  they  employ  highly  cooled,  free -jet  expan¬ 
sions  of  helium  containing  a  very  small  quantity  of  the 
species  ol  interest  (such  as  iodine)  and  have  obtained  highly 
lesolved  speetroscopie  data  on  vibrational  and  rotational 
slimline  ot  very  small  clusters  Impressive  progress  on  the 
study  ol  metal  clusters  has  been  made  using  matrix  isola 
lion  techniques  (inel.il  clusters  condensed  m  a  matrix  of 
solid  aigon  for  example)  with  laser  spectroscopy,  particu¬ 
larly  at  the  University  of  lotonto  14 

Very  recently  laser  spectroscopic  studies  of  small 
metal  clusters  m  the  gas  phase  have  been  undertaken  m  a 
nearly  ideal  experimental  configuration  at  the  University  in 
Hern.  Switzerland.1  5  Hus  technique  involves  the  formation 
ol  a  metal  cluster  beam  and  selective  electronic  excitation 
by  cluster  size  using  a  tunable  laser  A  second  high-intensity 
laser  pulse  then  strikes  the  sample  lom/uig  only  those 
species  which  have  been  electronically  excited  file  entire 
expel  iment  is  performed  in  the  neighborhood  of  the  inlet 
optics  to  a  quadrapolc  mass  spectrometer  This  extra- 
oidmaiy  technique  provides  spectroscopic  data  for  clusters 
ol  a  single  size  In  a  sense,  it  represents  the  ultimate  cluster 
expel  iment 

I  he  additional  ability  to  use  sub-nanosecond  and 
picosecond  (that  is  less  than  10  4  to  10  1  2  sec)  laser  pulses 
ot  known  and  of  variable  frequency  is  being  used  in  the 
study  ol  the  dynamics  of  many  systems,  including  cluster 
systems  Hus  very  powerful  technique  wdl  provide  the 
ability  to  study  the  dynamic,  stalc-to-state  changes  in 
clusters 

Summary  and  prospects  for  the  near  future 

flic  existence  of  matter  in  small  aggregates  or  clusters 
is  seen  lo  arise  in  a  great  variety  of  important  scientific 
studies  and  engineering  applications.  With  the  advent  of 
new  quantum  theoretical  methods  to  calculate  the  physi¬ 
cal  properties  of  systems  with  up  lo  .10  atoms,  the  extreme¬ 
ly  rapid  advance  in  laser  physics  and  optics,  the  recent 
me  oi  pot  niton  of  nucro-processcrs  into  multi-channel 


detection  instrumentation  with  the  capability  loi  single 
photon  or  single  particle  pulse  counting,  and  the  pm 
gressively  easier  and  less  expensive  access  to  compute-!-  lot 
both  theoretical  and  experimental  uses,  a  near  explosion  ol 
activity  in  this  fascinating  and  important  area  ot  small 
cluster  research  has  begun 
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ABSTRACT 

Silver  is  evaporated  into  an  inert  gas  flowing  at  pressures 
between  0.1  to  3  torr  and  nucleates  to  form  clusters  in  a  size 
range  of  40  to  110  X.  This  two-phase  mixture  flows  through  a 
two  stage  molecular  beam  into  a  modified  electron  microscope. 
Electron  diffraction  patterns  are  taken  in  the  cluster  beam 
and  analyzed  using  the  bulk,  face  center  cubic  unit  cell  in 
conjunction  with  a  twelve  beam  multiple  scattering  program. 
Deviations  from  the  bulk  structure  are  seen  and  are  larger 
for  smaller  cluster  size.  Several  possible  interpretations  are 
discussed.  The  effect  of  different  carrier  gases  on  oven  per¬ 
formance  and  cluster  size  is  also  presented. 
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INTRODUCTION 


The  objective  of  this  study  is  to  obtain  information  on 
the  structure  of  small  clusters  of  metal  atoms,  prepared  by 
evaporation  into  a  low  pressure  inert  gas,  by  analysis  of  their 
electron  diffraction  patterns.  Potentially  large  effects  of 
supporting  substrates  upon  the  clusters'  structures  were  avoided 
by  employing  the  quenching  gas  as  a  carrier  to  transport  the 
clusters  through  a  differentially  pumped  gasdynamic  system 
(forming  a  "molecular"  beam)  into  a  high  vacuum  region  where 
the  electron  diffraction  pattern  can  be  recorded^.  A  neces¬ 
sary  subsidiary  goal  is  to  discover  conditions  which  produce 
sufficiently  dense  metal  cluster  beams  so  that  electron  dif¬ 
fraction  patterns  can  be  observed 

EXPERIMENTAL 

The  metal  cluster  generator  has  been  described  previous¬ 
ly^  and  is  schematically  illustrated  in  Fig.  1.  This  oven 
cluster  source  has  been  mounted  into  one  part  of  the  "diffrac¬ 
tion  chamber"  of  an  Hitachi  HU-llA/B  which  is  operated  as  an 
electron  beam  source  (condenser  lenses  only  used) .  The  elec¬ 
tron  microscope  has  been  modified  with  the  addition  of  a  coarse 
aperture  just  above  the  cluster  source  and  additional  pumping 
capacity  from  a  15  cm  diffusion  pump  with  1,000  /sec.  It  is 
connected  to  the  diffraction  chamber  directly  opposite  the 


oven.  Electron  diffraction  patterns  are  recorded  on  glass 
photographic  plates. 

The  metal  vapor  source  used  was  a  small  helix,  hand-wound 
tungsten  wire  (0.38  mm  diameter)  into  which  one  or  more  small 
ingots,  formed  from  99.999  %  silver,  are  inserted.  The  amount 
of  silver  in  each  charge  is  estimated  to  be  in  the  range  0.1  - 
0.25  g.  The  gas  inlet  tubing  was  shortened  so  that  it  now 
discharges  through  the  heater  rather  than  between  heater  and 
as  shown  in  Fig.  1.  Pressures  measured  at  various  points, 
with  gas  flow  on  but  power  oven  off,  are  recorded  in  Table  I 
together  with  the  pressure  ratios  across  the  copper  inlet 
tube  and  across  nozzles  and  Except  for  the  lowest 

pressure  with  helium,  which  is  outside  the  range  of  operating 
conditions  actually  used,  these  pressure  ratios  are  low  enough 
to  cause  "choked",  sonic,  constant  mass  flow  rate  through 
both  and  N^.  However,  even  though  there  is  a  low  pressure 
ratio  across  the  inlet  gas  tube,  the  flow  is  everywhere  sub¬ 
sonic.  For  steady  isentropic  flow  from  the  tube  exit  to 
the  larger  tube  exit  area  (1.5  mm  diameter  vs.  0.81  mm  for  ) 

requires  subsonic  flow  at  the  tube  exit  with  a  Mach  number 
l-otween  0.14  and  0.2.  With  viscous  effects  this  Mach  number 
is  even  lower.  The  flow  transit  time  for  a  monatomic  gas  is 
calculated  to  be  ca.  0.5  sec  for  the  oven  chamber.  The  second 
orifice,  ,  measures  C.51  mu  in  diameter.  With  both 

orifices  choked  only  one-fifth  to  one-eighth  of  the  carrier 
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gas  exits  via  A  somewhat  higher  proportion  of  the  metal 

clusters  ought  to  pass  through  and  the  electron  beam,  be- 

4 

cause  of  their  much  greater  "molecular"  weight  (4  x  10  to 
2  x  106)  . 

Diffraction  patterns,  if  observed  at  all,  are  commonly 
visible  for  only  about  10  seconds,  and  never  for  more  than 
about  90  seconds.  We  were  clearly  not  operating  the  oven  under 
steady-state  conditions,  although  they  were  constant  over  the 
time  span  of  an  exposure  for  film  recording  of  a  pattern  (i.e., 
fraction  of  a  second).  The  oven  gas  pressure  and  temperature, 

Pq  and  Tq,  along  with  other  experimental  information,  are 
summarized  in  Table  II. 

Exposed  plates  are  developed  in  the  usual  manner  and 
their  degree  of  blackening  is  quantified  by  scanning  across 
a  diameter  of  the  diffraction  pattern  with  a  recording 
microdensitometer.  Patterns  from  three  plates  are  presented 
in  Fig.  2.  Peak  spacings  (diffraction  ring  diameters),  peak 
widths,  and  peak  heights  are  all  measured  from  these  densito¬ 
meter  records.  Calibration  of  spacings  and  minimum  peak  widths 
is  obtained  by  tracing  the  patterns  obtained  from  a  commercial 
aluminum  diffraction  standard.  Unit  cell  dimensions  are  de¬ 
rived  from  ring  radii  using  the  expression 
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a  =  LA 


A 


2  2  2 
+  k  +  l 


/  (r(l  -  3r2/8L2)  ) 


(1) 


where  L\  is  a  "camera  constant"  determined  from  the  aluminum 
standard,  A  =  0.053551  8  for  the  50  Kev  electrons,  and  L  = 
47.13  cm.  The  term  in  (r/L)  is  a  small  flat-plate  correction. 
Cluster  diameters,  D,  are  estimated  from  peak  widths  using 


D  =  LA  / 


(2) 


where  W  is  the  observed  peak  width  and  Wq  is  the  "instrument 
broadening"  taken  to  be  the  peak  width  found  for  the  aluminum 
standard.  Peak  heights  were  taken  to  be  proportional  to  the 
diffracted  power  in  each  peak  and  are  corrected  for  the  non¬ 
linear  response  of  the  plates.  Estimated  standard  deviations 
are  assigned  to  each  intensity,  oQbs ,  and  multiplied  by  the 
same  correction  factors.  It  should  be  emphasized  that  only 
diffraction  patterns  strong  enough  to  be  seen  by  eye  on  the 
phosphor  screen  were  recorded. 


RESULTS  AND  DISCUSSION 
Source  performance 

In  addition  to  varying  gas  pressures  and  filament  tempe¬ 
ratures,  in  search  of  optimum  conditions  for  the  production 
of  high  densities  of  smaller  clusters,  we  also  explored  the 
effect  of  changing  the  carrier  gas.  Argon  and  helium  show  dif¬ 
ferences  reflecting  their  differing  atomic  masses  and  colli- 

(2  3) 

sion  cross  sections  '  .  Past  experience  with  nucleation 


studies  indicate  that  a  larger  number  of  smaller  clusters 


should  result  from  a  more  rapid  cooling  of  the  hot  metal 
vapor  The  chemically  inert  molecular  gas,  sulfur  hexa¬ 

fluoride  SFg ,  with  its  many  internal  degrees  of  freedom,  i.e., 
its  much  larger  heat  capacity  per  molecule  should  be  a  more 
effective  third  body  in  the  early  stages  of  cluster  formations. 
Since  silver  has  a  relatively  small  affinity  for  oxygen,  C02 
was  also  successfully  used  as  a  carrier  gas. 

A  greatly  increased  production  of  silver  clusters  is  ob¬ 
served  using  molecular  gases.  This  occurs  at  substantially 
lower  filament  temperatures  (i.e.,  lower  vapor  pressures)  as 
indicated  in  Table  II.  Silver  deposits  on  the  oven  liner 
were  collected  on  microscope  grids.  Electron  diffraction  pat¬ 
terns  for  these  samples  verified  the  metallic  silver  structure 
with  no  evidence  for  compounding  due  to  possible  molecular 
gas  decomposition. 

An  argument  was  made  in  an  earlier  paper ^  for  a  rela¬ 
tionship  between  cluster  size  and  the  product  of  oven  pressure, 
PQ,  times  the  temperature  of  the  evaporating  metal  sample,  T  . 
The  results  are  presented  in  Fig.  3. 

Both  employed  the  same  cluster  generator  seen  in  Fig.  1 
and,  whereas  the  prior  investigation  used  one  carrier  gas  (Ar) 
and  several  metals  (In,  Bi ,  and  Pb) ,  this  work  is  for  one  metal 
(Ag)  and  several  carrier  gases  (Ar,  He,  C0_  and  SF.. )  .  It  is 

Z  b 
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seen  that  all  the  results  with  the  filament  source  are  self- 

consistent  with  the  heavier  gases  producing  a  given  cluster 

size  as  progressively  lower  P  T  .  The  results  for  Pb  in  a 

o  ora 

boron  nitride  crucible  required  a  greater  p0Tom  for  a  given 
cluster  size,  but  operated  for  a  longer  time  before  the  metal 
sample  was  depleted. 

Diffraction  Analysis 

From  the  electron  diffraction  patterns,  one  easily  obtains 
crystalline  unit  cell  dimensions  and  crystallite  sizes  from 
the  diffracted  rings'  diameters  and  breadths,  respectively. 
However,  more  detailed  atomic  scale  information  is  contained 
in  the  relative  intensities  of  the  rings.  This  information  may 
be  extracted  by  matching  the  intensities  calculated  from  a 
model  to  the  observed  intensities.  However,  the  calculations 
leading  from  model  to  diffracted  intensities  are  more  complex 
in  the  case  of  electron  scattering  than  for  X-rays  or  neutrons 
because  of  multiple  scattering  or  "dynamical"  effects.  These 
effects  are  large  enough  to  cast  doubt  upon  any  such  analysis 
based  upon  simpler  kinematic  (or  the  first-approximation, 
Blackman- formula ^ )  calculations^^.  Accordingly,  a  computer 
program  has  been  written  which  calculates  the  diffracted  in¬ 
tensities  in  as  complete  and  correct  a  manner  as  time  and  bud¬ 
get  permitted.  An  example  for  silver  is  shown  in  Fig.  4,  nor¬ 
malized  to  the  (3 1 1 ) reflection.  If  the  scattering  were  kinema- 
tical  the  curves  would  all  be  horizontal  with  the  value  at 


=  0.  Variations  from(31l) (also  varying)  by  a  factor  of  2 
to  3  can  be  seen  at  large  Dfc. 

The  corrected  peak-height  data  obtained  from  the  plates 
listed  in  Table  II  are  used  as  imput  to  the  computer  program. 
The  best  fitting  cluster  diameter  D  and  root  mean  square 
amplitude  of  thermal  vibration  U  are  included  in  Table  II. 
Also  listed  are  the  (dis ) agreement  factors  and  R  ,  for  the 
multislice  and  kinematic  calculations  respectively.  These 
are  defined  as 


R  ’  1  '  E(W  ^alc’2  '  %bs  '  "‘WW2’ 


(3) 


In  every  case  R^  is  an  improvement  over  R^,  but  it  must  be 
pointed  out  that  R  is  calculated  using  the  value  of  U  opti- 

mi zed  for  Rw. 

M 

If  this  model  and  calculation  are  to  be  deemed  success¬ 
ful,  strong  correlation  between  several  quantities  in  Table 
II  should  be  seen.  For  example  higher  temperature  clusters 
should  have  greater  U  and  a  for  a  given  size,  and  and 
should  agree.  The  above  do  not  seem  to  correlate,  although 
the  goodness  of  fit  R^  does  vary  with  size  Dg  in  a  systematic 
way  as  shown  in  Fig.  5.  The  best  agreement  occurs  for  the 
largest  clusters. 
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The  variation,  between  the  experimental  and  theoretical, 
12-beam  multiple  scattering  calculations,  for  each  diffraction 
peak  is  seen  in  Fig.  6.  The  goodness  of  fit  is  clearly  better 
as  cluster  size  increased  from  top  to  bottom  (these  three 
experiments  lie  on  the  upper  dotted  curve  in  Fig.  5) .  The 
most  prominent  and  consistant  features  are  the  larger ( 111 ) peak 
and  smal ler ( 400 ) similiar  to  Ref.  5  (Pb  clusters). 

A  variety  of  interpretations,  other  than  the  bulk  fee 

structure,  have  been  considered  qualitatively  but  have  one 

or  more  features  that  do  not  fit  the  observations.  They  include 

(5) 

the  combined  fee  plus  liquid  model  ,  stacking  faults  and 
twinning  between  fee  regions,  decahedral  (pentagonal  bipyra¬ 
mid)  and  icosahedral  sphere  packing  models^'7^,  and  the  pre¬ 
sence  of  vacancies. 

We  have  obtained  improved  agreement  between  observed  and 
calculated  intensities,  but  it  appears  that  a  simple  cubic- 
closest-packing  model  of  the  silver  metal  clusters  studied 
here  is  not  correct.  It  may  be  that  a  partially  liquid  or 
amorphous  structure  is  correct,  but  this  has  not  been  demons¬ 
trated  conclusively.  Single-phase  (homogeneous )  models  other 
than  fee  can  be  tested  with  the  multi-beam  computer  program 
by  replacing  one  subroutine. 
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TABLE  I 

Oven  Operating  Pressures 


Inleta) b) c) d) e) 

Oven 

First 

Second 

Gas 

P 

oo 

torr 

P  /P 
o  oo 

P 

o 

torr 

P,/P 

1  o 

Pumping 

Chamber0^ 

P1 

torr 

Vh 

Pumping 

Chamberd)e) 

P2 

torr 

He 

2.99 

0.10 

0.30 

0.53 

0.16 

<  5xl0~5 

^  1 .6xl0-5 

He 

7.22 

0.14 

1.01 

0.36 

0.36 

$  5xl0~5 

He 

9.92 

0.16 

1.56 

0.32 

0.51 

<  5xl0~5 

Ar 

3.07 

0.14 

0.44 

0.40 

0.18 

<  5xl0~5 

Ar 

5.62 

0.18 

1.03 

0.30 

0.31 

^  5xl0~5 

C°2 

2.88 

0.17 

0.48 

0.38 

0.19 

^  5xl0~5 

♦ 


a)  Gas  supply  at  room  temperature  (22°C)  and  oven  power  off. 

b)  Measured  at  inlet  of  supply  tube  10  cm  load  and  1.5  mm  inside  diameter. 

c)  The  opening  or  nozzle  between  the  oven  and  the  first  pumping  chamber 
has  =  0.81  mm. 

d)  The  opening  between  the  first  and  second  pumping  chambers  has  =0.51  mm 
with  a  separation  between  nozzles  of  4.3  mm  (see  Ref.  1  for  additional 
details) . 

e)  Measured  over  the  15  cm  baffle  and  diffusion  pump. 


SUMMARY  OF  EXPERIMENTS  AND  THEORETICAL  COMPARISONS 
(Silver  Clusters  in  Various  Carrier  Gases) 
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Fig.  1 


Fig.  2 


Fig.  3 


Fig.  4 


Fig.  5 


FIGURES 


Schematic  diagram  is  shown  for  the  metal  cluster  gene¬ 
rator  with  TC  -  thermocouple ,  W  -  prism  mirror  window, 

H  -  heater  for  metal  vaporization,  S^  and  S^  ~  sets  of 
screws  for  configuration  adjustments,  CG  carrier  gas, 
and  N2  -  flow  orifices  or  nozzles. 

Microdensitometer  tracings  of  Experiments  16,  22  and 
30  (see  Table  II)  are  plotted  with  O.D.  the  optical 
density,  log10(I/IQ)  and  S  =  4n\  *sin  (0/2).  Plate  30 
is  displaced  0.5  O.D.  for  clarity  with  Miller  indices 

shown. 

Experimental  cluster  size  D  ,  as  a  function  of  oven 
pressure  and  temperature  is  determined  using  Eq .  (2) 

for  several  metals  and  carrier  gases  with  F-filamcnt 
source  and  C-crucible  source. 

Theoretical  Ag  cluster  Peak  intensities  are  normalized 
to  the (31l) peak  and  plotted  as  a  function  of  Diameter  . 

The  goodness  of  fit  is  seen  to  improve  as  size  n 


increases . 
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SULFUR  HEXAFLUORIDE  EXPANSIONS  IN  AN  ARGON  CARRIER  GAS 
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Sang  Soo  Kim  and' Gilbert  D.  Stein 


April  19,  1980 

ABSTRACT 

Nucleaticm  of  sulfur  hexafluoride  in  an  argon  carrier  gas  has  been 
studied  in  very  small,  supersonic,  adiabatic  nozzle  expansions.  Light 
scattering  is  used  to  detect  the  onset  of  condensation  and  high  energy 
electron  diffraction  is  used  to  ascertain  cluster  structure  for  a  range 
of  operating  conditions  where  mixed  clusters  and/or  binarv  nucleation  car. 
occur.  For  the  limited  range  in  SF,  mole  fraction,  v  =  0.01  to  0.03, 

O  'O 

noth  the  thermodynamic  and  cluster  diffraction  data  indicate  that  the  SF^ 
nucleates  homogeneous lv  with  Ar  accruing  ont:  the  condensed  SF,  after  the 

h 

Ar  sunersa tura ces . 
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Studios  ol  sulfur  hexafluoride  t^F^)  cooled  in  adiabatic  expansions 

have  been  undertaken  in  numerous  research  laboratories  due  to  its  relevance 

with  regard  to  isotope  separation  schemes  employing  UF^  and  as  the  now- 

standard  example  for  the  study  and  illustration  of  multi-photon  absorption 

phenomena.  We  have  evolved  a  multi-faceted  approach  to  the  studv  of 

phase  change  involving  gasdvnamic  studies  and  theoretical  comparisons  using 
.  .  .  1 

the  governing  equations  ol  motion  ,  laser  scattering  experiments  to  monitor 

■1 

nucleation  and  growth  ,  molecular  beam  investigations  to  ascertain  the 
important  operating  parameters  with  regard  to  making  intense  cluster  beams  \ 
and  high  energy  electron  diffraction  for  probing  cluster  si.:e,  temperature 
and  structure.  The  purpose  or  this  eommun ieat ion  is  to  relate  an  unusual 
diffraction  experiment  and  the  conditions  under  which  we  believe  that  mixed 
clusters  of  SF  and  Ar  are  formed. 

P 

The  experimental  apparatus  consists  of  a  conventional,  three-stage, 
differentially  pumped  molecular  beam  which  enters  a  scattering  chamber  where 
it  is  crossed  bv  a  -0  KeV  electron  beam.  The  diffraction  patterns  are  ob¬ 
tained  using  film  in  a  variable  shutter  arrangement  known  as  a  rotating 
sector.4’"5  However  small  Laval  nozzles  an  used  in  place  of  the  usual  ori¬ 
fice  free  jet  source.  Nozzle  7  is  made  from  glass  capillary  tubing  and  has 

a  throat  diameter  D  =  0.13  mm,  exit  diameter  P  =3.0  mm,  and  a  length  of 

e 

2.5  cm  (see  Ref.  1  for  details).  Previous  experiments'5  were  used  to  ascertain 

the  range  of  operating  parameters  required  to  produce  SF,,  cluster  beams. 

Thev  showed  that  5F.  mole  fractions  in  the  range  =  0.03-0.125  were  best 

b  S  Fn 

suited  for  these  studies. 

Diffraction  patterns  are  shown  in  Fig.  lA  for  v  =  0.1  ipure  Ar) ,  '.'Hg, 

^  *  ry 

0.032,  and  0.125  and  denoted  as  Exps.  A  through  !)  respec r i ve  1  v .  The  Ar  pattern 

i— f) 

is  Lace  centered  cubic  (FCC)  having  a  unit  cell  parameter  a  =  5. A  A.  The 

Miller  indices  for  bodv  centered  cubic  (BCC)  SF  are  shown  in  Exp.  D  vielding 

b 

a  unit  cell  dimension  of  a  =  5.S  X.  Ail  these  experiments  are  conducted  with 
the  same  nozzle  and  operating  conditions  with  only  mole  traction  varied.  Tin 


.-ottering  parameter  is 


(V-/\ )sin(Q/2)  where  \  is  the  electron  wavelength 


and  g  is  the  scattering  angle.  The  pattern  ror  Fxn.  Ti  has  on  1 v  the  110  Desk 
ter  SF  while  Kxn.  C  has  some  features  of  i  snnern'-s  i  r  ion  or  linear  comb  i  na  •  i  on 


nr  Exp.  A  and  ,  esoecialle  at  small 
'•’near  ’o  oe  rhe  case. 


However  i: 


’ther  ang  1  os  this  doe  s 


1  order  •  o  understand  the  process  or  tormar  r  t.  ror  ;He  condensed 


’nse  t  1 1  run  Lea  r  t  or 


lid  semt’o  1 


and  Ar  (open  symbols)  are  shown  in  phase  diagram  Fig.  lB.  The  =  0.01 

co  0.03  onsets  are  shown  here  for  comparison  with  Exp.  B  and  C.  The  cooling 
rate,  for  nozzle  expansions,  increases  with  decreasing  throat  size.  Onsets 
of  condensation  plotted  here  for  both  large^and  small  nozzles^  (e.g.  Nozzle  7) 
have  been  obtained  using  light  scattering.  Note  that  two  of  the  small  nozzle 
onsets  for  SF^  cross  the  Ar  vapor-solid  equilibrium  line.  The  corresponding 
Ar  partial  pressures  reveal  that  the  Ar  is  also  supersaturated,  S  >  1  (S  = 
(pv/p00)T  where  pv  is  partial  pressure  of  the  species  in  question  and  pm  is 
the  equilibrium  pressure  at  the  same  temperature  T) . 

Lines  of  constant  undercooling  (&Tp  =  (T^-T)  where  T  is  the  temperature 
of  the  gas  and  T^,  the  vapor-solid  equilibrium  temperature  at  the  same  pressure) 
are  included  for  the  Ar  equilibrium  line.  The  onsets  of  Ar  nucleation  are  for 
a  large  nozzle7  in  a  helium  carrier  gas  having  yAr  =  0.01  to  0.06,  plus  one 
pure  Ar.  Note  that  they  all  fall  along  the  line  AT^  =  20. 

The  conditions  at  the  exit  of  Nozzle  7  for  Exp.  A- D  are  shown  as  the 
partially  solid  symbols  and  corresponding  open  symbols  for  theoretical  solutions 
of  the  equations  of  motion  with  phase  change  and  viscous  effects.^  Note  that 
exit  conditions  for  Exp.  B  and  C  occur  near  the  line  &T^  *  20,  but  this  is  not 
an  onset  condition  since  the  clusters  have  grown  to  a  size  large  enough  to 
produce  strong  diffraction  patterns.  The  theoretical  solution  for  Exp.  C 
predicts  nearly  all  of  the  SF^  to  be  condensed  at  the  nozzle  exit  (ge  =  a>o 
or  Ce  =  1.0,  see  Table  I).  Within  the  nozzle  at  SF6  onset  the  Ar  is  un¬ 
saturated  (S^r  =  0.042)  and  thus  the  nucleation  process  is  homogeneous,  i.e. 
pure  SF^  clusters.  By  the  time  the  expansion  has  reached  the  ATp  =  20  line, 

Ar  has  condensed  onto  the  SF^  clusters  to  such  an  extent  that  it  is  visible 
in  the  diffraction  pattern.  The  homogeneous  nucleation  of  Ar,  as  distinct 
cluster  species  from  that  of  the  SFg,  must  approach  or  exceed  the  £Tp  =  20 
line  before  its  nucleation  rate  can  rise  to  a  value  high  enough  to  provide 
enough  surface  for  efficient  depletion  of  the  Ar  vapor  phase.  Heterogeneous 
nucleation  of  Ar  onto  the  SF&  can  potentially  begin  shortly  after  the  Exp.  C 
expansion  crosses  p^  for  Ar.  (The  energy  of  formation  for  a  cluster  on  a 
heterogeneous  site  is  always  lower  than  that  for  homogeneous  nucleation.) 

The  situation  is  similar  for  Exp.  B  except  that  it  has  a  lower  mole  fraction 

■,  ,  Lower  fraction  condensed  C„  =  0.379  and  therefore  an  even  smaller  mass 

-o  e 

fraction  condensed  at  the  exit,  ge  =  0.0244  compared  to  0.1075  of  Exp.  C. 

A  closer  examination  of  Exp.  C  shows  that  the  stroneest  reflection  :or 
3F  ,  the  110  peak,  is  still  verv  prominent  suggesting  that  the  SIX  crvsta: 


-4- 


structure  is  still  largely  intact.  The  strongest  Ar  reflection.  111,  is 

easilv  discernible  although  unresolved  from  the  SF,  200  and  is  substantially 

o 

smaller  than  the  pure  Ar  111  peak  of  Exp.  A.  Many  peaks  from  Exp.  A  and  D 
seem  to  appear  in  C  as  well  but  not  as  complete  peaks  or  as  linear  super¬ 
positions.  In  particular  the  220  peaks  do  not  superimpose  well.  Thus  there 
are  some  differences  between  that  of  pure  Ar  and  pure  SF^  crystalline 
structures.  Since  the  unit  cell  parameters  typical  of  the  pure  species  Jo 
not  match  up  closely  one  may  not  expect  epitaxial  growth  of  the  Ar  onto 

the  SF  .  Perhaps  there  is  a  transition  laver,  amorphous  in  nature,  sepa- 
b 

rating  the  two.  The  diffraction  pattern  for  Exp.  B  shows  only  one  small 
SFj,  peak  (110)  with  perhaps  just  a  hint  of  the  SF^  200  unresolved  with  the 
Ar  200  and  111.  This  is  consistent  with  the  other  patterns  since  the  SF(j 
mass  fraction  condensed  is  much  lower  than  Exp.  0. 

One  feature  of  these  expansions  which  will  in  general  discourage  the 
formation  of  Ar  condensate  onto  the  SF^  clusters  is  due  to  the  SF(:)  latent 
heat  which  must  be  carried  away  bv  Ar  collisions  before  the  cluster  temper¬ 
ature  can  approach  that  of  the  gas.  Estimates  of  cluster  temperature  ob¬ 
tained  from  electron  diffraction  data^  are  higher  than  that  of  the  gas  bv 
30  to  50°C.  The  higher  the  mole  fraction  the  greater  this  difference  is 
expected  to  be.  Thus  even  though  the  Ar  can  supersaturate  at  some  point 
in  these  expansions,  its  condensation  onto  the  SF^  clusters  is  not  insured. 
Perhaps  this  is  the  reason  that  these  mixed  clusters  are  seen  so  seldomlv 
in  these  experiments. 

The  results  of  the  work  presented  here  mav  be  summarized  as  follows: 

(lb  The  mixed  cluster  diffraction  pattern,  Exp.  shows  some 

features  of  the  crystalline  structure  from  both  pure  Ar  and 

oure  SF,  and  some  features  which  can  be  ascribed  to  neither. 
b 

Exp.  B  shows  only  one  SF0  peak. 

(IT  The  thermodynamic  history  for  these  supersonic  adiabatic 

expansions  reveals  that  the  SF„  clustered  homogeneously  at 

some  point  within  the  nozzle,  crossed  the  Ar  vapor-solid 

equilibrium  line  and  most  probably  condensed  Ar  onto  the 

surtaco  of  the  SF  clusters.  An  undercooling  !'  =  2o  i  . 

"  P 

required  :  or  the  oncet  ot  homogeneous  Ar  nuele.it  iou. 

1  11  The  SF  clusters  are  warmer  than  the  ambient  :  i  ;  .  iud  ano 

a  si  ,;iu*  ic.mr  traction  the  o'ndonsed  nhu  .«•  tus:  e  Ar  :: 

its  tit  :  r  ict.  ion  m.viks  ire  to  v  '.coy  iveyo  ■>.  e  ‘ho  Si 

*  ■'  e  r  a  ::::o  ' ;  co"  i  i  :  :  ,  -e  -o  -  r  •  :  -  r  ■ :  e  ■  .  -  r—a  1  :.v  :  m  i  x  i  ■ 
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(4)  For  the  case  of  the  small  nozzles,  in  which  their  operating 
conditions  are  adjusted  so  as  to  produce  SF^  nucleation  near 
the  exit,  it  is  possible  to  have  the  Ar  in  the  range  ATp  a  20 
at  which  point  Ar  nucleation  could  also  be  initiated.  Thus 
the  possibility  of  a  binary  nucleation  process  arises.  The 
growth  of  this  heterogeneous  cluster  could  produce  diffraction 
patterns  radically  different  from  those  reported  in  this 
investigation. 
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CAPTIONS 


Figure  1.  Electron  diffraction  patterns  are  taken  in  a  molecular  beam 

from  clusters  formed  in  SF^-Ar  mixture  expansions  with  mole 

fractions  that  varv  from  0.0-0.125.  The  molecular  beam 

conf iguration  (D .  the  skimmer  diameter,  D  the  collimator 
^  c 

diameter,  stagnation  pressure  and  temperature  is  the 
same  for  all  these  experiments.  The  vertical  lines  with 
arrows  up  locate  the  SF^  BCC  peaks  and  arrows  down  those  of 
Ar  FCC  structure.  In  the  phase  diagram  pv  is  the  partial 
pressure  of  either  the  SFb  or  Ar  as  appropriate.  The  solid 
symbols  are  for  SFh  onsets  of  nucleation  and  the  open  symbols 
for  Ar  onsets.  The  SF0  data  is  for  =  0.01  to  0.03  for 

a  varietv  of  nozzles  having  different  throat  areas.  The  Ar 
partial  pressure  for  the  lower  three  SFb  onsets  da  )  indicate 
that  the  Ar  is  also  supersaturated.  The  adiabatic  expansion 
line  is  shown  for  Exp.  A-D  with  the  experimental  and  theoretical 
exit  conditions  indicated.  Lines  of  constant  undercooling, 

ATp  =  10  and  20°C  are  shown  for  the  Ar  vapor-solid  equilibrium 
line.  The  two  solid  lines  with  ~  =  1  and  7  =  10^  3re  theoretical  SF 
onsets  computed  for  the  small  nozzles",  using  Che  classical 
nucleation  theory  with  7  a  convection  factor. 
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ABSTRACT 

The  condensation  of  SF.  in  Ar  and  He  carrier  gases  has  been  studied  in 

o 

diverging  Laval  nozzles  with  throat  diameters  in  the  range  of  0.1  to  0.2  mm. 
The  clustering  is  detected  using  laser  scattering,  molecular  beam  intensity 
measurements  and  high  energy  electron  diffraction.  Although  the  nozzles  have 
cooling  rates  in  the  10 '  °C/sec  range,  the  correction  for  nucleation  lag  time 
results  in  changes  which  are  less  than  157„  from  the  steady  state  rate  theorv 
For  the  3"  SF^  in  Ar  the  onset  mav  enter  the  Ar  super sa turat i on  region  with 
the  attendant  possibility  for  binarv  nucleation.  Even  though  these  expansion 
are  not  isentropic,  due  to  nozzle  wall  boundary  Livers,  they  are  nevertheless 
much  more  effective  cluster  sources  than  the  "isentropic"  free  iet  expansions 
The  onset  of  condensation  is  consistent  with  earlier  results  for  much  larger 
nozzles  and  agrees  with  the  classical  Volrrer  theor-.-. 


I  INTRODUCTION 


The  study  and  use  of  Laval  nozzles  for  adiabatic,  supersonic  expansion 
ot  a  condensable  vapor  is  important  both  for  many  technological  applications 
and  fundamental  investigations  into  the  dynamics  of  condensation  via  homo¬ 
geneous  nucleation.  Systematic  ’  tvestigations  date  back  60  years  to  the 
studies  by  Stodoia  of  steam  condensation  in  turbine  nozzles.^  Since  then 
many  other  investigators  have  made  contributions  to  this  field  for  a  wide 
range  of  condensable  vapors,^  ^  Recently  there  has  been  particular  interest 
in  condensing  sulphur  hexafluoride,  SF^ ,  expansions  due  to  its  similarities 
with  UF^  and  a  variety  of  proposed  methods  for  cas  phase  isotope  separation. 
The  use  of  very  small  Laval  nozzles  (throat  diameter  D*  in  the  range  0.025- 
0.25  mm)  as  cluster  sources  for  molecular  ‘'earns  has  been  pioneered  by  the 
Karlsruhe  group  with  current  applications  including  fueling  and  heating  of 
Tokamak  fusion  mach ines^“ ’ ^  and  cluster  research. ^  Investigations 
into  the  design  and  use  of  small  Laval  sources  has  been  undertaken  in  our 
laboratory^  in  order  to  understand  and  optimize  their  design  for  studies 
of  the  physical  properties  of  small  clusters  using  both  high  energy  electron 

diffraction  and  tunable  laser  spectroscopv.  The  use  of  free  jet  sources 

18-20 

tor  tne  above  mentioned  research  techniques  has  been  well  established. 

The  purpose  of  the  research  reported  here  is  to  explore  the  SF  nucleation 

6 

process  in  these  small  nozzles  and  compare  ollt  he3  tte  ore  t  ica  1  pre- 

lictions  of  the  nucleation  theory  as  incorporated  into  the  governing  equations 

in  these  small  nozzles 

ot  supersonic,  adiabatic  flow.  The  cooling  rates  are  in  an  interesting 
range,  intermediate  between  those  of  large  Laval  nozzles  in  which  the  steady 
state  nucleation  theorv  is  valid,  and  those  of  free  jet  molecular  beam  source 
ir  which  the  cooling  rate  is  considered  >v  many  researchers  to  ‘-e  sc  ripiJ 
t  ’■  a '  the  -.ticle.it  ion  the  or-,  is  no  longer  va  ltd.  H.e  e:  ec’s  of  nucle.it  ion  i  ■ 
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t 1 me  on  the  steady  state  nucleation  theory  and  its  impact  on  measurable 
properties  will  be  discussed. 

DESCRIPTION  OF  EXPERIMENT 

The  experimental  arrangement  is  shown  schematically  in  Fig.  1.  The 

SF  is  supplied  from  high  pressure  bottles  of  known  mole  f raction  XsKiK  in 
b 

either  an  argon  or  a  helium  carrier  gas.  (The  SF,  mole  fractions  required 

o 

have  previously  been  determined^).  With  measurement  of  the  stagnation  pressure 
p  ,  temperature  T  and  mass  fraction  ,  the  initial  state  of  the  gas 
before  expansion  is  thus  determined.  The  nozzle  is  installed  in  the 
stagnation  chamber  with  a  static  pressure  tap  at  the  nozzle  exit.  The 
nozzle  contours  have  been  measured  and  the  data  fitted  to  a  seventh  order 
polynomial,  see  Table  I  for  details.  An  argon-ion  laser  tuned  to  514.5  nm 
with  a  0.1  watt  CW  beam  is  focused  over  the  nozzle  centerline  about  1.25  mm 
beyond  the  nozzle  exit.  A  tube  with  2  apertures  of  1  mm  diameter  is  con¬ 
nected  to  a  fiber  optic  which  exits  the  vacuum  chamber  and  is  terminated 
at  an  RCA  7850  photomultiplier  tube.  The  nozzle  exit  diameter  is  2.5  mm 
(2  mm  for  Nozzle  7)  and  the  two  apertures  are  located  so  as  to  cut  off  a 
scattering  volume  from  the  laser  beam  of  about  1  mm  length,  i.e.  a  volume 
well  within  the  nozzle  exit  flow  field.  When  light  scattering  is  not  being 
used  the  double  aperture  is  moved  out  of  the  way  and  the  molecular  beam 
skimmer  moved  up  near  the  light  scattering  position  thus  skimming  off  the 
centerline  flow.  After  a  second  stage  of  differential  pumping  the  cluster 
beam  enters  a  scattering  chamber  where  it  is  crossed  by  a  40  KeV  electron 
beam  for  diffraction  studies.  The  cluster  beam  intensity  is  monitored  with, 
an  ionization  ;:age  detector.  The  electron  diffraction  patterns  are  measured 

initio  channel,  electron  scintillator,  pulse-count  j  ae  detection  svs'en-. 

'!  .•  !u- '  or  beam  i>-  modulated  with  a  chopper  located  between  skimmer  and. 
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collimator  and  synchronous  detection  (i.e.  up-down  counting)  is  used 
to  eliminate  background  scattering.  The  single  channel  detector  is 
programmed  to  accumulate  signal  for  a  given  number  of  chopper  cycles  and 
is  then  stepped  to  the  next  increment  in  scattering  angle  9  and  the  process 
repeated  until  the  entire  diffraction  pattern  is  taken.  The  total  electron 
beam  flux  or  current  is  monitored  continuously  using  a  Faraday  cup  inside 
the  beam  trap.  The  primary  function  of  the  diffraction  apparatus  is  to 
record  patterns  for  structure  analysis  of  clusters,  however  in  the  work 
reported  here  it  is  used  to  estimate  mean  cluster  size  and  temperature. 

LIGHT  SCATTERING  AND  MOLECULAR  BEAM  MEASUREMENTS 

With  the  incoming  laser  beam  polarized  perpendicular  to  the  90° 
scattering  detector  (i.e.  its  electric  field  along  the  flow  direction) 
the  Rayleigh  scattering  intensity  has  been  measured  as  the  function  of  p^ 
for  a  range  of  SF&  mole  fractions  in  Ar,  12.5%  SFfo  in  He,  and  pure  Ar. 

This  experimental  sequence  has  been  used  for  several  nozzles,  three  of  which 
are  reported  here  and  two  are  shown  in  Fig.  2.  The  differential  scattering 
cross  section  per  unit  volume  (dc/dTi)^  is  defined  as  the  average  differential 

scattering  cross  section  per  particle  (dc/dT.)  times  the  cluster  concentration 

3  v 
per  cm  ,  N: 

_ 

(dc/AT)  =  N  (dc/AT)  =  N  1  (dc/d.TI  f  (r) dr  (1) 

V 

o 

where  the  average  differential  scattering  cross  section  is  obtained  by 

integrating  over  a  normalized  size  distribution  f(r) .  The  cross  section 

per  unit  volume  is  related  to  the  incoming  and  scattered  light  intensities. 

2  i 

I  and  I  respectively  and  the  scattering  geometrv  bv  t I  / 1  >  I4R  /”dJ)  - 
y' 1  ■  y  0  o  a 

(dr.d  t  _  where  d  >  is  the  diameter  of  the  aperture  furthes"  from  the  laser  hear. 
an<i  R  is  the  distance  from  this  aperture  to  the  centerline.  The  different:  .  : 


scattering  cross  section  per  cluster  is  the  usual  Rayleigh  or  dipole 

2  A  3  2  o  2  a 

expression  (da/dC.)  ~  a  k  =[  r  (tn  -l)/(m  +2)1  (2~r/\)  where  or  is  the 

polar izah i 1 i tv ,  k  is  the  wave  number  (2v/\) ,  m  is  the  index  refraction, 

and  \  is  the  light  wavelength.  Hie  operational  definition  for  the  onset 

of  condensation  is  taken  as  the  intensity  as  it  rises  above  the  background 

bv  0.1  na ,  which  translates  to  a  value  for  (da/dT.)  of  2.8  x  10  ^  cm 

v 

This  condition  for  onset  is  shown  in  Fig.  2. 

With  the  light  scattering  apparatus  moved  out  of  the  way,  the  nozzle 

can  then  be  moved  toward  the  skimmer  such  that  the  molecular  beam  intensity 

is  a  maximum  and  the  relative  distance  is  close  to  that  for  the  light 

scattering  (3.25  mm  from  the  exit).  The  unclustered  beam  intensities  for 

this  molecular  beam  geometric  configuration  is  just  discernible  above  the 

measured  background  ionization  gage  pressure.  As  Pq  is  increased  beyond 

1.5  bar  (1  bar  =  105  Pa  ^  1  atm)  the  intensity  rises  dramatically  as  a 

21  22 

result  of  clustering  in  the  nozzle.  ’  Since  these  intensity  curves  are 
nearly  linear  (See  Fig.  3)  the  criteria  for  onset  is  taken  as  that  pressure 
obtained  from  a  linear  extrapolation  back  to  zero  intensity.  This  criteria 
is  virtually  the  same  as  taking  a  small  intensity  rise  above  background  as 
done  with  the  light  scattering. 

It  should  be  noted  however  that  the  unclustered  beam  intensity  can 
be  seen  above  the  nac  kg  round  •.an.:!  1 2  •  ho'tzft  'he  si  znal-to-noise  level  is 

low.  In  the  laser  ierect;  ui  met  ho.;  be  in  clustered  flow,  i.e.  molecular 
•  c  u  •  r  o  r ;  r..’ ,  cannot  o  do  toe  :  :  :>■  :  >  *!>e  .  •  ra"  1 :  ah  r  intensity  seen  by 

■ 1  i>  n(  emu  It  i  r»l  :<-r  .  ■;  r  o  '  •  ,ru  ;  i  r  »<n  me  f  hod  is  more  sen  si  five 

■  i ",  ’ :  ■ .  •  !  i  a  -  :  ,•  i  •  t  e  r  :  ;  irr  »  •  u  •  o.-t  :  nc  'he  initial  cluster 

y.V  a •>«>!>  .  j gt.  ••• 

;  ■  u  *  *  '  i*  •  :  1 


pro  airc  in  the  r.iivv 
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of  1-8  b.ir  (high  densities,  really)  the  collision  mean  free  path  is 

much  less  than  the  nozzle  diameter  .with  the  Knudsen  number  Kn  «  \ ,  / D 

fp 

-3-4 

in  the  range  of  10  -10  even  at  the  nozzle  exit.  Thus  these  flows  are 

continuum  or  collision  dominated.  Typical  Reynolds  numbers  for  these 
4  17 

flows  are  in  the  10  range  so  that  boundary  layers  are  expected  to  be 

thin.  Nevertheless  since  these  nozzles  have  very  small  inlets  of  the  order 

of  0.1  mm  and  very  gradual  divergence  angles,  the  viscous  boundary  layer 

can  extend  to  the  centerline  thus  resulting  in  fully  developed  flow.  (Only 

in  the  case  of  pure  St'  expansions  is  there  evidence  of  an  isentropic  core.  > 

Although  isentropic  flows  are  far  simpler  to  deal  with  analv t ically , the 

importance  of  a  more  complete  understanding  of  these  small  nozzle  flows 

necessitates  a  non-isentropic  analysis  including  the  effects  of  viscosity 

and  heat  conductivity.  Sufficient  gasdynamic  measurements  have  been  made 

so  that  this  more  complicated  analysis  can  be  made  and  compared  to  experiment. 

The  details  of  the  nozzle  geometries  used  here  are  given  in  Table  I  and 

tvpical  ranges  in  nozzle  exit  properties  are  shown  in  Fig.  4  for  Nozzle  7. 

SF6 

Two  sets  of  curves  are  presented  spanning  the  range  of/mole  fraction  from 

12.57.  to  07.  in  Ar.  The  exit  Mach  number  increases  with  leveling  out 

for  argon  to  a  value  of  5.5.  At  the  same  time  the  exit  temperature  drops 

with  increasing  p  to  values  as  low  as  30°K  and  exit  static  pressures  are 
o 

in  the  range  of  1-2  torr  (133-266  Pa). 

Cooling  rates  in  geometrically  similar  nozzles  scale  with  the  nozzle 

throat  diameter.  For  a  given  throat  size  the  cooling  rate  is  dependent  op. 

the  variation  in  nozzle  cross  sectional  area  with  position.  The  larger 

the  special  gradient  in  area  the  more  rapid  is  the  cooling  rate.  Tvpic.il 

cooling  rates  are  shown  in  Fig.  5  for  these  cases:  a  nozzle  with  1  cm 
d  i  ame  ti  r  . 

throat  jaUXSa  ,  Nozzle  7  iron  this  work  having  a  throat  diameter  of  0.12' 


-.r. .  ind  a  tree  iot  expan  ;  i  on  wire  l 


ts  rj;  roat  d  i  ire  r  e  r  also  0.12  3  mm. 
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maximum  cooling  rate  always  occurs  near  the  nozzle  throat,  x  =  0,  and  is 

denoted  by  an  asterisk.  Atypical  onset  of  condensation  is  also  noted. 

The  thermodynamic  conditions  at  onset  are  most  important  with  regard  to 

the  nucleation  process  since  it  is  near  maximum  here.  The  cooling  rate  at 

onset  for  the  small  Laval  nozzle  is  about  an  order  of  magnitude  higher  than 

the  larger  nozzle.  (There  is  some  variation  with  specific  heat  ratio  Y 

but  it  is  not  significant  for  the  purposes  of  this  discussion.)  The  free 

jet  expansion  is  not  controlled  by  solid  boundaries  as  are  those  of  the 

Laval  nozzles  and  it  thus  represents  a  maximum  cooling  rate  for  a  given  gas 

and  throat  size.  In  Fig.  5  onset  occurs  at  a  cooling  rate  40  times  that  of 

the  Laval  nozzle  with  the  same  size  throat.  If  the  cooling  rate  gets  too 

high  a  variety  of  non-equilibrium  processes  set  in.  For  example  in  free  jet 

expansions  studies  have  been  made  of  non-equilibrium  rotational,  vibrational, 

23 

and  even  translational  effects.  Thus  it  is  not  surprising  that  there  may 

be  non-steady  state  effects  for  the  nucleation  process.  For  most  nozzle 

flows  with  throats  2  1  cm  and  starting  conditions  £  STP  the  steady  state 

nucleation  theory  obtains,  and  there  are  negligible  effects  due  to  nucleation 
2  small  used  in  the  present  experiments 

lag  time.  The  /nozzles  /Ki£3f&  have  cooling  rates  high  enough  that  a  time- 

dependent  nucleation  rate  must  be  considered. 

These  flows  may  be  assumed  to  be  adiabatic  with  regard  to  the  inter¬ 
face  between  the  gas  flow  and  the  solid  nozzle  walls.  However,  from  one 
point  to  another  within  the  flow  the  temperature  gradients  mav  be  high  enough 
that  heat  conductivity  from  one  part  to  another  is  not  negligible.  The 
effects  of  viscosity,  for  all  of  the  experiments  reported  here  reach  to  the 
centerline  of  the  flow  at  some  point  in  the  nozzle.  These  combined  effects 
result  in  exit  Mach  numbers  which  are  lower  than  the  ideal  case,  with  pressure 
and  temperatures  whicii  are  higher  than  for  isentropic  expansions.  The  cool  in 
t  the  gas  is  nevertheless  quite  substantial,  producing  verv  h :  gh  super  a  t  u.-  i  • 


and  resulting  in  ellicient  cluster  nucleic  ion  and  growth,  i,e.  compared 
to  an  isentropic  free  jet  source  at  the  same  p  and  D* . 


THEORETICAL  CLUSTER  FORMATION  AND  GROWTH 

The  effects  of  boundarv  layers  and  their  growth  to  the  centerline 
has  been  treated  in  a  previous  paper  with  the  comparison  to  a  variety  of 
gnsdvnamic  measurements.^7  Hie  governing  equations  for  fluid  flow,  conser¬ 
vation  of  mass,  momentum  and  energy;  coupled  with  an  equation  of  state;  a 
nucleation  rate  law;  and  a  droplet  growth  equation  are  used  to  predict  all 
the  thermodynamic  properties  throughout  the  nozzle  including  cluster  size 
and  concentration  as  nucleation  proceeds.  The  nucleation  rate  expression 
employed  is  that  of  the  so  called  classical  theorv, 

Jss  =  r4>2<S)l\  eXP<-iG*/kT)  <2' 

3 

where  J  is  the  number  of  critical  size  clusters  formed  per  second  per  cm  , 

k 

and  AG  the  Gibbs  free  energy  of  formation  of  the  critical  size  cluster 

k  *2 

iG  =  4-n-r  ct/3.  The  critical  radius  is  given  bv  the  Gibbs -Thompson-He  Imho  1 1 

k 

equation  r  =  Ccv^'kT.fn  (p^/  P—)  •  The  saturation  ratio  is  defined  as  p^/p^, 

Pv  is  the  partial  pressure  of  the  condensable  vapor,  T  is  the  temperature, 

is  the  vapor-liquid  (or  vapor-solid-)  equilibrium  pressure  at  the  same  T, 

k  is  Boltzmann's  constant,  a  is  the  surface  tension,  m  is  the  mass  of  a 

single  vapor  molecule,  v  is  the  volume  of  a  condensable  molecule  in  the 

c 

condensed  phase  and  “  is  a  correction  factor  used  to  match  theorv  to 
experiment.  The  magnitude  of  the  nucleation  rate  is  ver"  sensitive  to 
surface  tension  (cubic  in  Jie  exponential  expression)  which  is  unknown  :r 
'he  solid  regime  below  the  SF  triple  point  (.123  Ri  but  has  been  estimate-.! 


from  published  data  above  the  triple  point 


a  -  ipnro-t  i  ma  tel  2t>  (  or  : 


'!  <  temperature  dependence  (or  tiie  data  above  rlie  'r:ple  point  is  incornori' 
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inl<»  t!io  extrapolation  below  222  K  giving:  o(t)  =  6 . 7  7  (/*-T/80)  .  The 
value  of  F  is  approximately  10^. 

The  theoretical  predictions  for  the  nucleation  rate  in  Nozzle  7  is 
shown  in  Fig.  6  for  a  range  of  starting  pressures  and  a  mole  fraction  of 
0.03  SF,  in  Ar.  The  time-dependent  nucleation  theory  is  shown  as  the  solid 

V 

lines.  As  p^  is  increased  the  point  of  maximum  nucleation  moves  upstream 
in  the  supersonic  portion  of  the  nozzle  and  rises  to  a  higher  maximum  value 
and  falls  off  more  quickly.  As  p^  is  decreased  the  maximum  nucleation  rate 
decreases  by  more  than  two  orders  of  magnitude  and  spreads  out  in  time  or 
position  through  the  nozzle  by  about  one  order  of  magnitude,  thus  resulting 
in  a  net  cluster  production  or  concentration  which  is  down  by  an  order  ot 
magnitude.  The  calculations  3re  carried  out  for  one  value  of  boundary 
layer  displacement  thickness  at  the  nozzle  exit  and  a  constant  viscous 
dissipation  factor  (A  »  0.6  and  =  .0002  in  Ref.  17).  Thus  the  family 

of  curves  do  not  follow  the  exact  process  in  the  nozzle  as  a  function  of  p  , 
although  it  is  close,  but  does  illustrate  the  trends  qualitatively,  if  not 
completely  quanta tively. 

The  mass  fraction  condensed  is  g  in  grams  SF  /gram  of  mixture,  and 

is  plotted  in  Fig.  7  for  the  same  flow  conditions  as  Fig.  6.  When  all  oi 

the  SF,  is  condensed  out  a  =  a)  =  0.0987  and  this  limit  is  reached 
6  max  o 

before  the  exit  of  the  nozzle  when  p^  r-  A.  Note  the  very  rapid  rise  in  g  a 

high  p  ,  with  the  onset  of  condensation  occurring  sooner  in  nozzle  (lower  x ) 
o 

as  p  is  increased.  For  the  growth  of  the  clusters  bovond  the  critical 
o 

size  the  vapor  condensation  coefficient  is  taken  as  unify  as  usually  done 

tor  condensation  when  the  vapor  is  a  small  mole  traction  in  a  carrier  .:u  . 

I  ;  one  de  sines  the  onset  of  condensation  for  these  im  hul.it  inns  as  that  po : 

wln-re  ■'  is-  equal  to  1"  of  g  =  r  then  the  onset  o:  condensation  at  the 
'  max  o 

:  eg/ ie  exit  will  occur  tor  p  lust  below  1.8  bar. 
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The  two  most  useful  features  of  a  droplet  size  distribution  are  the  mean 
size  and  the  droplet  concentration.  These  quantities,  designated  r^  and 
at  the  nozzle  exit,  are  shown  plotted  in  Fig.  8.  As  expected  the  droplet 
concentration  increases  with  p^  due  to  the  much  higher  peak  nucleation  rate 
as  seen  in  Fig.  6.  As  p^  increases  past  1.5  bar,  nucleation  at  the  nozzle 
exit  rises  to  a  value  where  it  begins  to  significantly  deplete  the  vapor 


phase  and  the  concentration  N 


ins  to  increase.  The  concentration 


N  increases  with  SF.  mole  fraction  as  seen  in  previous  results.  The 
e  6 

concentration  N  is  the  meaningful  quantity  when  comparing  to  light  scattering 
6 

data  since  the  detector  views  the  scattering  from  a  fixed  volume  of  gas, 

independent  of  the  pressure  or  density  in  the  flow  at  that  point. 

The  cluster  concentration  for  37,  SFg  is  also  shown  plotted  as  a  function  ot 

lusters  per  gram  of  mixture,  N  .  Bevond  p  =  4  bar  N  decreases  even  though  the  peak 

ge  '  ro  ge 

nucleation  rate  increases  bv  a  factor  of  4  as  p  goes  from  4  to  7  bar.  The 

ro 

nucleation  rate  is  shut  off  more  quickly  at  higher  densities  due  to  more 

efficient  droplet  growth.  The  total  number  of  clusters  formed  per  cubic 

centimeter  (the  area  under  che  J  vs.  x  Ficure)  still  increases  with  p  but 

-  o 

not  as  rapidlv  as  the  flow  densitv  increases.  Thus  N  decreases.  The 

ge 

maximum  in  occurs  at  the  approximate  pressure  for  which  all  of  the  SF( 

;  -  predicted  to  condense  from  the  vapor. 

The  mean  cluster  radius  at  the  exit  increases  with  increasing  n  and 

o 

i Iso  with  increasing  mole  fraction  o:  Si'  .  Note  also  that  even  while  tie 


:  hii’i't  ;  crow  to  significant  size  betore  'lie  onset  o :  condensation,  comm 
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C rat  ions  in  this  early  part  of  the  flow  are  so  low  that  they  are  not  seen 

in  the  light  scattering  apparatus.  From  the  electron  diffraction  experiments 

it  has  been  determined  that  the  cluster  structure,  like  the  bulk,  is  body- 

centered  cubic  with  two  molecules  per  unit  cell  and  a  unit  cell  dimension  oi 

about  5.8  Angstroms.  Thus  the  droplet  radius  can  be  converted  to  number  of 

molecules  per  cluster,  G,  and  is  also  shown  in  Fig.  SA .  Notice  at  the  higher 

values  of  p  the  cluster  size  varies  from  5  x  10^  for  3%  S !•',  up  to  1.2  x  10 
ro  6 

for  the  12.5%.  These  are  very  large  clusters  so  that  if  clusters  in  the 
range  of  50  to  500  molecules  per  cluster  are  to  be  studied  p ^  must  be  near 
2  bar. 

NON-STEADY  STATE  NUCLEATION  EFFECTS 

As  mentioned  above  most  condensation  experiments  in  Laval  nozzles  with 

throat  diameters  =  1  cm  have  cooling  rates  in  the  10  *  C/ sec  range  where 

2 

corrections  to  the  steady  state  nucleation  theory  are  negligible.  The 
problem  of  nucleation  lag  time  arises  when  the  change  of  state  of  the  gas 
becomes  sufficiently  fast  that  the  cluster  size  distribution  does  not  adjust 
quicklv  enough  to  establish  a  local  equilibrium  distribution  (it  being  a 
metastable  equilibrium  for  the  supersaturated  state!.  Therefore  the  steady 

state  nucleation  rate,  (Eq.  2)  must  be  modified  to  account  for  this  phenomena. 
32 

Abraham  *"  carried  out  detailed  numerical  solutions  for  supersaturated  water 

28  30 

and  tound  that  the  treatment  bv  Collins'  and  Andres  and  Boudart  give  the 

tlosest  predictions  to  his  computer  solutions.  Collins'  prescription  is  used. 

in  this  analysis.  The  nucleation  rate  therefore  becomes  t ime -dependent  and, 

expressed  in  terms  of  the  steady  state  rate,  is 

J  =  .1  (1  -exp(-r/-  >  >  O' 

t  s  s  r. 

we. re  •  is  time  and  -  is  the  nucleation  lag  time  given  bv  the  expression 


T  -  (m*2,  ^4^  fcf  —  <*> 

n  p  \  m  a 

v 

where  p  is  the  density  of  the  condensed  SF,  .  Solutions  to  the  equations 
c  6 

of  motion  including  nucleation  and  growth,  described  in  the  section  above, 

have  been  carried  out  using  both  rates,  Eq.  2  and  Eq .  3.  Computer  solutions 

using  the  steady  state  nucleation  theory  only  are  shown  in  Figs.  6-8  as 

dotted  lines.  The  time-dependent  nucleation  rate  (Fig.  b)  rises  more  slowly 

than  the  steady  state  by  a  factor  of  two,  peaks  out  slightly  downstream, 

rises  to  a  higher  maximum  value,  and  falls  off  further  downstream.  The 

increase  in  maximum  rate  for  the  time-dependent  theory  varies  from  12  to  15 

to  07  as  Pq  varies  from  7  to  3  to  l.S  bar.  At  p^  =  1.8  bar,  there  is  so 

little  mass  condensed  (g  =  0.001  at  the  nozzle  exit)  that  the  two  nucleation 

rates  coincide  past  the  peak.  The  physical  explanation  for  the  qualitative 

to  be  that 

differences  between  the  two  rate  expressions  is  be lieved  3dXSK.'€fXXJ2fJOt3B{JiXl 

as  the  nucleation  rate  is  building  up,  the  steady  state  rate  is  higher  at 

each  position  and  clusters  grow  more  rapidly  than  if  thev  were  nucleated  further 

downstream  via  the  time-dependent  rate.  This  increase  in  cluster 

growth  rate  depletes  the  vapor  quicker  and  thus  shuts  off  the  nucleation  more 

rapidly  than  in  the  time-dependent  case.  This  results  in  a  peak  nucleation 

rate  which  occurs  earlier  in  the  nozzle  and  at  lower  maximum  values. 

This  effect  can  also  be  seen  when  examining  theoretical  predictions 
for  g  as  shown  in  Fig.  7.  For  the  steadv  state  nucleation  case  g  is  somewhat 
higher  than  for  the  time-dependent  case  in  agreement  with  the  relative  ditier- 
ence:,  cluster  growth  rate.  The  differences  are  largest  where  z  changes  most 


rapidlv  and  this  occurs  at  high  flow  densities,  i.c.  higher  p^.  However, 
ter  p  g  '■*  bar  the  theorv  predicts  complete  vapor  condensation  bv  the  tim< 
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the  flew  reaches  Che  nozzle  exit  and  thus  the  difference  between  the  two 

nucleation  theories  would  not  appear  in  any  measurement  at  this  point. 

About  the  largest  difference  at  the  exit  for  the  two  theories  will  appear 

for  those  p  where  g  is  increasing  most  rapidlv,  i.e.,  in  the  range  of 
o  e 

2.4  to  2.6  bar.  The  difference  in  gg  for  2.6  bar  is  0.70  versus  0.675,  a 

difference  of  about  3.5%.  This  is  not  considered  a  significant  difference. 

The  theoretical  solutions  can  also  be  used  to  predict  mean  cluster  size 

(there  is  always  a  size  distribution)  and  cluster  concentration  throughout 

the  nozzle  and  in  particular  at  the  exit.  Several  cases  are  shown  in  Fig.  8. 

The  droplet  radius  increases  rapidly  from  1  to  3  bar  because  also  rises 

rapidly  in  this  region.  Beyond  that  point  the  peak  nucleation  rate  still 

increases  with  p  providing  a  net  increase  in  N  .  As  g  approaches  uu  the 
ro  e  e  o 

higher  cluster  concentration  tends  to  make  the  average  cluster  size  smaller. 

However,  vapor  density  also  increases  with  pQ,  providing  more  vapor  to 

condense  on  the  existing  clusters.  For  the  3%  case,  the  two  tendencies 

nearly  cancel  out,  resulting  in  first  a  small  decrease  in  size  and  then  a 

net  increase.  In  Che  case  of  the  higher  mole  fractions  N  levels  out  beyond 

e 

4  bar.  The  slight  drop  and  rise  in  mean  radius  is  also  seen  for  these  mole 

fractions,  however  beyond  4  bar  the  droplet  r3dius  grows  more  rapidly  than 

the  3%  case  due  to  the  leveling  off  of  .  Two  tilings  should  be  noted  at 

this  point.  First,  the  cluster  concentration  for  the  steady  state  nucleation 

case  is  lower  than  for  the  time-dependent  case  as  seen  for  3%,  and  the  mean 

cluster  size  is  larger.  This  is  completely  consistent  with  the  picture 

when  discussing 

described  for  this  phenomena  HXi/Figs.  6  and  7.  Secondlv,  decreases 

bevond  4  ‘vir  since  g  =  ju  ,  even  though  M  is  increasing  in  this  range, 
e  o  o 

IT  is  shows  that  the  density  of  the  gas  is  increasing  faster  than  the  increase 

:n  cluster  concentration.  Here  again  N  tor  the  steadv  state  rate  is  lover 

go 

t!  in  *he  t  ime— leoenden'  case  above  2.2  bar.  tticv  cross  over  below  this  pei:  ■ 


as  does  N  .  The  greatest  difference  between  the  two  theories  results  in 

a  77.  difference  in  r  and  a  127.  difference  in  N  .  The  difference  in  light 

e  e 

scattering  cross  sections  (da/df})v>  is  also  127.  or  less,  all  of  which  can 
be  considered  as  relatively  minor  changes. 

ONSET  OF  NUCLEATION 

The  onset  of  condensation  using  light  scattering  and  the  molecular 

beam  measurements  are  presented  in  Fig.  9.  The  vapor-solid  equilibrium 

line  for  SF.  is  denoted  as  p  and  provides  an  indication  of  the  amount  of 
o  0® 

undercooling  attained.  The  equilibrium  line  for  Ar  has  been  included  to 

indicate  where  it  is  possible  to  supersaturate  the  carrier  gas  in  addition 

to  the  SF^  vapor.  Data  from  two  larger  nozzles,  with  cooling  rates  as  shown 

2  2 

in  Fig.  5  and  throat  dimensions  of  1.16  x  1.10  cm  and  0.70  x  0.40  cm,  are 

included.  In  contrast,  the  nozzles  for  this  work  with  typical  throat  areas 
-4  2 

of  10  cm  and  higher  cooling  races,  have  nucleation  onsets  which  overlap. 
The  undercooling,  on  average,  increases  with  decreasing  p^  as  expected, 
although  the  scatter  is  larger.  The  detectable  onset  of  nucleation, us ing 
light  scattering,  depends  to  some  extent  on  the  vapor  pressure  at  onset, 
p^>  (°n  vapor  density,  really)  but  translates  to  a  condensate  mass  frac¬ 
tion  of  g,  =  0.0  laj  . 

k  o 

Two  theoretical  onset  curves,  based  on  computer  solutions  and  onset 

condition  g,  ■  0.01  u>  ,  are  shown  for  comparison  to  the  data  The  left-most 
k  o 

curve  was  calculated  using  the  time-dependent  nucleation  theory,  i.e.  the 

c'assical  theory,  Eq.  2  substituted  into  Eq.  3.  This  theoretical  curve 

moves  away  from  the  vapor-solid  equilibrium  line  as  p  decreases  in 

35 

agreement  with  earlier  results  and  is  due  primarily  to  the  increase  in 
*  with  T  but  also  reflects  an  increase  in  nucleation  lag  as  seen  in  Fig.  n . 

As  p  becomes  less  than  0.1  torr, the  theoretical  onset  curves  awav  from 
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the  SF^  equilibrium  line  more  rapidly  and  eventually  crosses  the  Ar 

5 

vapor-solid  equilbrium  line.  Using  the  correction  factor  F  =  10  the 
onset  curve  is  seen  to  be  moved  to  the  right  closer  to  the  equilibrium  line 
and  passing  through  the  data  at  higher  pressure.  The  data  are  seen  to 
follow  the  theoretical  curve  as  pv  is  lowered  although  there  are  a 

few  data  points  further  from  the  theory  on  each  side.  The  two  lowest 
onsets  occur  in  the  supersaturated  region  for  Ar  and  so  could  in  principle 
be  formed  in  a  binary  nucleation  process  involving  both  SF^  and  Ar. 

The  uncertainty  in  the  points  of  onset  derive  primarily  from  two 
features, both  of  which  are  less  severe  as  onset  p  increases.  The  first 
is  an  error  in  temperature  resulting  from  the  adiabatic  assumption  for  the 
centerline  flow  due  to  wall  boundary  layer  effects.  The  second  is  due  to 
some  additional  expansion  beyond  the  nozzle  exit.  The  first  affects  both 
the  light  scattering  and  molecular  beam  measurements .  The  second  feature 
affects  the  molecular  beam  measurements  only  since  the  light  scattering 

volume  is  within  the  exit  "Mach  cone"  at  all  times  (1.25  mm  from  exit)  and 

2 

thus  does  not  sense  the  expansion  occurring  beyond  the  end  of  the  nozzle. 

The  centerline  core  sampled  by  the  molecular  beam  skimmer  (located 
3.25  mm  beyond  the  nozzle  exit)  falls  outside  the  Mach  cone  and  thus  the 
flow  expands  as  the  free  jet,  beyond  the  nozzle  exit,  adjusts  to  the  lower 
ambient  pressure  (0.065  +  0.010  Torr) .  Gasdvnamic  measurements  in  the 
nozzle  and  the  free  jet  reveal  that  the  centerline  flow  has  completely 
expanded  to  the  ambient  pressure  by  the  time  it  reaches  the  skimmer.  Thus 
the  expans  ion  ,  for  the  molecular  beam  measurements,  occurs  first  in  the 
nozzle  and  then  in  the  free  expansion  from  nozzle  exit  to  skimmer  entrance. 
Tin  skimmer  interf  ’one  with  the  flow  is  exnect'd  to  h.>  small  li.e. 
within  the  error  bars  of  the  results  shown  in  Fie.  9)  and  the  flow  i. 
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as sumed  "frozen"  beyond  the  skimmer  entrance  .  The  nozzle  exit  condi¬ 
tions  (as  with  the  light  scattering)  are  determined  with  measured  proper¬ 
ties  listed  in  Table  II.  However,  the  onset  of  condensation  for  the 
molecular  beam  experiments  occurs  at  the  skimmer  entrance.  An  isentrope, 
obtained  with  the  specific  heat  racio,  y,  for  the  gas  mixture  using  room 
temperature  properties,  is  fitted  through  the  exit  state  to  determine 
both  the  conditions  for  onset  and  for  crossing  the  vapor-solid  equilibrium 
line.  Since  the  pressures  at  exit  and  skimmer  are  known,  the  temperature 
determination  requires  only  that  the  flow  be  adiabatic,  not  reversible 
adiabatic  (i.o.,  isentropic).  Here  again  the  incertainty  in  assuming 
adiabatic  flow  is  within  the  error  bars.  Recall,  however,  that  the  molec¬ 
ular  beam  detection  of  onset  is  more  sensitive  than  the  light  scattering 
experiments.  Thus  the  mass  fraction  condensed  at  onset  may  be  considerably 

lower  chan  the  g,  =*  O.Oliu  condition  for  the  light  scattering  data  and  the 
k  o 

theoretical  curves  in  Fig.  9.  Therefore  comparison  of  the  molecular  beam 
onsets  to  an  extrapolation  of  the  theoretical  curves  would  not  be  warranted, 
and  their  location  slightly  "to  the  right"  of  such  an  extension  is  qualita¬ 
tively  expected. 


CONOLCS IONS 


The  following  are  our  conclusions  derived  from  this  investigation  of 
small  nozzle  gasdvnamics  and  cluster  nucleation: 

1.  The  onset  of  SF^  nucleation  in  these  small  nozzles,  having  an  order 

o 

magnitude  or  more  increase  in  cooling  rates,  is  virtually  the  same  as 
for  the  larger  nozzles  in  the  region  where  their  onset  pressures 
overlap . 

2.  As  the  onsets  move  to  vapor  pressures  below  1  torr  it  is  possible 
for  the  nucleation  process  to  occur  in  the  supersaturated  region 
of  the  carrier  argon  gas  thus  introducing  the  possibility  of 
binary  condensation. 

3.  Theoretical  calculations  on  the  effect  oi  nucleation  lag  times 
for  nozzle  cooling  rates  in  the  107  to  108  °C/sec  range  is  not 
an  overriding  feature  for  these  flows.  Estimated  variation  in 
physical  properties  between  the  time-dependent  and  steady  state 
nucleation  processes  range  from  0  to  157„. 

d.  Prediction  of  the  onset  of  condensation  with  experimental  measure¬ 
ments  is  satisfactory  if  a  nucleation  correction  factor  7  ~  10^ 
is  used  in  the  computer  solutions,  along  with  previously  estimated 

parameters  to  account  for  boundary  layer  effects  as  determined 

,  17 

iron  gasdynamic  measurements. 

5.  The  onset  of  condensation  of  SF.  in  a  helium  carrier  gas  is 

h 

greatly  delayed  (i.e.  higher  p  )  over  that  in  argon.  This  is  caused 

(i.e.  higher  kinematic  viscos  i tv  1 
primarily  by  the  greater  vi scous  dis s i pa t ion /occurr i ue  in  these 

nozzles  due  to  wall  effects.  In  spite  of  this  higher  entropy 

nroduc  t  i  or.  ,  those  nozzle  a  can  still  produce  copious  amounts  of 

large  ,'K  .-lu-tor;  and  in  this  sense  are  much,  more  e;  ticient  than 

:  roe  et  mole  cul  ir  ream  source-..  The  improvement  in  nozzle  design 

.  a  tod  -v  this  work .  for  use  with  helium  expansions,  has  beer. 
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incorporated  into  a  series  of  new  nozzles.  Investigations  into 
noble  gas  clustering  in  helium  carrier  gas  expansions  will  be 
reported  soon. 

6.  These  beam  sources  are  capable  of  producing  clusters  in  a  size 
range  all  the  way  from  dimers  and  trimers  up  to  values  in  excess 
of  106  molecules  per  cluster,  with  cluster  densities  at  the  nozzle 

•  •  r  •,  ,  ,„11  -3 

exit  in  a  range  of  1  to  4  x  10  cm 

7.  The  experimental  measurements  for  SF.  onset  using  both  molecular 

o 

beam  and  light  scattering  techniques  are  consistent.  The 
molecular  beam  measurements  are  more  sensitive  in  the  sense  that 
they  reveal  the  clustering  process  XXJtSLXXlM  sooner  than  the 
particular  arrangement  for  the  light  scattering  as  used  here. 

The  ionization  gage  measurements  show  the  clustering  emerging 
from  unclustered  beam  measurements  as  pq  is  increased,  whereas 
molecular  light  scattering  is  not  observed  above  the  background 
stray  light  level. 
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FIGURE  CAPTIONS 

Fig.  1  The  experimental  arrangement  is  shown  in  which  the  incoming  flow 

contains  a  fixed  mole  fraction  of  SF^  supplied  from  high  pressure 
bottles.  Flow  enters  a  stagnation  chamber  containing  nozzle  N 
from  which  p^,  T^  and  nozzle  exit  pressure  p^  are  measured.  The 
laser  beam-LB  crosses  the  flow  centerline  normal  to  the  page  and 
light  scattering  is  measured  at  right  angles  through  a  double  slit, 
fiber  optic-FO  connection  to  a  photomultiplier  tube-PMT.  A  skimmer 
Sc  and  Collimator  C  sample  the  flow  in  a  standard  molecular  beam 
configuration.  The  molecular  beam-MB  is  collected  in  an  ionization 
gage-IG.  The  molecular  beam  with  clusters  is  crossed  with  a  40 
KeV  electron  beam-EB  which  is  collected  in  a  Faradav  cup-FC  which 
is  housed  within  a  beam  trap-BT.  The  total  beam  current  is  moni¬ 
tored  with  an  amntece r-A  and  the  electrons  scattered  at  an  angle  c 
which  is  adjustable  is  collected  through  a  slit  onto  a  scinti  1  lator-S  i 
coupled  to  a  fiber  optic,  photomultiplier  and  pulse  counting 
electronics . 

Fig.  2  The  differential  light  scattering  cross  section  per  unit  volume 

at  the  laser  locat ion- (dc  'dT )  for  40  decree  scattering  is  shown 

v 

as  a  function  of  s tagnat ion  pressure p  for  a  varietv  of  SF,  mixtures. 

o  b 

Pure  argon  expansion  is  included  here  as  a  lower  limit  (0  SF,  mole 

o 

fraction).  Figure  A  is  for  Nozzle  6  and  B  for  Nozzle  7.  Note  the 
scale  change  in  the  ordinate. 

Fi  '  '•  The  nuclearion  experiment  s  usiiv  the  molecular  heam  ionization 

cave  do' e  c  t  or  are  plo'ted  as  a  function  of  p  and  converted  to  in 
molecular  flux  a:  •  he  electron  ‘'earn  intersect  ion. 


Fig.  A 


Fig .  5 


Fig.  6 


F  i  g 
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Typical  flow  properties  at  the  exit  of  Nozzle  7  are  shown  plotted 
as  a  function  of  .  They  include  the  static  pressure  p  ,  temper¬ 
ature  T  and  Mach  number  M_ .  Thev  cover  the  span  in  mole  fractions 
from  0 %  to  12.57..  Note  that  the  exit  pressures  are  in  the  1-2  torr 
range  while  temperature  approaches  30  K. 

Characteristic  cooling  rates  are  shown  for  Laval  nozzles  with  throat 

diameter  equal  to  1  and  .0125  cm  and  for  a  free  jet  of  the  same 

diameter  as  the  small  Laval  nozzle.  The  cooling  rate  at  the  throat 

is  denoted  with  *  and  typical  onset  cooling  rates  indicated  by  an 

arrow.  The  steady  state  nucleation  theory  is  satisfactory  for  the 

lower  cooling  rate  but  must  be  modified  for  the  upper  two. 

The  nucleation  rate  J  has  been  calculated  for  3%  SF,  in  Ar  for 

o 

Nozzle  7  using  appropriate  boundary  layer  parameters  as  determined 
from  gasdynamic  measurements.  The  steady  state  Jgg  is  compared  to 
the  time-dependent  rate  .  Due  to  the  nucleation  lag  time  the  time- 
dependent  rate  builds  up  more  slowly  than  the  steady  state,  rises  to 
a  higher  peak  value  and  decreases  downstream  of  the  steadv  state. 

Also  as  p^  is  decreased  J  begins  further  down  the  nozzle  and  is 
spread  out  over  a  longer  flow  length. 

The  mass  fraction  condensed  g  (grams  condensed  per  gram  of  mixture'! 
is  plotted  as  a  function  of  position  in  the  nozzle  for  the  same  . 
flow  conditions  as  Fig.  6.  For  p  >  3  bar  virtually  all  the  vapor 
is  n rod i c r ed  to  be  condensed  out. 

The  mean  particle  radius  at  the  exit  -  r  plus  cluster  concert  t  r;i  t  i  n- 

ho  *  h  D'-r  unit  "olume  -  N  and  per  uni'  gram  of  mixture  -  N  nr- 

e  \*  <  - 

p !  '-red  for  several  SF  mole  tractions  as  a  function  ol  n  .  1  oi 

h  o 

V.  . .  mnli  fractions  the  onset  of  "u.- 1  eat  t  on.  occurs  at  lower  p 
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antl  tlie  clusters  grow  to  larger  mean  sizes  while  the  number 
density  per  unit  volume  is  inversely  proportional  to  mole  fraction. 
The  steady  state  nucleation  theory  predicts  a  slightly  larger  size 
and  smaller  concentration  for  the  37,  mixtures.  Note  that  cluster 
concentration  per  gram  of  mixture  goes  through  a  maximum  at  = 

4  bar. 

Fig.  9  The  onset  of  SF&  nucleation  has  been  plotted  in  a  phase  diagram 
where  p  is  the  partial  pressure  of  the  SF^.  The  equilibrium 
lines  for  both  SF,  and  Ar  are  shown.  Previous  data  from  "large" 

D 

nozzle  i  are  compared  with  these  results  and  the  theoretical  pre¬ 
dictions  are  shown  as  the  two  solid  lines  .  =  1  and  10  .  The 

svmbols  with  circles  around  them  denote  the  molecular  beam  data. 
Note  that  three  of  the  lower  onset  points  cross  the  Ar  equili¬ 


brium  line. 
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Fig.  1  The  experimental  arrangement  is  shown  in  which  the  incoming  flow 
contains  a  fixed  mole  fraction  of  SF6  supplied  from  high  pressure 
bottles.  Flow  enters  a  stagnation  chamber  containing  nozzle  N 
from  which  p0,  Tq  and  nozzle  exit  pressure  pe  are  measured.  The 
laser  beam-LB  crosses  the  flow  centerline  normal  to  the  page  and 
light  scattering  is  measured  at  right  angles  through  a  double  slit, 
fiber  optic-FO  connection  to  a  photomultiplier  tube-PMT.  A  skimmer 
Sc  and  Collimator  C  sample  the  flow  in  a  standard  molecular  beam 
configuration.  The  molecular  beam-MB  is  collected  in  an  ionization 
gage-IG.  The  molecular  beam  with  clusters  is  crossed  with  a  40 
XeV  electron  beam-EB  which  is  collected  in  a  Faraday  cup-FC  which 
is  housed  within  a  beam  trap-BT.  The  total  beam  current  is 
monitored  with  an  aiwneter-A  and  the  electrons  scattered  at  an 
angle  9  which  is  adjustable  is  collected  through  a  slit  onto  a 
scintillator-Si  coupled  to  a  fiber  optic,  photomultiplier  and 
pulse  counting  electronics. 


1'2  3  4  5  6 


P0  (B«r| 


P0  (8«fi 

The  differential  1  i -th c  scatterinr  cross  section  per  unit  volume 
at  the  laser  location- (dc'/dT.) for  90  decree  scattering  is  .shown 
as  a  function  of  stagnation  pressure  p  for  a  variecv  of  SF 
mixtures.  Pure  arpor.  expansion  is  included  here  as  a  lower  lim. 
(0  S r ^  mole  fraction!.  Figure  A  is  for  Nozzle  r>  and  3  lor  Nozzi 
Note  the  scale  chance  in  the  ordinate. 


4  Typical  flow  properties  at  the  exit  of  Nozzle  7  are  shown  plotted 

as  a  function  of  pQ.  They  include  the  static  pressure  pe ,  temperature 
Tg  and  Mach  number  M^ .  They  cover  the  span  in  mole  fractions  from 
07„  to  12.5'1.  Note  that  the  exit  pressures  are  in  the  1-2  torr  ranee 
while  temperature  approaches  30  X. 


X  (cm) 
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The  nucleacion  race  J  has  been  calculated  for  3Z  SF ^  in  Ar  for  No 
7  using  appropriate  boundary  layer  parameters  as  determined  from 
gasdvnamic  measurements.  The  steady  state  Jss  is  compared  to  the 
time-dependent  rate  J  .  Due  to  the  nucleation  lag  time  the  time- 
dependent  rate  builds  up  more  slowly  than  the  steady  state,  rises 
a  higher  peak  value  and  decreases  downstream  of  the  steady  state. 
Also, as  p  is  decreased  from  7  Bar  J  begins  further  down  the  note 

rO  .  ’ 

anc  is  spread  out  over  a  longer  flow  Lencth. 


The  mass  fraction  condensed  g  (grams  condensed  per  gram  of 
mixture)  is  plotted  as  a  function  of  position  in  the  nozzle  for 
the  same  flow  conditions  as  Fig.  6.  For  p0  >  3  bar  virtually 
all  the  vapor  is  predicted  to  be  condensed  out. 


Torr 
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The  onset  of  SF^  nucleation  has  been  plotted  in  a  phase  diagram 
where  Pv  is  the  partial  pressure  of  the  SF^.  The  equilibrium 
Lines  tor  both  SF^  and  Ar  are  shown.  Previous  data  from  "large" 
nozzles^  are  compared  with  these  results  and  the  theoretical 
predictions  are  shown  as  the  two  solid  lines  ~  =  1  and  10-1.  The 
svmbols  with  circles  around  them  denote  the  molecular  beam  data. 
Note  chat  three  of  the  lower  onset  points  cross  the  Ar  equilibrium 
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ABSTRACT 

An  investigation  of  diverging  supersonic  nozzles  with  throat  diameters 
in  the  range  of  0.025  to  0.25  mm  has  been  carried  out  using  gasdynamic 
measurements  for  a  variety  of  gases  and  mixtures  with  sulfur  hexafluoride. 

Due  to  the  small  nozzle  size,  boundary  layers,  although  thin,  may  consitute  a 
major  fraction  of  the  flow  or,  indeed  at  low  pressures,  may  be  fully  developed. 
Nevertheless  these  nozzles  have  been  shown  to  be  very  efficient  for  production 
of  clusters  condensed  from  the  expanding  gas.  They  produce  orders  of  magni¬ 
tude  increases  in  molecular  beam  intensities,  relative  to  the  conventional, 
"isentropic"  free  jet  sources.  The  effects  of  gas  properties  and  nozzle 
design  on  performance  have  been  analyzed  and  compared  to  theoretical  solutions 
to  the  governing  equations  of  motion  including  condensation  due  to  cluster 
nucleation  and  growth. 


PACS  Numbers 

1. )  47 . 55 . Ea 

2. )  47 . 55 .kf 

3. )  82 .60 .Nh 

4. )  36 . 90 ,+f 


’  r-  -  ■  t  :  id  re :  ■  Mi  Jw<  s  >.  o  C*<rp .  .  ij  ie.i  >  .  ill:,'';1-  "  ’ !l’  • 

’r  add r  1  s  C  l"  -vl  vania  Lahorar  or'-' ,  3a  1  m.  Ma-i>ac!a>«.- 1  ts 


INTRODUCTION 


The  study  of  small  Laval  nozzles  has  been  undertaken  in  our  laboratory 
as  part  of  a  research  investigation  into  the  properties  of  microscopic 
atomic  or  molecular  particles.  These  nozzles  constitute  the  source  or 
first  stage  of  a  molecular  beam  apparatus  which  will  provide  a  continuous 
stream  of  clusters  in  a  high  vacuum  environment.  It  is  well  known  that 
when  a  gas  supersaturates  in  an  adiabatic  expansion  the  amount  of  under¬ 
cooling  that  occurs  before  the  onset  of  condensation  depends  on  the  cooling 
rate  in  the  expansion.  The  slower  the  expansion,  the  smaller  the  amount 
of  undercooling  for  cases  where  the  expansion  is  continuum  or  collision 
dominated.  For  supersonic  expansions  the  cooling  rate  scales  with  the 
throat  dimension.  Thus  larger  nozzles  will  form  droplets  earlier  in  the 
expansion  (i.e.  at  lower  supersaturations  and  higher  temperatures)  than  for 
small  nozzles.  Free  jet  or  uncontrolled  expansions  into  a  tiigh  vacuum 
environment  represent  the  most  rapid,  in  fact  the  upper  limit,  cooling 
rate  attainable  in  an  adiabatic  gas  expansion.  Thus  it  is  expected  that 
a  free  jet  expansion  will  produce  far  less  condensate  than  a  controlled 
Laval  nozzle  expansion  of  the  same  size.  This  is  in  fact  the  case  as 
has  been  demonstrated  previous ly . 1  ^ 

The  use  of  small  nozzles  has  become  attractive  for  a  number  of  im¬ 
portant  technological  applications.  One  method  of  uranium  isotope  sepa¬ 
ration  involves  the  adiabatic  cooling  of  U Ff  to  as  cold  a  temperature  as 

o  6  7 

possible  for  use  of  tunable  spectroscopic  selection  of  isotopes.  ’  There¬ 
fore  snail  nozzles  are  potentially  useful  in  that  thev  delay  nucleation 
due  to  their  rapid  cooling  rates  and  yet  maintain  relatively  high  cluster  densities 
compared  /W’"\uo  uncontrolled  free  jet  expansions.  Another  application 
which  is  being  implemented  in  Germany  involves  a  small  Laval  nozzle  ex¬ 
pansion  tn  produce  clusters  of  hydrogen  to  be  used  for  fusion  machine 

8  0 

t  tie  lie and  heating.  *'  A  prototype  ol  this  scheme  is  currently  under 
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i  evaluation.  A  third  use  for  these  small  nozzles  in  conjunction  with  a 

high  vacuum,  molecular  beam  type  of  configuration,  would  be  to  form  small 
metal  or  semiconductor  clusters  for  production  of  thin  film  devices. ^ 

The  prospect  of  producing  thin  films  with  unusual  or  more  uniform  properties 
is  the  goal  here.  Finally,  there  are  a  number  of  applications  for  small 
nozzles  in  basic  research.  Primarily  they  can  be  used  as  an  efficient  source 
of  small  clusters  in  order  to  study  the  physical  properties  of  the 

clusters  themselves  or  to  use  them  in  the  study  of  atoms  or  molecules  adsorbed 
onto  them. 

It  is  perhaps  ironic  that  the  development  of  molecular  beam  sources  has 

gone  through  a  couple  of  full  circles.  The  initial  sources  in  the  1930's 

were  orifice  or  free  jet  expansioij)|  at  such  low  pressures  that  they  were 

effusive  in  nature.  Then  in  1951  it  was  suggested  that  the  source  pressure 

be  increased  and  that  the  orifice  be  replaced  with  a  hypersonic  nozzle  to 

make  use  of  the  so-called  Mach  focusing  in  order  to  increase  beam  intensity. 

to  some  problems  with  regard  to  viscous  effects  and  flow  blockage,  in  the 

particular  configuration  considered,  this  technique  was  dropped  in  favor 

of  high  pressure  free  jet  expansions  which  eliminated  the  viscous  effects 

but  retained  the  advantages  of  Mach  focusing,  This  method  has  been  used 

up  to  the  present  time  as  a  beam  source  for  a  great  variety  of  basic  and 
13-15 

applied  research.  Recently  with  the  advent  of  greatly  increased 

interest  in  the  study  of  small  aggregates  in  the  gas  phase,  the  Laval 
nozzle  has  reemerged  as  an  important  device.  The  design,  however,  is 
substantially  different  from  that  originally  proposed. 

As  an  examp le,  cons ider  the  enormous  increase  in  cluster  beam  intensity 

w  i  t  H 

from  a  Laval  nozzle  as  compared  /  a  free  jet  source  (with  the  same  throat 
diameter) .  The  relative  beam  intensity  of  a  pure  argon  cluster  beam  is 
shown  in  Fig.  1.  By  the  time  the  source  pressure  is  raised  to  p  =3  bar 
(1  bar  =  lO^Pa  s  Newton/m  =  0.99  atmosphere  =  7fc0  torr)  the  beam  intensity 
is  off  scale  for  the  ionization  gage  detector  used  in  th i s  experiment.  On 
the  other  hand,  eve  i  at  p  =11  bar  the  free  jot  a.  am  intensitv  is  so  low 
that  it  was  not  usable  for  our  electron  beam  diffraction  experiments.  Hi  us 
it  becomes  obvious,  both  from  theoretical  consideration  )•-,  well  as  experi¬ 
mental  evidence,  that  tile  controlled  expansion  ot  a  small  Laval  nozzle  i  ■ 
fir  superior  to  that  of  the  tree  jet  with  respect  ui  roakinu  cluster  neurits . 


Due 


That  this  groat  enhancement  in  beam  intensity  exists  has  been 
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demonstrated  in  several  laboratories.  However,  in  order  to  understand 
the  dynamics  of  the  expansion  in  some  detail,  to  be  able  to  predict  when 
condensation  will  begin  for  a  given  gas  or  gas  mixture,  and  to  determine 
how  to  design  nozzles  for  future  applications,  a  series  of  gasdynamic 
studies  has  been  undertaken.  For  example,  it  is  known  that  solid  bounda¬ 
ries  will  grow  boundary  layers  which  for  normal  flow  situations  would  be 
thin  at  these  pressures.  However  even  thin  boundary  layers  (of  the  order 
of  0.1  mm)  can  constitute  a  significant  fraction  of  the  nozzle  radius  or 
even  the  entire  radius  (i.e., fully  developed  flow).  Nevertheless,  even 
though  there  are  viscous  effects  in  small  nozzles,  there  are  conditions 
in  which  the  clustering  is  still  far  superior  to  the  "inviscid”  free  jet 
expansions.  How  important  and  under  what  flow  conditions  are  viscous 
effects  important  or  dominant  in  the  small  nozzle  expansions?  Is  there 
an  isentropic  core  in  these  expansions  and  if  so  under  what  circumstances? 

Is  it  possible  to  design  nozzles  which  use  helium  as  a  coolant  or  carrier 
gas  when  the  viscous  effects  in  helium  flows  are  much  greater  than  for 
heavier  species,  that  is, species  with  much  greater  collision  cross  sections? 
Is  it  possible  to  optimize  the  design  parameters  such  as  throat  size, 
divergence  angle,  nozzle  length  and  contour  to  produce  clusters  of  given 
concentrations  and/or  size?  Since  pumping  capacity  in  high  vacuum  systems 
is  always  a  problem,  how  small  a  nozzle  can  one  use  and  still  get  sufficient 
beam  densities? 

■XPFRIMENT.AL  CONFIGURATION 

The  apparatus  built  to  test  nozzles  is  shown  in  Fig.  2.  It  is 

connected  to  the  laboratory  foreline  manifold  which  is  pumped  by  two 

mechanical  pumps  in  series  with  a  large  enough  capacity  so  that  the 

pressure  downstream  of  the  nozzle  exit, p^, is  in  the  ranee  of  0.01-0.1  torr . 

The  nozzles  are  axisvmmetric  and  made  from  glass  capillarv  tubing  with 

i  ntranc  diameters  D  =  0.021-0.1  mm  with  exit  diameters  D  =  1-4  mm. 

o  e 

.4  ver  1 1  <li  J  fiT'-nr  nozzle  extensions  are  used  to  test  the  eifeet  of  nozzle 
lenerh  on  performance.  A  static  pressure  tap  is  installed  in  these  ex- 
'  n-ior.s  in  order  to  avoid  drilling  holes  in  the  class.  The  static  pressure 
■r  'Tie  nozzle  exit,  p1  ,  is  measured  using  l  It1  torr  capacitance  manometer. 
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Stagnation  chamber  pressure  p  and  temperature  T  are  also  measured. 

o  o 

Gas  is  supplied  from  bottles  either  as  pure  vapor  or  as  mixtures  of  known 
mole  fraction.  Downstream  of  the  nozzle  exit  a  Pitot  tube  is  mounted  on 
a  pressure  transducer  and  manipulated  with  x-y-z  motion  capability.  Thus 
the  total  pressure,  Pq2>  can  be  measured  as  a  function  of  position  either 
along  the  flow  direction  or  normal  to  it.  Tq  is  always  near  room  temper¬ 
ature  and  p  is  varied  from  1  to  11  bar.  Since  it  is  known  from  the 
o 

molecular  beam  experiments  that  condensation  will  occur  in  these  nozzles 
as  Pq  is  increased,  an  argon  ion  laser  and  associated  optics  are  included  in 
order  to  detect  the  onset  of  nucleation  just  beyond  the  exit  of  the  nozzle. 

Results  from  four  nozzles  will  be  discussed  in  this  paper  and  their 
contours  are  shown  in  Fig.  3.  The  diameter  has  been  expanded  by  a  factor 
of  8  to  highlight  the  differences.  Nozzle  7  looks  like  a  cut-off  version 
of  nozzle  number  4  with  respect  to  the  downstream  contour  but  has  a  larger 
throat  diameter.  Nozzle  1  has  the  lowest  divergence  angle  near  the  throat, 
9^,  and  Noz2le  6  the  greatest.  Nozzle  1  opens  up  more  rapidly  than  the 
others  toward  the  exit.  Nozzle  7  has  the  smallest  exit  diameter,  2  mm, 
and  Nozzle  6  is  the  shortest.  In  order  to  get  an  idea  of  the  actual 
contour  Nozzle  7  is  also  plotted  on  a  1:1  scale  showing  that  the  expansion 
is  indeed  very  gentle.  Two  angles  are  used  to  characterize  each  nozzle 
defined  as  follows: 


_  _  -1  f dD 

0T  -  tan 


dx/ x=o 

1  /D  N 

-1  '  e  o', 


9n  3  tan  •  L 


(1) 

(2) 


where  L  is  nozzle  length. 

Note  that  none  of  these  nozzles  has  a  converging  section.  This  is 
done  mainly  to  eliminate  boundary  layer  growth  through  the  minimum  section 
frhroatl  of  the  nozzle.  Also,  since  thev  have  relativolv  high  exit  Mach 
numbers,  conditions  at  the  entrance  up  to  Mach  1  will  not  significantly 
influence  conditions  at  the  exit.  Three  different  exit  configurations  an 
shown  in  Fig.  3.  They  h3ve  the  effect  of  reducing  nozzle  length  going 
: rom  A  to  0,  with  the  additional  length  beinR  5  mm,  0 . 5  mm,  and  0.  The 
•eenie r r i ca 1  contours  of  the  nozzles  are  measured  using  precision  drill  rod  - 
as  plug  -age-..  These  data  are  used  in  a  computer  program  to  provide  a  best 
;  i :  using  i  seven! h -ordt  r  polvt-.omia  1 
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D(x)  =  T.  a.x  , 


(3) 
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where  the  coefficient  and  other  pertinent  dimensions  for  the  nozzles 
are  given  in  Table  I . 


MEASUREMENTS  AND  GASDYNAMIC  INTERPRETATIONS 

The  basic  experimental  procedure  is  as  follows: 

(1)  Install  a  particular  nozzle  and  nozzle  exit  extension. 

(2)  Choose  a  particular  gas  mixture  to  be  used  for  a  complete 
set  of  measurements. 

(3)  Run  at  constant  pressure,  measure  pQ,  T  ,  p^,  p^  on  centerline. 

(4)  Record  Pitot  pressure  traces  as  a  function  of  x,  y,  or  z. 

(5)  Take  intermittent  measurements  of  light  scattering  at  nozzle  exit. 

(61  Change  pq  and  repeat  3-5  above. 

(7)  Change  gas  or  gas  mixture  and  repeat  2-6. 

(8)  Change  nozzle  and/or  nozzle  exit  extension  and  repeat  2-7. 

It  would,  of  course,  be  desirable  to  have  measurements  along  the  complete 
nozzle  contour,  but  this  would  have  entailed  an  enormous  design  and  in¬ 
strumentation  program.  Also  since  the  primary  interest  is  in  the  flow 
out  of  the  nozzle  exit,  it  is  the  properties  at  this  point  which  are  of 
greatest  concern.  The  data  are  summarized  in  Table  II. 

The  first  level  of  analysis  makes  use  of  the  above  listed  measurements 

in  conjunction  with  straightforward  one-dimensional  gasdynamic  analysis. 
The  single  most  important  parameter  for  characterizing  flow  conditions  is 
the  Mach  number  at  the  nozzle  exit.  The  Mach  number  is  calculated  in  three 
wavs.  If  one  assumes  that  there  is  an  isentropic  core  then  p^  and  p^  will 
vield  a  Mach  number  via  the  well  known  relation 


—  =  1+ 

P„ 


2  s/ 


v 

Y-l 


(4) 


where  v  is  the  ratio  of  specific  heats  and  M^  is  the  Mach,  number  at  the  exit 
usinv.  static  pressure  measurement.  The  second  method  ot  calculatinc  an  exit 
Mach  >'.umher  is  to  us<-  the  Pitot  pressure  measurement,  lb-re  the  flow  must 
also  a.-  issnpi  d  to  be  isentropic  ton  the  c-  n ter  1  i no )  so  that  the  total 
pr-ssure  ahead  of  th-  Pitot  probe  is  the  same  a*  p  .  With  th i s  assumption 
Mac)'  number  can  ">e  obtained  with  th--  supersonic  Pitot  equation: 
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2  v-1  2  v-1 

/2yMp-  Y+l'  / 2+ (Y - 1 ) M~ 

v  y  +  1  /'  V  ,  ..M2  / 

(y+l)M 


(5) 


The  third  method  for  calculating  exit  Mach  number  is  to  make  use  of  the 
so  called  Pi  cot -Rayleigh  equation.  This  requires  only  local  measurements, 
chat  is,  the  static  pressure  at  the  exit  and  the  centerline  Pitot 
pressure  p^9  at  the  same  location  in  the  nozzle  as  the  static  pressure  tap. 
This  method  does  not  rely  on  any  previous  history  of  the  centerline  flow, 
that  is,  it  does  not  have  to  be  isentropic.  Since  the  nozzle  angular 
gradients  are  very  small,  treating  the  flow  as  one  dimensional  is  a  good 
assumption  and  thus  p^  is  also  the  centerline  static  pressure.  The  Pitot- 
Ravleigh  Mach  number  is  obtained  using: 


1 

!i_  _  l'2yli''r+l\'1 

p02  "  V  Y  +  1  ) 


( _ 2 _ 

V(Y+1)mJ/ 


(b> 


The  exit  Mach  number  calculated  in  the  three  ways  described  here  are 
shown  plotted  in  Fig.  4  for  Nozzle  7  as  a  function  of  p^  for  argon  expansions. 
Results  for  the  long  nozzle  extension  A  /aishown  as  the  dotted  curves  and 
the  short  nozzle  extension  B  as  solid  lines.  First  of  all, note  that  both  M 

s 

and  M  are  greater  than  the  Pi Cot-Rayleigh  Mach  number  This  is  a 

reflection  of  the  fact  that  the  centerline  or  core  flow  is  not  isentropic. 

When  comparing  the  results  of  the  long  and  short  nozzle  extensions  we  see 
that  the  static  Mach  numbers  agree  but  the  Pitot  and  the  Pitot-Ravleigh  do 
not.  It  it  is  assumed  that  the  additional  nozzle  length  results  in  in¬ 
creased  viscous  dissipation,  then  the  larger  values  for  M  and  the  smaller 
values  of  M^  are  qualitativelv  consistent.  It  is  interesting  that  the  static 

Mach  number  M  is  the  same  in  both  cases.  In  Fig.  4  (bl  we  sec  the 
s 

r •.•suits  using  the  same  nozzle  for  sulphur  hexafluoride,  $F  .  Here  the  Mach 
numbers  computed  in  all  three  ways  are  in  fairly  close  agreement.  This  can 
'<■  interpreted  to  mean  that  there  is  relatively  little  viscous  dissipation 
on  the  centerline  or  core  flow,  i .  e .,  that  the  core  flow  is  isentropic. 

!  t  is  also  noteworth v  to  point  out  ttiit  all  of  tin-  Mach  number  cure  - 
it.  r.urlv  tilt  hevond  n  -  1  '->.ir .  This  i  .  norma  1 1  v  i  n  t  •  rnr"  t  ed  to  mean 
•.hi’  d.ir"  I  .veri  i •  I:.  nozzle  are  not  ci'.ann :  i  •  with  t  nervasi nr 

or.”-  nr.  md  •  hu  •.  t  h  i :  tin  r  •  ■  i  ;  indeed  an  isentropic  core  with  size  no! 
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varying  appreciably  with  pressure.  However,  since  the  three  Mach 
numbers  do  not  agree,  except  in  the  case  of  SF^,  there  must  have  been 
some  entropy  production  between  the  stagnation  chamber  and  the  nozzle 
exit.  Beyond  PQ  =  3  bar  there  is  some  entropy  increase  in  the  flow  due 
to  latent  heat  release  due  to  condensation,  except  in  the  case  of  pure  SF^ 
and  He.  The  discrepancy  between  the  Mach  number  based  on  geometric  area 
ratio,  and  the  other  Mach  numbers  is  due,  of  course,  to  the  boundary 
layer. 

Ln  order  to  obtain  a  more  complete  picture  of  the  nozzle  flow, 

Pitot  traverses  perpendicular  and  parallel  to  the  flow  direction  are  taken. 

Pitot  traverses  across  the  flow  just  beyond  the  nozzle  exit  (i.e.  along 
the  v  directions  with  the  z  =  0,  see  Fig.  3)  are  shown  in  Fig.  5  for  two 
of  the  nozzles  fitted  with  Che  long  nozzle  extension  A.  The  Pitot  tube 
has  an  0D  of  0.5  mm  and  an  ID  of  0.25  mm  and  is  set  0.76  mm  from  the 
nozzle  exit.  Figure  5A  shows  the  results  for  both  helium  and  argon 
expanded  through  Nozzle  6.  The  first  thing  to  be  noticed  is  that  wherever 
there  are  flat  sections  of  the  profile,  they  constitute  relatively  small 
fractions  of  the  exit  diameter.  For  the  argon  the  first  evidence  for  a 
flat  part  of  the  Pitot  (and  therefore  the  velocity)  profile  appears  at 

=  2,5  bar.  The  helium  on  the  other  hand  must  be  expanded  from  a  stagnat ion  pressure 
of  nearly  8  bar  before  a  level  part  of  Che  Pitot  trace  is  seen.  The  level 
section  is  interpreted  as  a  flat  velocity  profile,  i.e.  a  region  outside 
the  boundary  layer.  It  is  also  seen  that  the  flat  velocity  profile  disap¬ 
pears  for  the  argon  at  higher  pressures.  This  is  believed  due  to  the  massive 
condensation  occurring  at  these  pressures  and  the  possibility  that  the  droplets 
are  large  enough  to  traverse  the  stand-off  shock  wave  and  enter  the  probe 
before  complete  evaporat'on. 

The  Pitot  traverse  measurements  at  the  exit  for  Nozzle  U  are 

shown  in  Fig.  5B.  The  most  important  difference  between  Nozzle  d  and  the 

others  is  that  its  throat  diameter  is  smaller  bv  a  factor  of  2  or  more. 

Thus  viscous  effects  can  reach  the  centerline  from  the  wall  more  readily. 

Here  it  is  sc’-n  that  even  starting  pressures  as  high  as  10  bar  produce  no 

flat  pt'otile.  Also  notice  that  the  magnitude  of  the  Pitot  pressure  is 

with 

around  5  torr  as  compared  /  120  tore  in  Nozzle  t ,  in  enormous  difference. 

!'v  t.  the  ireoo  must  no  expanded  up  to  nearlv  8  Bar  before  a  flat  velocity 
prof  il(  >;  obtained.  He  call  that  the  downstream  or.ib  it  ~t  pressure,  in 
"i  •  3,  is  in  the  range  of  n.  1-0. 01  torr  and  so  t.vre  is  an  additional 
exports t.  '■  be'.-ond  the  nozzle  exit.  Also,  as  the  tree  jet  continues  down¬ 
s’  r.  cl  the  nozzle  exit,  there  v.'l  bo  a  Wine  zone  with  the  ambient  gar. 
t’: Jt  traces  alone  the  nozzle  cet.t'-rline  show  a  v.’ r  i ..  C  i  or  in  pressure 
•  *'  r’:<  lower  PtS'ip’1  i  or  ,.raon  and  at  all  p  re  asur.  s  for  h-lium. 

'  ■  '  ’•  '  i’ifo1  pr  m-  •  I  1  r  ■  shown  at  sever  if  x  loc.iti.t-,  iron  i'e 
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■  xit  of  the  nozzle  downstream.  The  flow  is  argon  through  Nozzle  1, 
exit  configuration  3  (i.e.  .Nozzle  1-B)  with  the  exit  defined  as  0.5  mm 
beyond  the  end  of  the  glass  nozzle  (see  Fig.  3).  Many  of  the  Pitot 
profiles  with  this  configuration  show  pressure  bumps  or  "horns"  on 
either  side  of  the  centerline.  This  type  of  feature  has  been 
observed  in  supersonic  boundary  layer  flow  as  the  Pitot  probe 

approaches  the  solid  boundary. ^  Under  these  circumstances  the  shock 
wave  standing  off  the  Pitot  probe  can  interact  with  the  boundary  layer 
causing  separation  upstream  of  the  probe  thus  changing  the  flow  field 
seen  by  the  Pitot  tube.  This  explanation  is  consistent  with  all  of  the 
Pitot  traverses  that  we  have  made  in  the  course  of  these  studies.  The 
horns  have  been  seen  only  when  the  probe  is  within  the  nozzle  or  in  the 
neighborhood  of  the  exit,  i.e., close  to  solid  boundaries.  It  also  disap¬ 
pears  in  nozzle  flows  where  pQ  ^  3  bar.  As  the  probe  moves  out  of  the 
nozzle  the  horns  get  smaller  and  then  disappear  by  the  time  the  probe  is 

1.5  mm  downstream  fi.e.,x/D  =  0.6).  Notice  also  that  the  centerline  Pitot 

e 

pressure  is  dropping  as  the  probe  is  moved  further  downstream. 

The  symmetry  of  the  horns  is  related  to  the  direction  the  probe  is 
moving,  that  is^way  from  or  toward  the  solid  boundary.  When  the  probe 
direction  is  reversed  the  asymmetry  of  the  Pitot  profile  is  reversed. 

The  Pitot  profiles  are  not  taken  as  continuous  traverses  across  the  flow> 
but  as  discrete  points  waiting  a  sufficient  time  at  each  point  to  insure 
that  no  time  lag  appears  in  the  measurements. 

An  extensive  set  of  Pitot  profiles  has  been  taken  for  a  great 

variety  of  flows  in  several  of  the  nozzles.  An  example  is  provided  in 

Fig.  7  for  Ar  and  Ar-SF^  mixtures  in  Nozzle  7-B.  The  Pitot  probe  is 

located  0.5  mm  into  the  nozzle,  at  the  plane  between  the  glass  nozzle  and 

the  exit  extension.  (See  Fig.  2  or  3).  The  effect  of  the  SF  on  the 

0 

boundary  layer  growth  can  be  seen.  At  p^  =  3  bar  the  flat  portion  of  the 
profile  is  quite  email  in  the  Ar  flow  and  gets  larger  as  the  percentage 
oi  SF  is  increased.  Thus  the  SF.  reduces  the  boundary  lavor  thickness. 

b  h  ‘ 

(.Although  the  flows  presented  here  were  all  steady,  under  some  circumstance 
especially  with  tilth  SF^  mole  frictions,  the  flow  became  unsteady  when 
probed  by  th.-  Pitot  'ube.)  As  expected,  the  increase  ir.  the  size  of  the  eo 
t  low  with  starting  pressures  is  evidenced  in  ail  thru  sets  of  data. 

Because  the  P.  tor  probes  mist  bo  smaller  •  iian  t!  enarxeterisr ic  f  lot- 
d  tamt  I  or  and.  since  the  i :  y  in  the  flow  near  the  exit  is  filling  rapid!'. 


-10- 


chere  is  some  concern  as  to  whether  there  might  he  significant  corrections 

Co  the  one-dimensional  Pitot  probe  theory  assumed  for  l-.qns.  (5)  and  (6). 

It  is  known  that  due  to  viscous  effects  Pitot  tube  measurements  in  the  flow 

field  will  be  different  than  that  of  the  true  total  pressure  if  the  Reynolds 

number  of  the  probe  Re  is  less  than  100  (based  on  probe  diameter  and 

properties  just  ahead  of  the  probe  but  behind  the  detached  shockwave).  In  the 

range  100  >Re^>10  the  measured  value  faLls  to  about  857,  of  the  theoretical 

value  and  then  rises  above  the  theory  for  Re  <  10.  In  order  to  test  tli  is  effect  it 

P 

is  necessary  to  probe  a  flow  of  known  properties.  The  obvious  choice 

for  this  investigation  is  the  use  of  free  jet  orifices.  The  theory  for 

free  jet  expansions  has  been  verified  many  times  and  with  a  variety  of 

1 8 

different  types  of  measurement.  Pitot  data  have  been  taken  with  two 
different  orifices  by  different  investigators  more  than  a  year  apart. 

hey  are  converted  to  Mach  number  using  Eqn.  (5)  and  compared  in  Fig.  8 

19 

to  the  theoretical  Mach  number  given  by  Ashkenas  and  Sherman: 


AV  d  J  2Vy-i;/laV  d  .• 
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where  A  and  x^/D  are  constants  dependent  on  y  with  the  following  values: 
for  y  =  5/3,  A  =  3.26,  xq/D  =  0.075;  for  y  =  1.4,  A  =  3.65,  x  /D  =  0.40. 

The  earlier  data  are  from  an  orifice  with  throat  size  was  0.155  mm  for 
argon,  helium  and  nitrogen  while  the  later  data  are  for  argon  only  using  an 
orifice  of  0.65  mm  diameter.  The  error  bars  on  the  earlier  data  are  due  to 
an  uncertainty  in  the  probe  location  with  regard  to  the  exit  of  the  orifice. 

It  is  seen  jicwever  ,  that  the  data  fit  the  theory  quite  well , but  are  slightly 
lower  than  the  theory  except  for  the  helium  data  without  error  bars  (i.e., 
large  orifice  experiments.)  This  small  deviation  mav  be  due  to  condensation 
effects.  All  data  presented  in  Fig.  8  are  for  high  enough  starting  pressures 
P  so  that  Re  >  100,  i.e.  no  pressure  corrections  are  required.  At  lower  p 

op  O 

the  deviation  between  measurements  and  theory  can  be  reconciled  with  the  Pitot 
correction.^-'  Thus  the  Pitot-Ravleigh  Eqn.  (6)  can  be  used  with  confidence  if 
Rep  >  100  and  all  cf  the  results  presented  in  this  paper  meet  this  criterion. 

The  exit  Mach  number  has  been  determined  for  numerous  gases  and  nozzles. 

The  tvo  most  important  factors  affecting  the  value  of  the  Mach,  number  are  dynamic 

viscosity  u,  and  specific  heat  ratio  v.  Increasing  a  and  decreasing  v  will 

lower  the  Much  number.  Also  as  o  is  increased  latent  hear  due  to  ccndon- 

o 

cation  ’..ill  decrease  the  Mach  number  slirrhtlv.  The  r-  -t-lts  for  Nozzle 
■'.re  -resented  in  Fin.  9A  and  show  a  Mach  number  olatenu  with  increasin'.' 


pressure  for  A:,  Ar-dr,  mixtures,  ar.-j  bF 

6 


At  tb"  hither  'treasures  :b 


Mach  number  tor  Ar, 


3  and  it".  SF  -Ar  mixtures  all  occur  at  a  value 


about  .  .  5.  IT.  us  in  this  range  the  cour." ernclin  •  e: f-c-  %  ol  deert-a: 
and  decreasing  viscosity  with  i  ncreasin.:  :'F  :rol.  ma  ,’i  ,  plus  V 
ol  Iieat  addition  jus'  cancel  out.  The  pur  >  SK  would  •  1  i ; • .  1  • 


lower  Mach  number  on  the  basis  of  y  reduction  only,  but  the  reduced 

viscous  effects  in  these  flows  result  in  a  thinner  boundary  layer  thus 

raising  the  Mach  number.  The  net  result  is  a  decrease  in  Mach  number 

to  about  5.0.  The  12.5%  SF^-Ar  appears  to  be  dominated  more  by  the 

y  and  heat  addition  effect  than  by  the  effects  of  lower  viscosity,  resulting 

in  an  exit  Mach  number  of  4.  The  increased  viscous  effects  in  the  helium 

flows  are  evident  by  their  much  reduced  Mach  numbers  in  spite  of  the  fact 

that  helium  has  a  high  y.  Neither  the  He  nor  the  12%,  SF,-He  mixture 

o 

reaches  a  plateau. 

Consider  now  the  results  of  Nozzle  1  relative  to  Nozzle  7.  Careful 
inspection  of  the  contours  in  Fig.  3  reveal  that  Nozzle  1  has  a  slower 
area  expansion  near  the  throat  but  a  more  rapid  increase  in  the  latter  part. 

Thus  small  changes  in  viscous  effects  between  Ar  and  the  3  and  6%  Ar-SF. 

D 

mixtures  produce  larger  differences  in  Nozzle  1  due  to  the  slower  opening 
throat  area.  This  results  in  higher  exit  Mach  numbers  with  increased  mole 
fraction  of  SF^  up  to  6%  (See  Fig.  9B) .  However,  the  helium  mixtures  perform 
consistently  better  in  Nozzle  1.  The  small  section  near  the  throat  is  about  as 
detrimental  as  it  is  in  Nozzle  7,  however  the  downstream  section  with  a  more  rapid 
opening  will  reduce  the  viscous  dissipation  in  this  region.  The  viscous  effects  in 
the  argon  expansions  are  minimal  in  this  section  for  both  nozzles.  The 
12.5%,  SF^-He  mixture  reaches  a  Mach  number  plateau  above  4  bar.  Thus  the 
data  provide  some  indication  of  what  nozzle  features  influence  gas  ex¬ 
pansions  with  different  viscous  effects  (i.e,  Reynolds  number  Re  =  puD/p,  where 
o  is  density)  and  specific  heat  ratios. 

In  addition  to  the  Mach  number,  other  dimensionless  parameters  are 
calculated.  Reynolds  number  is  used  to  characterize  viscous  effects, 

Knudsen  number  (and  its  inverse  Kn  is  an  indication  if  and  when 
rarefied  gasdvnaraics  is  of  importance,  and  Prandtl  number  Pr  indicates  the 


relative  importance  of  thermal  to  viscous  boundary  layers.  Prandtl  number 


is  Pr  =  j.c  / k  where  c  is  the 
P  P 

thermal  conductivity  and  does 


constant  pressure  specific  heat  and  k  is  the 
not  vary  greatly  through  the  nozzle  or  with 


chances  in  p  .  Tvpical  values  for  Ar  are  0.65  to  0.7  and  for  SF,  0.44 

°  -1  6 
to  0.48.  Representative  values  for  Re  and  Kn  as  a  function  of  p^  tor 

.'.ozzl'  7-B  are  shown  in  Fig.  10.  Revnolds  and  Knudsen  numbers  (Kn  = 


\  'D  whet'  \  is  the  mean  free  path)  range  over  two  orders  of 


nnnnituJ: .  This  is  true  for  all  SF. -Ar  .but  He  and  SF  -He  have  about  an 

O  h 
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order  of  magnitude  greater  variation.  There  is  less  than  a  factor  of 
2  change  with  pressure. 

If  relationships  can  be  established  for  boundary  layer  effects,  they  could  then 

be  used  to  predict  many  of  the  nozzle  characteristics  without  having  to 

conduct  gasdyr.amic  experiments.  One  such  correlation  is  given  in  Fig.  11. 

The  boundary  layer  thickness  6  is  estimated  here  using  the  Pitot  pressure 

traverses  across  the  exit  of  the  nozzle.  Core  flow  is  taken  as  that  part 

of  the  profile  which  is  flat  and  the  remainder  is  asstlmed  to  be  boundary 

layer.  The  displacement  thickness  5^  is  estimated  using  to  get  an 

o  2 

isentropic  area  A  and  A-/A  =  r  /(r-6.)  .  The  results  provide  the 
I  G  X  cl 

magnitude  of  the  boundary  thickness,  with  5  >  5^  as  it  should  be,  and 

5  =  r  (fully  developed  flow  at  the  exit)  for  Re  <  4,000.  Pitot  and  static 

pressure  data  for  several  nozzles  and  gases  are  converted  to  Mach  number 

which,  when  compared  to  geometric  Mach  number,  is  related  to  the  displace- 

e 

ment  thickness.  When  plotted  as  a  function  of  Re 2  (related  to  laminar 
boundary  layer  growth)  times  Pr  (ratio  of  thermal  to  viscous  boundary  layer 
thickness)  the  results  fall  along  lines  having  a  monotcnic  variation  with 
kinematic  viscosity,  v  =  p/p  (see  Fig.  12). 

COMPARISON  TO  SOLUTION  OF  GOVERNING  EQUATIONS 

A  complete  two-dimensional  solution  to  the  equations  of  motion  with 
viscosity  and  heat  conduction  along  with  nucleation  and  growth  of  the 
condensed  phase,  with  its  attendant  latent  heat  release,  would  constitute  an 
enormous  undertaking  both  in  terms  of  man-years  and  computer  funds.  The 
authors  feel  this  is  not  warranted  in  terms  of  the  goals  of  this  work 
which  require  an  understanding  primarily  of  the  centerline  flow.  A  one- 
dimensional  approach,  incorporating  some  of  the  two-dimensional  features 
i or  momentum  and  energy  transport,  has  been  developed.  It  includes  the 
important  phvsical  phenomena  for  the  flow  process  but  this  compromise 
requires  chat  several  features  be  charactei ized  by  parameters.  They  are 
adjustable  but  must  be  chosen  in  a  self-consistent  and  physically  meaningful 
wav.  The  particular  approach  chosen  here  is  not  unique  in  the  sense  that 
is  the  only  wav  to  pursue  a  simplified  description  ot  this  process.  It 
has,  however,  provided  some  theoretical  guidance  and  insight  with  modest 


•iac  of  computer  tire. 
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Detailed  computations  for  the  case  of  37.  SF^Ar  expansions  have  been 
carried  out  for  Nozzle  7.  The  one-dimensional  steady  state  equations  for  conservatio, 
of  mass,  momentum  and  energy  have  been  employed  with  an  additional 
term  for  viscous  effects  in  the  momentum  equation  and  a  term  for  heat 
conductivity  within  the  gas  in  the  energy  equation; 


pAu  =  m 

Continuity, 

(8) 

dphpudu  =  dPv 

Momentum, 

(9) 

2 

d(h-hy)  =  dq+dPk 

Energy, 

(10) 

where  p  is  density,  A  is  nozzle  cross  sectional  area,  u  is  velocity,  m  is 
mass  flow  rate  and  is  constant  for  steady  flow,  p  is  pressure,  h  is  enthalpy 
and  q  is  the  latent  heat  release  per  unit  mass  due  to  phase  change,  i.e. 
dq=Ldg  where  L  is  the  latent  heat  for  the  condensible  (SF&  here)  and  g  is 
the  mass  condensed  per  mass  of  mixture  (carrier  gas,  condensible  vapor  and 
condensate).  The  quantities  dP^  in  the  momentum  equation  (9)  and  dP^  in  the 
energy  equation  (10)  are  "Production"  terms  due  to  the  transport  of  momentum 
(viscous  effects)  and  energy  (heat  conduction)  to  the  centerline  flow. 

Density  and  enthalpy  are  for  the  mixture  including  the  condensed  phase  if 
present,  i.e.  p^p^+py+p^  and  h=h£+hy+hjc  =  cp^T+cp  T+cT.  The  subscript  i  is 
for  the  inert  carrier  gas  (Ar  or  He) ,  v  is  for  the  condensible  vapor  (SF^) 
and  k  for  the  condensed  phase.  The  specific  heat  at  constant  pressure  is 
Cp  and  c  is  the  condensed  phase  specific  heat.  The  perfect  gas  equation 
of  state  is 

p  *  (l-g)pRT/^  State,  (11) 

where  T  is  temperature  and^g  is  the  molecular  weight  of  the  gas  mixture 
which  changes  as  condensation  begins  due  to  vapor  depletion,  thus  (l-g)p=p  . 

The  boundary  conditions  at  the  solid  wall  are  the  no  slip  velocity 
condition  and  no  heat  exchange,  i.e.  an  adiabatic  wall.  However,  due  to 
temperature  gradients  in  the  boundary  layer,  energy  transport  occurs  within 
the  gas  resulting  in  heat  flow  from  the  periphery  inwards.  For  those 
situations  where  the  boundar.'  layers  grow  together  making  a  fully  developeJ 
flow,  the  cwo  transport  terms  are  brought  into  use  in  the  governing  equations. 

The  term  dPv  in  Eqn.  (9)  is  then 

_  '} 

dPv  u  ~  d>:  (12 

:r~ 

('see  e.c.  Ref.  20  Mqn.  2-2ua)  unJ  dp^  in  Kqn.  H1  '*  bt  o  s 
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(see  e.g.  Ref.  20  Eqn.  7-5)  where  p,  is  Che  dynamic  viscosity  and  k  is  the 
thermal  conductivity. 

Since  only  the  centerline  flow  is  of  interest  in  molecular  beam 
applications,  solutions  of  the  entire  two  dimensional  flow  field  are  not 
desired  and  would  represent  a  substantial  effort  to  attain.  Hie  viscous 
and  heat  conduction  "production"  terms  of  Eqns.  (12)  and  (13)  are  essential 
for  a  physically  meaningful  description  of  the  centerline  flow.  An  integral 
approach  is  used1  in  which  the  production  term  in  each  volume  element 
dV  =  2nrdrdx  is  summed  over  the  flow  volume  cross  section  slice  between  x 
and  x  +  dx  (see  Fig.  3)  and  divided  by  the  slice  volume  to  give  an  average 
value  (for  either  dPv  or  dp^) 

_  f,rw 

dp  =  j  2nrdpdrdx  :  2nrdrdx  .  (14) 

.i  .  , 

o 

In  order  to  evaluate  Eqn.  (14)  the  functions  u=u(r),  T=T(r),  u=u(T)  and 
k=k(T)  must  be  specified. 

Data  for  the  temperature  dependence  of  the  viscosity  and  the  thermal 
conductivity  of  Ar2^  and  SF ^22  are  fitted  with  linear  functions.  The 
velocity  distribution  is  assumed  to  be  parabolic,  where  fully  developed, 
going  from  zero  at  the  wall  to  a  maximum  at  the  centerline  (subscript  c): 

u(x,r)  =  uc(x)(l-r2/r2)  (15) 

where  rw  is  the  radial  distance  to  the  wall.  The  temperature  profile  is  also 
assumed  to  be  parabolic  rising  from  the  value  at  the  centerline  Tc(x)  to 
an  adiabatic  wall  temperature  Taw.  The  use  of  a  temperature  recovery  factor 
of  the  form  r  =  (Taw-Tc) / (T0-Tc )  is  frequently  used  in  viscous  flow  where 
T0  is  the  stagnation  temperature  =  295k.  For  a  Prandtl  number  near  unity  (0.7) 
r  =  0.86  giving  Taw  =  0.14  Tc(x)  +  25 8.  The  temperature  distribution  then 
becomes 

T(x,r)  =  Tc(x)  +  AT(r2/rJ)  ,  (16) 

where  AT  =  Taw-Tc. 

The  average  value  of  the  transport  or  "Production"  term,  Eqn.  (14), 
for  viscositv  and  heat  conduction  can  now  be  evaluated  and  they  are 

_  - 2u  dx 

dP.  =  - o —  V^c  +  “  dT  -T  > 

r~ 
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and 


dpk  =  pir  L  2  \  c  +  I  37  ^;+  ~  i  Uc+  g  ^  (18) 

c  c  rw  rw 

whore  dy/dT  and  dk/dT  are  obtained  from  the  straight  line  curves  fitted 
through  the  temperature  dependent  y  and  k  and  averaged  for  the  particular 
mole  fraction  of  Ar  and  SFg  at  that  x  position. 

Since  the  radial  gradients  dT/dr  and  du/dr  are  small  near  the  flow 
centerline, the  production  terms  dPv  and  dp^  in  Eqns .  (9)  and  (10)  should  be 


a  small  fraction,  call  it  F 


we  can  write 


and 


dP 


dP, 


of  the  volume  averaged  term  over  all  r.  Thus 

dP.,  » 


a 


u. 


(19) 


The  use  of  the  integral  technique  eliminates  the  need  for  the  actual  velocity 
profiles  so  that  dp^fdp^)  is  not  sensitive  to  the  details  of  the  velocity 
(temperature)  distribution.  Typical  values  of  F^  used  to  calculate  the 
transport  properties  on  the  centerline  are  in  the  range  0.0001  to  0.0005. 

The  point  at  which  fully  developed  flow  begins  and  ends  can  be  varied  in 
the  program. 

It  should  be  stressed  here  that  the  radial  distribution  of  properties 
are  used  only  to  compute  the  value  of  the  production  terms  dPv  and  dp^ 

(Eqns.  (17)  and  (18)  substituted  into  Eqn.  (19))  which  appear  in  Eqns.  (9) 
and  (10).  The  solution  to  the  governing  equations  of  motion,  Eqns.  (8), 

(°),  (10) are  carried  out  with  the  use  of  an  equation  of  state,  Eqn.  (11), 
plus  droplet  nucleation  and  growth  equations  described  bv  Eqns.  (21)  and 
(22).  They  are  solved  as  a  closed  set  of  one -d imens ional  (x)  flow  equations 
with  phase  change  and  transport  terms  for  viscosity  and  heat  conduction. 

The  boundarv  laver  effects  are  handled  in  (he  usual  wav“*  with  a  displacement 
fhiekness  !t|  which  crows  from  zero  thickness  at  the  nozzle  entrance  to  a  valut 
at  the  exit,  'i,  characterized  hv  a  non-dimensional  thickness  '  =  «i.'r  = 

' 1  ‘  Wf 

2  ■,)/!>  ,.  The  values  for  the  centerline  properties  used  in  calculating  dp  and 
d.P'  are  taken  as  those  obtained  from  the  one-d  imens  ioua  <  solutions. 

rhr  variation  ot  displacement  thickness  witli  \  is  provided  for  using  an 
exponential  growth  equation  giving  idle  equivalent  mvzl.  diameter  as 

V 

iii  >  =  2v,  ,ix  l  -  )  —  (2''i 


J 


in  ■  ) 


t  2>"i 


-lf>- 


whoro 


nii  .i.iMiaai 


where  xe  is  the  distance  from  nozzle  entrance  to  exit  and  N  is  varied 

typically  from  1  (i.e.  linear)  to  2. 

In  addition  to  the  gasdynamlcs,  equations  to  describe  the  droplet 

nucleation  and  growth  are  also  Included .23-25  quC  to  tjie  re i a 1 1  ve  1  y  high 

cooling  rates  for  gases  passing  through  these  small  nozzles,  a  correction 

to  the  steady  state  nucleatioo  theory  is  incorpora  ted  The  steady 

state  nucleation  theory  is  applicable  for  nozzle  cooling  rates  of  up  to 

106  °C/ sec.  Cooling  rates  of  10^  °C/scc,  typical  of  free  jet  expansions, 

are  thought  to  be  too  high  to  be  dealt  with  using  the  nucleation  rate  theory. 

7  8 

The  nozzles  in  this  paper  have  cooling  rates  in  the  range  of  10-10  where 
adjustments  to  the  steady  state  theory  are  expected  to  apply  (Correction'* 
by  a  factor  of  2-4  at  most). 

The  steady  state  nucleation  theory  is  a  function  of  the  thermodynamic 
variables  obtained  in  conjunction  with  the  solution  to  the  governing 
equations  (8)- (11)  , 


'Pv\2/?.t\ 


VkT 


w.1  rrK  -^G 


.rrnv 


(2D 


J  Is  the  number  of  critical  size  (*)  clusters  formed  cm*J  sec*1  and  AG'  = 

*2 

4rrr  o/3  is  the  Gibbs  free  energy  of  formation  of  a  critical  size  cluster, 
i.e.  the  size  which  is  large  enough  to  become  a  stable  droplet  of  condensed 
SF^.  Here  py  is  the  SF^  vapor  pressure,  k  is  Boltzmann's  constant,  a  is  the 
SFb  surface  tension,  m  is  the  mass  of  an  SFf,  molecule  and  ^  is  the  volume 
of  one  SF(j  molecule  in  the  condensed  phase.  The  critical  radius  is  r  = 

2ovk  'kT  Zn  S,where  S  is  the  SF^  saturation  ratio,  S=(pv',pVi=l-j.  with.  pVo)  the 
vapor-liquid  equilibrium  pressure  at  temperature  1. 


Since  the  cluster  size  Is  small  compared  to  the  mean  free  path 
of  the  gas  and  since  the  condensate  (SF^  is  a  small  fraction  f 37.)  of 
the  argon  carrier  gas^the  droplet  growth  law  is  that  obtained  from  ele¬ 
mentary  kinetic  theory  with  a  condensation  coot' :  1  cii-nt  o  near  1.0. 

Thus  the  growth  law  is 

dO/dt  =  a.1-'  ,  <2 

where  0  is  the  droplet  area  it-.'mr"  and  J  is  th-'  SK()  impingement  rate  pet 
".nit  area,  -p,.  '2i— .v.kD  .  Taking  proper  a  .-court  o’  tie  nutf'er  ot  droplet 


torme-.i  c .  'dx.pl u:  the  growth  •»;  preview  lv  nucleated  clusters,  the  tuts 
tract  tor  condense  •'  can  *'o  computed,  g -g ;  t;1  ,  and  since  dq-id*.  h’qr.s .  (21''  -  ... 
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(22)  provide  dq  and  g  in  Eqns.  (3)  and  (4). 

To  summarize:  Eqn.  (19)  provides  the  production  or  transport  terms 
when  boundary  layer  effects  penetrate  to  the  flow  centerline,  Eqn.  (20) 
provides  the  effective  nozzle  area  A,  Eqns.  (21)  and  (22)  provide  dq  and 
g,  and  h=h(T).  Therefore  Eqns.  (8) -(11)  can  be  considered,  after  the  above 
substitutions,  as  4  equations  in  4  unknowns,  p,  p,  T  and  u. 

Due  to  the  high  cooling  rates  in  these  nozzles  , and  densities  typicallv  falling 

less  than  1%  of  the  starting  densitv,  vibrational  relaxation  for  the  SF 

33  6 

molecule  is  included  using  published  relaxation  times.  The  program 
calculates  all  the  thermodynamic  variables,  gasdynamic  properties  including 
velocity  and  Mach  number,  details  of  the  phase  change  including  nuc lea  Cion 
race,  number  of  droplets  formed,  and  size.  The  primary  concern  here  is 
to  match  the  gasdynamic  properties,  i.e.  p^  and  P^t  These  data  for  the 
37,  SF -Ar  in  Nozzle  7  are  plotted  as  the  circles  in  Fig.  13.  As  p  is  increased 

D  O 

the  boundary  layer  thickness  8,  the  displacement  thickness  8d,  the  centerline 

production  dPv  and  dP^,  and  the  length  of  centerline  viscous  effects  x* 

are  all  changing.  The  calculation  procedure  holds  these  parameters  constant 

while  p  is  varied  and  thus  the  theoretical  results  represent  "planes" 
o 

crossing  the  plane  of  the  experimental  data  in  a  hvpi rspuee. 

Upon  examining  many  sets  of  computer  results  using  a  systematic 

variation  it:  parameters  a  reasonable  and  consistent  picture  of  the  nozzle 

flow  is  obtained.  The  Pitot  pressure  data  seen  in  fig.  13  (b)  are  intersected 

with  lines  of  constant  It  is  found  that  variation  of  N,  F  and  x*  al! 

have  virtually  no  effect  on  p^*  Thus  Pitot  pressure  is  determined  or.ee  A 

is  fixed  for  a  given  inlet  p  .  The  variation  of  4  from  J.05  to  0.5  as  p 

o  o 

Increases  from  3  to  8  bar,  shown  indirectly  in  Fig.  13(b)  and  explicitly  in 
rig.  14(a),  is  qualitatively  correct  since  as  the  density  increases  the 
transport  of  momentum  and  energy  from  the  nozzle  wall  is  decreased  giving 
thinner  boundary  la-’ors .  Below  3  bar  viscous  effects  (F  )  are  important  as 
well  as  \  for  tin’  values  ot  p  . 

Static  pressure  p^  is  much  lower  than  and  cl>. . e.  -os  in  most  ot  the 

listed  pa  rare  tors  produce  small  fractional  chain--*'-  in  the  static 
pressure.  i.fbe  accuracy  of  the  experiment  il  date,  is  -  .02  tetr.) 
variation  in  ooundat"  later  growth,  X  in  I  m.  1 l’1'!  ,  n-ailts  in  a 

c:  luge  o',  slope  ter  pi  vs  rq,  ’• n  l  :  •; .  1’  ~ae  best  mat  eh  to  the  s  lore  ot  tin- 

data  occurs  for  N=1.3  ■••h-.ch  was  then  used  :  or  ill  tin-  z.i.ail.i.ion-  pro-cu:  d 
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hore.  The  magnitude  of  tor  a  given  p0  is  then  determined  by  fixing  F^. 

The  Pitot  data  below  2  bar  fall  below  the  calculated  values.  This  is  believed 
to  be  due  in  part  to  a  drop  in  measured  Pitot  pressure  from  that  of  the 
true  Pitot  pressure  at  low  probe  Reynolds  numbers  (i.e.  less  than  100)  which 
are  encountered  in  the  flows  at  low  pQ.  Even  when  corrected  they  fall  below 
the  curves  for  6  =  0.5-0. 6  rWg.  In  this  range  the  viscous  effects  are  quite 
prevalent  since  boundary  layers  get  thicker  at  lower  He  and  at  1  bar  an 
increase  of  F  to  0.0008  still  does  not  match  the  Pitot  data.  An  upper 
value  for  F,  used  here  is  set  at  that  value  which  drives  the  Mach  number 

subsonic . 

In  addition  to  the  theoretical  comparison  to  pressure  measurements 
shown  if.  Fig.  13  the  results  for  A  and  M^  are  summarized  in  Fig.  14.  Dis¬ 
placement  thickness  A  varies  from  0.5  to  0.65  with  exit  Mach  number  leveling 
off  at  5.4.  The  computed  results  are  numbered,  with  the  conditions  for  the 
program  given  in  Table  III.  Since  the  program  is  run  with  discrete  choices 
for  parameters,  the  fit,  although  close,  is  not  expected  to  be  perfect  in 
every  case.  The  somewhat  lower  values  of  A,  compared  to  those  obtained  from 
measurement,  is  qualitatively  correct  since  A  is  based  on  boundary  layer 
displacement  thickness  while  the  measurements  in  Fig.  11  are  velocity  boundary 
layers  and  thus  6  >  5^  as  expected. 

DISCUSSION 

Flows  at  higher  pressures  have  thinner  boundary  layers,  anc[  smaller 
centerline  viscous  dissipation  occurring  onlv  in  the  early  part  of  the 
expansion.  Tne  boundary  layer  becomes  fully  developed  shortly  inside 
the  nozzle  entrance  where  the  diameters  and  are  verv  small.  Further 
downstream  where  the  nozzle  angle  opens  up  more  rapidly  we  assume  that  a  oounda 
layer  edge  reappears,  leaving  a  core  flow  with  very  little  viscous  dissipation, 
i.e., nearly  isentropic  but  with  a  value  of  entropy  greater  than  the 
starring  value.  This  is  consistent  with  the  Pitot  traverses  normal  to 
the  flow  at  high  pressures  which  have  flat  profilis  near  the  center, 

lndic.it  l:-.?,  flat  velocity  profiles  and  thus  little  viscous  effects  in 
thi;  region.  A:  the  pressure  decreases  the  Reynolds  number  will  decrease 
md  the  boundary  lover  increases.  Urns  the  point  ni  wh  ici  a  boundary  Inver 

•  ■■i.e  reform4  v:I!  occur  progressively  further  down  tee  nozzle. 

r. suiting  in  a  prof  il-  w;th  a  flat  section  tiiat  »a  a  smaller  portion  oL  tile’  t  i . 

•  'or  p  j  "•  2  bar  tie  t  low  remains  fully  develop)  J  along  tv.  entire  noz.’.l 
except  '.or  a  v-ev  short  inlet  r  tiou.  Here  t1  theeret  tea!  calculat  i 
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allow  the  possibility  fur  the  flow  actually  to  go  from  supersonic 
to  subsonic  again  due  to  viscous  effects.  This  is  similar  to  heat  addition 
in  condensing  supersonic  flows.  As  the  nozzle  opens  up  downstream  the 
flow  can  again  become  supersonic. 

It  should  again  be  emphasized  that  the  amount  of  dissipation  seen  in 
these  flows  is  relatively  small,  that  is, in  the  range  of  0.02-0.17,  of  the 
average  dissipation  at  any  given  nozzle  crossection  where  the  flow  is 
fully  developed.  This  picture  is  consistent  with  all  of  the  gasdynamic 
data.  In  addition  the  phase  change  droplet  nucleation  and  growth  calcu¬ 
lations  in  this  program  predict  condensation  in  agreement  with  both 

measurement  of  laser  light  scattering  and  ionization  gage  beam  detection 

3d 

in  the  molecular  beam  configuration.  (There  is  experimental  and  theoretica 

evidence  for  condensation  approaching  1007.  of  the  condensable  species  in 

3 

the  higher  pressure  runs.  ) 

However  for  some  flow  conditions  and  especially  with  light  carrier 

gases  such  as  helium,  viscous  dissipation  can  bo  prevalent  to  such  an 

35 

extent  that  nozzles  become  less  efficient  than  tree  jet  expansions. 

Nozzles  have  been  designed  so  that  nucleation  of  small  mole  fractions  of 

35 

noble  gases  in  a  helium  carrier  gas  can  bo  used  for  beam  experiments. 

In  pure  SF^  expansions  the  boundary  layers  remain  thin  even  in  the  very 

small  entrance  sections  of  several  of  the  nozz’es  so  that  an  isentropic 
core  was  maintained  all  the  way  to  the  exit.  Analysis  of  the  data  using 

one-dimensional  gasdynamics,  along  with  the  calculation  of  thermodynamic 
properties  and  dimensionless  quantities,  has  provided  a  means  for  designing 
small  nozzle  sources  to  encourage  or  indeed  discourage  condensation. 
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TABLE  II 
Gasdynamic  Data 


Nozzle  Gas  Pq  P1  P02  \ 

Bar  Torr  Torr 


Ar 

3 

1.459 

53.84 

4.98 

4 

1.663 

60.65 

4.95 

5 

1.866 

68.93 

4.98 

6 

2.067 

75.40 

4.95 

7 

2.27 

81.87 

4.92 

0.125  SF. 

3 

2.13 

20.78 

2.69 

-He 

4 

2.54 

80.43 

4.94 

5 

2.92 

93.37 

4.96 

SF6 

3 

1.244 

34.79 

4.88* 

4 

1.562 

43.06 

4.84 

5 

1.82 

51.68 

4.87 

6 

2.121 

58.87 

4.86 

7 

2.421 

66.78 

4 , 84 

0.125  SF 

2 

1.31 

23.65 

3.71 

-Ar 

3 

1.584 

33.15 

4.12 

4 

1.866 

45.93 

4.35 

5 

2.113 

56.0 

4.51 

6 

2.303 

63.18 

4.60 

7 

2.63 

66.78 

4.42 

0.03  SF6 

2 

1.236 

36.99 

4 . 58 

-Ar 

3 

1.452 

56.43 

5.23 

4 

1.743 

68.53 

5.2b 

5 

1.881 

75.0 

5.30 

6 

2.107 

82.47 

5.25 

7 

2.332 

87.42 

5.14 

0.0625 

2 

1.235 

2S.68 

4.11 

SFb 
-  A  r 

3 

1.487 

51.68 

3  .ru 

4 

1.690 

6  e  .  0  h 

5.3 

’  Aocau-f  •  v  for  SF,  i  s  a  function  of  temperature  the  first  iteration  lor  J.j, 
is  obtained  usinz  yfTn)  which  is  then  used  to  estimate  exit  tempe ra tu r>.  • 

md  v  ( I  thin  used  to  net  the  final  value  of  M„ . 


TABLE  II 

Gasdynamic  Data  (Cont.) 


Nozzle 

Gas 

Po 

Bar 

P1 

Torr 

P02 

Torr 

5 

1.922 

79.00 

5.49 

6 

2.174 

87.62 

5.44 

7 

8 

2.655 

92.85 

5.06 

9 

2.898 

102.0 

5.08 

He 

3 

2.80 

10.92 

0 . 50 

4 

3.18 

29.11 

2.43 

5 

3.43 

79.0 

3.91 

6 

3.62 

158.  b 

5.42 

7 

3.82 

255.55 

6.72 

8 

4.03 

352.09 

7.69 

■4 

Ar 

2 

0.485 

1.37 

1.25 

3 

0.683 

1.97 

1.27 

4 

0.870 

2.66 

1.32 

5 

0.959 

3.44 

1.45 

6 

1.008 

7.12 

2.12 

7 

1 .06 

10.83 

2.57 

8 

1.11 

15.60 

3.03 

9 

1.18 

19.28 

3.28 

10 

1.24 

22.93 

3 . 50 

SF6 

•4 

0.575 

4.82 

2.68 

5 

0.630 

8.27 

’3 . 38 

b 

0. 713 

12.30 

3 . 8t- 

7 

0.749 

1 4 . 60 

4.10 

8 

0.865 

18.05 

-.24 

o.ohj  sr. 

•4 

C.840 

2.82 

1 .  '  ; 

-Ar 

0.935 

1 

1  .OK' 

"  <■ 

2  v 1 

TABLE  II 

Gasdynamic  Data  (Cont.) 


Nozzle 

Gas 

Po 

P1 

P02 

Bar 

Torr 

Torr 

7 

1.095 

10.57 

2.61 

8 

1.170 

13.45 

2.86 

6 

Ar 

2 

0.939 

33.0 

4.86 

3 

1.093 

44.5 

5.23 

4 

1.215 

55.28 

5.53 

5 

1.326 

59.59 

5.50 

6 

1.422 

65.34 

5.56 

7 

1.501 

71.09 

5.65 

SF6 

2 

0.694 

22.21 

5.19 

3 

0.789 

31.56 

5.76 

4 

0.900 

38.75 

5.97 

5 

1.025 

45.21 

6.03 

6 

1.146 

50.25 

6.02 

7 

1.211 

56.00 

6.17 

0.03  SF6 

2 

0.637 

22.21 

4.95 

-Ar 

3 

1.095 

44.50 

5.35 

4 

1.187 

54.56 

5.70 

5 

1.277 

61.75 

5.84 

6 

1.352 

68.93 

6.00 

7 

1.487 

73.96 

5.93 

0.0625 

2 

0.932 

16.46 

3.57 

SF6-Ar 

3 

1.115 

37.31 

4.95 

4 

1.229 

55.28 

5.75 

5 

1.382 

63.18 

5.79 

6 

1.433 

71.09 

6.04 

7 

1.525 

-9.0 

6.17 

0.125  SF 

2 

0.925 

13.62 

3.93 

-Ar 

3 

1.120 

2  8 .  hS 

4 , 4*4 
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TABLE  II 

Gasdynamic  Data  (Cont.) 


Nozzle 

G3  S 

Po 

P1 

P02 

“r 

Bar 

To  rr 

Torr 

4 

1.275 

38.75 

00 

5 

1.410 

48.81 

5.17 

6 

1.475 

56.0 

5.42 

7 

1.608 

63.18 

5.52 

He 

2 

1.581 

5.68 

1.45 

3 

2.115 

10.0 

1.70 

4 

2.376 

21.50 

2.41 

5 

2.527 

47.37 

3.52 

6 

2.666 

90.50 

4.77 

7 

2.775 

127.87 

5.57 

0.125  SF6 

2 

1.190 

7.12 

2.07 

-He 

3 

1.458 

27.25 

3.78 

4 

1.739 

45.93 

4.51 

5 

1.990 

55.28 

4.63 

6 

2.130 

61.75 

4.73 

7 

2.314 

6  7.50 

4  *  /  4 

Ar 

2 

1.401 

36.40 

4.16 

3 

1.753 

76.12 

3 . 40 

•4 

1.8  73 

80.4  3 

5.3' 

5 

2.083 

37.62 

3.32 

6 

2.288 

102.0 

S .  48 

7 

2.490 

109.18 

S .  *  3 

0.125  SFb 

3 

2 . 344 

18.62 

i  ^ 

.  4  _ 

-He 

■4 

2.834 

3  3.71 

2  . 

=> 

3.272 

■  Q.eS 

4  .0-. 

M 

’ .  t''J  1 

1  1  •>  .  2  ■ 

r3  .1  l 

0.03  .3:  ' 

j 

1  .  -3  i 

'  7 . 2  i 

\  .  1  * 

-Ar 


! 


1  .  -3  i 
l .  :oi 


TABLE  II 

Gasdynamic  Data  (Cont.) 


Nozzle 

Gas 

p0 

Bar 

P1 

Torr 

p02 

Torr 

«R 

4 

1.947 

80.43 

5.40 

5 

2.169 

89.06 

5.38 

6 

2.365 

96.25 

5.36 

7 

2.572 

105.59 

5.38 

8 

2.764 

113.30 

5.38 

0.063  SF, 

2 

1.402 

25.09 

3.59 

-Ar 

3 

1.751 

51.68 

4 . 64 

4 

2.003 

73.25 

5.17 

5 

2.231 

89.06 

5.41 

6 

2.423 

102.00 

5.56 

7 

2.658 

112.06 

5.56 

SF6 

2 

1.106 

27.96 

4.65 

3 

1.41 

39.46 

4.88 

4 

1.728 

48.09 

4.56 

5 

2.066 

56.0 

4.80 

6 

2.437 

64.62 

4.75 

7 

2.868 

71.09 

4.60 

8 

3.220 

79.0 

4.58 

He 

4 

3.641 

20.06 

1.85 

5 

4.000 

44.5 

2.69 

6 

4.245 

90.5 

3.77 

7 

4.467 

162.37 

4.94 

0.125  SF6 

2 

1.409 

24.37 

3.63 

-Ar 

3 

1.790 

37.31 

3.99 

4 

2.099 

44.5 

4.03 

5 

2.427 

54 . 56 

4.15 

6 

2.693 

66.06 

4 . 34 

8- 


TABLE  III 


COMPUTER  SOLUTIONS  FOR  NOZZLE  7-B,  31  SF6-Ar 

Run  No.  p  A  F  x 1  p.  p,.,,  M 

o  p,  p.  1  02 


torr 

torr 

torr 

1 

1 

0.6 

0.0006 

0.88 

1.03 

20.4 

3.64 

2 

1 

0.6 

0.0008 

0.88 

1.15 

20.2 

3.44 

3 

2.2 

0.6 

0.0002 

0.37 

1.46 

n.c. 

4.62 

4 

2.5 

0.58 

0.0002 

0.37 

1.63 

46.7 

4.39 

5 

3 

0.6 

0.0001 

0.37 

1.57 

60.  9 

5.10 

6 

3.5 

0.6 

0.0001 

0.37 

1.72 

71.2 

5.28 

7 

4 

0.6 

0.0001 

0.37 

1.88 

n.c. 

5.40 

S 

5 

0.6 

0.0001 

0.37 

2.21 

n.c . 

5.58 

9 

6 

0.55 

0.0003 

0.37 

2.16 

98.3 

5.19 

10 

7 

0.52 

0.0004 

0.37 

2.57 

101.8 

5.16 

11 

8 

0.52 

0.0004 

0.37 

2.78 

116.6 

5.31 

n.c 


not  calculated 


FIGURE  CAPTIONS 


Fig.  1 


Fig.  2 


Fig.  3 


Fig.  4 


Nozzle  7  is  compared  to  a  free  jet  source  in  a  molecular  beam 
conf igu'ation.  The  relative  beam  intens ity,  Ig ,  primarily  measures 
cluster  intensity  and  shows  the  enormous  advantage  of  the  nozzle 
over  the  free  jet.  (Values  of  Ig  greater  than  100  are  intense 
enough  for  experimental  use.) 

The  experimental  arrangement  permits  nozzles  and  exit  sections 
to  be  interchanged.  Light  scattering  is  used  to  detect  conden¬ 
sation,  M  -  mirror,  L  -  lens,  S  -  slit,  FO  -  fiber  optic,  PMT  - 
photomultiplier  and  BT  -  beam  trap.  The  Pitot  probe  has  x-y-z 
motion  capability. 

The  four  glass  nozzles  discussed  in  this  paper  are  shown  expanded 
8:1  in  diameter  in  order  to  see  the  qualitative  differences.  A 
1:1  contour  of  Nozzle  7  shows  the  actual  geometry.  Three  exit 
sections  A,  B,  and  C  were  attached  to  the  nozzles  to  test  the 
effect  of  nozzle  length  and  to  provide  a  static  pressure  tap  in 
the  case  of  exit  sections  A  and  B  (See  Table  I  for  additional 
details.)  The  upper  part  of  the  figure  shows  the  coordinates 
for  the  solution  to  the  equations  of  motion  plus  the  boundary 
layer  displacement  thickness. 

The  exit  Mach  number  for  Nozzle  7  is  plotted  as  dash  lines  for 

exit  A  and  solid  lines  for  exit  B.  The  Mach  number  M  is  based 

s 

on  the  static  pressure  measurement,  Eqn .  (4),  Mp  is  based  on  the 
Pitot  measurements  and  F.qn.  (5),  and  Mg  is  based  on  both  measure¬ 
ments  at  the  exit  and  Eqn.  (b).  The  geometric  Mach  number  Mg 
is  based  on  area  ratio  only,  i.e.,  no  viscous  effects,  and  is  the 
upper  limit.  The  correct  Mach  number  is  Mg  and  the  deviation 
from  the  other  Mach  numbers  is  due  to  viscous  effects.  Note  the 
closeness  of  the  three  Mach  numbers  in  the  case  of  SF^  indicating 
little  or  no  viscous  dissipation  in  the  "isentropic"  core.  The 
data  provide  increasing  Mach  numbers  in  the  order  Mg,  Ms ,  and  M 
at  anv  given  p  although  the  svmbols  overlap  tor  the  SF^  case. 
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Fig.  5  Pitot  traverses,  for  varying  stagnation  pressure,  are  shown  for 

Nozzle  6  (relatively  large  throat  size)  and  Nozzle  4  (small  throat 
size)  for  helium  and  argon  flows  both  using  exit  configuration  A 
shown  in  Fig.  1.  A  larger  diameter  nozzle  provides  more  intense 
beams  or  higher  exit  Mach  numbers.  Nozzle  4  is  virtually  of  no 

use  with  helium.  (The  nozzle  exit  diameter  is  D  . ) 

e 

Fig.  6  Pitot  traverses  across  the  flow  beyond  the  exit  of  Nozzle  1, 

Exit  B,  for  argon  are  shown  as  a  function  of  position  downstream 
of  the  exit.  Note  the  flat  core  section  close  to  the  nozzle  exit 
which  disappears  progressively  downstream.  Here  x  is  measured 
from  0.5  mm  beyond  the  exit  plane  of  the  glass  nozzle. 

Fig.  7  Pitot  profiles  for  Nozzle  7-B  as  a  function  of  pQ  are  shown  for 

argon,  3%, and  6.3%  SF^  in  argon.  For  one  gas  the  decrease  in  the 
core  size  with  decreasing  pressure  is  due  to  progressively  thicker 
boundary  layers,  i.e.  lower  Reynolds  numbers.  The  variation  in  core 
size  and  centerline  Pitot  pressure  from  one  mixture  to  another  is 
due  to  a  combination  of  variation  in  y  and  6.  The  measurements 
were  made  in  the  exit  plane  of  the  glass  nozzle  (i.e.,  x/Dg  =  -0.25 
in  Fig.  6). 


Fig.  8  Pitot  measurements  are  shown  for  two  free  jets  with  a  factor 

of  25  difference  in  throat  area  and  several  gases.  Mach  number 
from  Eqn.  (5)  is  compared  to  the  theoretical  values  given  by 
Eqn.  (7) .  This  theory  has  been  verified  with  numerous  types  of 
measurements  and  is  applicable  for  the  x/D  >  2. 

Fig.  9  Mach  numbers  are  plotted  with  properties  based  on  the  initial  mole 
fractions  for  mixtures  and  are  nearly  independent  of  pQ  above  3 
bar,  except  for  He  and  He-SF^,  while  SF^  is  constant  over  the 
entire  range.  Nozzle  1  shows  a  slightly  larger  variation  in  exit 
Mach  numbers  for  the  Ar  and  Ar-SF^,  mixtures.  The  He-SF^  mixture 
turns  over  toward  a  flat  value  in  this  nozzle  but  the  pure  helium 
does  not. 

Fig.  10  Tvpical  values  of  Revnold-  number  Re  and  inverse  Knud  sen  nuruor 
'Kn"l  for  Ar  and  SF,  in  Nozzle  7  arc  shown  with  ',f  -  the  throat 

D 

properties  and  c  -  the  exit.  The  characteristic  dimension  for 
noth  parameters  is  the  nozzle  diameter  and  shows  that  the  flaw 
are  still  collision  dominated  even  at  the  lower  pressure^.  A 
wider  variation  exists  for  lie-. 


Fig.  11  The  exic  boundary  layer  thickness,  6,  divided  by  nozzle  radius 
r  is  plotted  as  a  function  of  exit  Reynolds  number.  Pitot 
traverses  provide  6  using  exit  B  or  C.  The  upper  limit  of 
(6/r)g  =  1.0  represents  a  fully  developed  flow.  The  lower  curve 
is  the  displacement  thickness  ( 6^ / r)  obtained  from  M^. 

Fig.  12  The  ratio  M^/M^,  is  an  indication  of  the  boundary  layer  displace¬ 
ment  thickness  and  is  shown  as  a  function  of  Re?  Pr. 

Fig.  13  The  static  and  the  Pitot  pressure  data  for  3%  SF^-Ar  in  Nozzle  7 
are  compared  with  theoretical  solutions  shown  as  solid  lines. 

Fig.  1^  The  parameters  for  the  theoretical  solutions  of  Fig.  14  are 

shown  in  part  A.  Part  B  compares  the  exit  Mach  number  from  the 
computer  calculations,  that  include  the  variation  of  properties 
with  temperature  and  condensation,  to  that  obtained  from  the  data 
and  one-dimensional  gasdynamics  using  isentropic  exit  properties 
and  no  phase  change. 
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FIG. 


beam 


Nozzle  7  is  compered  to  a  free  jet  source  in  a  molecular  boa 
configuration.  The  relative  beam  intensity  1  ^  ,  primarily  m> 
cluster  intensity  and  shows  the  enormous  adv.-int.ic?  of  the  no 
over  th?  free  jet.  (Values  of  Ig  Rreater  than  100  art  inten 
enou'zli  tor  experimental  use.) 


PITOT  TUBE 


The  experimental  arrangement  permits  nozzles  and  exit  sections 
to  be  interchanged.  Light  scattering  is  used  tc  detect  cor.deii 
satior.,  M  -  mirror,  L  -  lens,  S  -  slit,  FO  -  fiber  optic,  PUT 
photomultiplier  and  L'T  -  beam  trap.  The  Pitot  probe  has  x-v-z 
motion  caoabilit'’. 


Geometric 

Throat 


Exit 


The  four  glass  nozzles  discussed  in  this  paper  are  shewn  expanded 
8:1  in  diameter  in  order  to  see  the  qualitative  differences.  A 
1:1  contour  of  Nozzle  7  shows  the  actual  geometry.  Three  exit 
sections  A,  B,  and  C  were  attached  to  the  nozzles  to  test  the 
effect  of  nozzle  Length  and  to  provide  a  static  pressure  tap  in 
the  case  of  exit  sections  A  and  B  <Sete  Table  I  fer  additional 
details;.  The  upper  part  of  the  figure  shows  the  coordinates 
for  One  solution  to  the  equations  of  motion, plus  the  boundarv 
layer  displacement  thickness. 


The  exit  Mach  number  for  Nor.zle  7  is  plotted  as  dashed  linos  for 
exit  A  and  solid  lines  Lor  exit  B.  The  Mach  n\.:nber  M_  is  based 
on  the  static  pressure  measurement,  liqn.  (<*),  M.j  is  based  on  the 
fitot  measurements  and  Eqn.  (5).  and  is  based  on  both  measure¬ 
ments  at  the  exit  and  Eqn.  (o).  The  geometric  ach  number  N. 
is  based  on  area  ratio  onlv,  i.e,  no  viscous  e fleets,  and  is  tin. 
upper  iiinit.  The  correct  Mach  number  is  M.,  and  the  deviation 
irom  the  other  Mach  numbers  is  due  to  viscous  eitocts.  Tote  the 
c loseness  of  the  three  Mach  numbers  in  the  case  of  Si'  indicating 

n 

: . :t  !<•  or  no  viscous  dissipation  in  the  '■  iscntronlc"  core.  1;  1 


data  provide  increasing  Mach  mmirors  in  the  order  M 
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Pitot  traverses,  for  varying  stagnation  pressure,  are  shewn  tor 
Nozzle  '•  (r  2  lat  ive  iy  large  throat  size)  and  Nozzle  4  (small  tnroat 
sizel  for  he  1 in"!  and  argon  tlovs  both  using  exit  ccnf iguration  A 
sntwn  in  Fig.  1.  A  Larger  diameter  nozz Le  provides  mor-  intense 
beams  or  higher  exit  Mach  numbt  r.; .  Nozzle  4  is  virtually  ot  r.o 
us-  with  helium.  (The  nozzle  exit  diameter  is  De.) 
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FIG.  6  Pitot  traverses  across  the  flow  beyond  the  exit  of  Nozzle  1.  Exit 
argon  are  shown  as  a  function  of  position  downstream  of  the  exit. 
Note  the  flat  core  section  close  tc  the  nozzle  er.it  which  J' sap- 
pears  progressively  downstream.  Here-  x  is  measured  ;  .-nm  0..’  mm  h 
the  exit  plane  of  the  glass  nozzle. 
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Pitot  profiles  for  Nozzle  7-d  as  a  function  of  n  are  shew", 
argon ,  3%  and  b .  37>  SFg  in  araoi-,.  For  one  pas  die  deer. if  ; 
core  size  with  decreasing  pressure  is  due  to  progressive!'.'  t 
boundary  layers,  i.e.  lower  Reynolds  numbtrs.  The  .ariati.n 
size  ard  centerline  Pitot  pressure  from  one  mixture  to  aticth 
due  to  •  combination  of  variation  in  y  and  "'ie  measurensi 

nado  in  the  ;xit  plant  of  th.  glass  nozzle  (i.t  ..  x/0t.  =  -0. 


FIC.  10  Typical  values  of  Reynolds  number  He  and  inverse  Knudsen  numbe 
Kn"^  for  Ar  and  SF,,  in  Nozzle  7  are  shown  with  *  -  the  threat 
properties  and  ■  -  the  exit.  The  characteristic  dimension  for 
both  Darn r.e ter s  is  the  nozzle  diameter  and  shows  that  the  flow 
are  still  collision  dominated  even  at  the  lower  pressures,  A 
wiser  variation  exists  for  lie. 


3%SFa-Ar 
6V«%  SFa-Ar 
12’/i%SF6-Ar 
12  Va%  SF6-  Ho 


NOZZLE  7 


Tlu:  exit  boundary  layer  thickness,  5,  divided  by  nozzle  radius 
r  is  plotted  as  a  function  of  exit  Reynolds  number.  Pitot 
traverses  provide  l  using  exit  B  or  C.  The  upper  limit  of 
('/r^  =  1,0  represents  a  fully  developed  flow.  The  lower  cur 
is  the  displacement  thickness  (i,;/r)c  obtained  from  IT  . 
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4  The  parameters  for  the  theoretical  solutions  of  Fis.  14  are 

shcvn  in  part  A.  Part  B  compares  the  exit  Mach  number  from  the 
computer  calculations,  that  include  the  variation  of  properties 
vith  temperature  and  condensation,  to  that  obtained  • rom  the  dat 
and  one -dimensional  pasdvnamics  usinc  isentropic  exit  properties 
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Abstract 

Over  the  past  several  years  an  investigation  into  the 
gasdvnamic  and  nucleation  properties  of  very  small  Laval 
nozzles  has  been  undertaken  in  our  laboratory  due  to  their 
great  potential  use  as  sources  in  cluster  molecular  beams.  A 
series  of  nozzles  designed  specifically  for  helium  carrier  gas 
expansions  has  been  tested  in  our  molecular  beam  apparatus  and 
used  to  study  the  condensation  of  the  noble  gases,  Ar,  Kr,  and 
Xe.  The  goal  of  producing  cluster  beams  with  densities  high 
enough  to  carry  out  high  energy  electron  diffraction  experi¬ 
ments  has  been  attained  for  these  gases  with  mean  cluster 
sizes  in  the  range  of  100-400  ator;  per  cluster.  The  onset  of 
nucleation  appears  to  correlate  with  the  product  of  nozzle 
diameter,  starting  pressure,  and  atomic  potential  well  depth: 
l.e.,  p0DQ  e/k. 

Introduction 

It  is  well  known  that  supersonic  nozzles  can  have  their 
contours  designed  to  control,  within  limits,  the  rate  of 
expansion  and  that  they  can  be  very  much  more  efficient  in 
nucleating  a  particular  gas  than  uncontrolled,  free-jet 
expansions  having  the  same  throat  diameter.^  *  Nozzles  made  o 
varying  throat  diameter,  exit  diameter,  and  nozzle  length  i  Jve 
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been  built  and  tested.  They  have  been  instrumented  so  that  the 
following  data  can  be  obtained:  stagnation  pressure  p  , 
temperature  and  TQ,  nozzle  exit  static  pressure  p. ,  pitot 
pressure  measurements  pQ2  near  the  nozzle  exit,  and  pitot 
traverses  parallel  and  normal  to  the  flow  direction. 

The  experimental  arrangement  is  shown  in  Fig.  la.  The 
nozzle  source  is  mounted  with  xyz  motion  capability  and  is 
shown  with  a  diverging  Laval  nozzle  installed.  If  pitot  (i.e. 
impact)  measurements  are  to  be  made,  the  skimmer  S  is  removed 
and  a  probe  and  pressure  transducer  mounted  so  as  to  separate 
the  first  and  second  pumping  chambers.  Replacement  of  the 
skimmer  re-establishes  the  standard  molecular  beam  configu¬ 
ration.  Hie  cluster  beam  is  crossed  by  a  40  keV  electron  beam 
for  diffraction  studies  and  Debye-Scherrer  patterns  are  taken 
with  a  single  channel,  electron  scintillator,  photon  counting 
dececrion  system  using  phase- sens itiveL detection  and  a  chopped 
molecular  beam. 

The  contours  for  the  nozzles  used  in  this  work  are  shown 
in  Fig.  lb.  The  throat  diameters  vary  from  0.05  to  0.1  mm 
and  exit  diameters  range  from  2  to  5  mm.  The  nozzles  are 
diverging  with  the  minimum  diameter  at  the  entrance.  The  sub¬ 
sonic  flow  ahead  of  the  entrance  is  not  important  for  the 
nucleation  process.  Their  contours  have  been  measured  and 
fitted  to  a  tenth  order  polynomial  as  shown  in  Table  I. 

Nozzle  4  has  the  smallest  divergence  angle  (<1°  total  included 
angle)  and  is  longer  than  the  other  nozzles.  The  noble  gas 
experiments  discussed  here  were  obtained  using  nozzles  11,  12 
and  13.  They  were  designed  to  determine  the  effect  of  throat 
diameter  DQ  and  nozzle  contour  on  their  performance  as  cluster 
sources.  Nozzles  12  and  13  have  the  same  contour  except  near 
the  nozzle  entrance.  Nozzle  11  diverges  more  rapidly  than  12 
and  13  near  the  inlet. 


Gasdynamics 

F.xtensive  gasdynamtc  measurements  using  a  variety  of  gases 
and  gas  mixtures  reveal  that  a  large  fraction  of  the  nozzle 
flow  lies  within  a  boundary  layer  with  some  vis^o^s  dissi¬ 
pation  occurring  all  the  way  to  the  centerline.  ’  Thus,  in 
order  to  correctly  determine  the  local  Mach  number  at  any  point 
in  the  nozzle,  one  must  measure  both  static  and  pitot  pressure 
since  the  total  pressure  in  the  flow  will  decrease  due  to 
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visccus  effects.  The  only  point  where  both  pressures  are 
measured  in  this  work  is  at  the  nozzle  exit.  The  Mach  number 
is  obtained  using  the  so-called  Rayleigh  supersonic  pitot 
equation. 


'.yM-l  R 


(1) 


where  y  is  the  local  value  of  the  specific  heat  ratio,  "lie 
exit  Mach  number  Mp  is  shown  in  Fig.  2  for  several  gas 
mixtures,  all  of  which  are  6  mole  7.  in  a  helium  carrier  gas 
(x  =  0.06).  This  Mach  number  is  seen  to  increase  rr.ona- 
tonically  with  p0 ,  to  values  as  nigh  as  10.  The  three  solid 
curves  are  for  nozzle  12  but  with  different  condensable 
specie^'.'  For  p0  <  4~bar  the  Mach  numbers  are'  nearly  the  same 
but  increase  with  molecular  weight,  i.e.,  from  Ar  to  Kr  to  Xe . 
Above  4  bar  the  curves  reverse  with  Ar  having  the  highest  Mach 
numbers  for  a  given  pQ.  For  Xe  a  comparison  for  the  three 
nozzles  11-13  is  shown.  Nozzle  12  gives  the  highest  Mach 
number  for  a  given  pQ  and  nozzle  11  (same  throat  size  but  more 
rapid  expansion)  is  second  highest.  Nozzle  13  (same  contour 
as  12  but  smaller  entrance  diameter  DQ)  exhibits  the  lowest 
Mach  numbers. 


With  the  system  arranged  in  its  molecular  beam 
configuration,  beam  Intensities  were  measured  using  an  ioni¬ 
zation  gage  (IG  in  Fig.  1).  The  nozzle-to-skimmer  distance,  xg 
was  varied  to  determine  its  effect  on  beam  intensity.  Pure 
argon  expansions  are  shown  in  Fig.  3.  Note  that  intensity 
variations  with  xgDe  are  not  large.  There  is  an  interesting 
variation  in  maximum  beam  intensity  going  from  a  downstream 
peak  location,  double  peak,  and  upstream  location  as  pQ  is 
increased,  analogous  to  some  free  jet  cluster  source  data.^>^ 
All  beam  experiments  were  carried  out  with  xg  =  2.5  mm.  This 
gives  a  value  of  xg/D^  **  0*5  for  nozzle  13  which  is  optimum  for 
the  high  pp  beams  as  seen  in  Fig.  3. 

Although  Laval  nozzles  are  much  more  efficient  cluster 
sources  than  free  jets  for  most  flows,  care  must  be  taken  in 
their  design  or  they  can  actually  produce  beam  intensities 
lower  than  a  free  jet  source.  Consider  for  example  the 
expansion  of  SF^  (\  =■  0.1Z5)  through  nozzle  4  shown  in  Fig.  4. 
When  expanded  in  an  argon  carrier  gas  the  beam  intensity 


increases  dramatically  near  pQ  =  2.5  bar  due  to  clustering. 
However,  for  helium  carrier  gas  (same  specific  heat  ratio  y) 
the  beam  intensities  remain  quite  low  even  out  to  pQ  =  8  bar. 
Thus  this  nozzle  is  of  little  use  for  helium  carrier  gas, 
condensation  experiments  with  TQ  near  room  temperature. 

Because  a  high  y  carrier  gas  is  required  for  condensation 
of  a  low  y  gas  in  an  adiabatic  expansion,  and  because  the  Laval 
nozzles  control  the  rate  of  expansion  but  cannot  limit  the 
maximum  Mach  number  in  these  expansions,  He  is  desirable  since 
it  is  unlikely  to  supersaturate  and  cause  mixed  condensable- 
carrier  gas  clusters.  The  use  of  pure  gas  expansions  is  also 
an  option,  provided  y  is  not  too  near  1,  but  does  not  allow  for 
the  possibility  of  controlling  the  cluster  size  distribution 
or  temperature  that  is  potentially  available  using  a  carrier 
gas.  In  addition  there  is  the  problem  of  the  high  cost  of  some 
pure  gases  such  as  Xe.  As  a  result  of  these  criteria  and 
prior  experience  with  nozzles  having  long,  small  divergence 
inlets,  e.g.,  nozzle  4,  several  shorter  nozzles  that  diverge 
more  rapidly  have  been  constructed.  ^ ^  _  _ 

There  are,  however,  a  number  of  problems  using  He  rather 
than  other  commonly  used  carrier  gases  such  as  Ar  or  N2. 

Boundary  layer  effects  are  more  severe  due  to  its  higher  kine¬ 
matic  viscosity.  In  addition  there  is  the  problem  of  reduced 
high  vacuum  diffusion  pump  performance.  (i.e.,  for  nozzles 
with  E0  >  0.1  mm  pressure  fluctuations  of  an  order  of  magnitude 
Or  more  occurred  in  the  second  and  third  stage  pumping  chambers 
when  p  >  2  bar.)  Thus  in  order  to  operate  at  sufficiently 
high  Pq  for  cluster  beam  generation,  it  is  necessary  to 
decrease  the  nozzle  throat  size  below  0.1  mm.  This  exacerbates 
the  problem  of  helium  boundary  layer  growth  so  that  nozzle 
design  must  be  carefully  considered. 

Nucleation 

Measured  molecular  beam  intensities  for  He-carried  Ar ,  Kr , 
and  Xe  with  xQ  *  0.06  are  shown  in  Fig.  5  for  nozzle  12  and 
Tq  =  295.  As  one  might  expect,  condensation  occurs  in  the  orctr 
Xe ,  Kr,  and  Ar  as  p  is  increased.  All  the  molecular  beam 
intensity  data  are  given  in  fixed,  arbitrary  units  (a.u.  in 
Figs.  3-7)  with  10  approximately  equal  to  2  x  10^®  atoms/ 

(cm^  s)  at  the  electron  beam  location.  Once  the  beam  intensity, 
I.,  exceeds  unity  (a.u.)  the  beam  density  is  sufficient  for 
electron  diffraction  measurements.  Because  all  three  of  these 
gases  are  noble  gases,  one  is  tempted  to  look  for  some  universa 1 
behavior.  Thus  normalizing  to  Xe,  these  three  curves  are 
replotted  as  dotted  lines  with  the  ordinate  equal  to  pQ  e/k 


where  e  Is  the  interatomic  potential  well  depth^  and  k 
Boltzmann's  constant.  When  replotted  the  onset  of  condensation 
occurs  at  the  same  value  of  p0  e/k  for  the  three  gases. 

A  plot  of  these  data  as  a  function  of  a  nondimensional 

or  reduced  pressure  p0/(e/o^)  or  Pv/(je/o'^),  where  a  Is  the 
radius  for  zero  interaction  potential®  and  py  is  the  conden¬ 
sable  vapor  partial  pressure,  provided  no  unification  of  the 
onset  point. 

For  a  given  species  and  mole  fraction  nucleaticn  results 
from  one  nozzle  to  another  will  be  compared.  The  results  for 
Ar,  x  =  0.25,  for  nozzles  11-13  are  shown  in  Fig.  6.  The 
solid  lines  are  for  the  data  plotted  as  a  function  of  pQ  in 
which  the  onset  of  nucleatiou,  for  the  nozzles  with  the  two 
largest  diameters,  occurs  first.  Nozzle  12,  with  a  lover 
e1tpatTSlon-rate  than  ~noz2le  11 .  exhfblts'more  rapid  growth  of 
the  condensed  ptiase  after  the  point  of  onset  since  the 
measured  beam  intensity  beyond  onset  is  due  entirely  to  the 
condensed  phase.  As  with  free  jet  cluster  sources,  these  data 
have  also  been  plotted  vs  pG  Dq  (dashed  lines)  such  that  the 
curve  for  nozzle  12  is  unchanged.  Recall  that  p  D  is 
proportional  to  the  total  number  of  binary  collisions  per  atom 
for  a  given  expansion.  Here  again  it  is  interesting  to  note 
that  detectable  condensation  onset  occurs  at  a  single  point  on 
the  pQ  Dq  scale.  Note,  however,  that  the  three  curves  would 
aJso  converge  to  a  common  onset  point  were  they  plotted  versus 
PoDo  wllich  f3  proportional  to  the  three-body  collision 
frequency  required  to  form  a  stable  dimer  in  the  initial 
cluster  formation  process.  These  results,  although  interesting 
are  not  universal,  as  seen  in  Fig.  7  for  Xe,  x  =  0.06 
expansions.  In  this  case  the  onset  of  nucleation  in  nozzle  11 
occurs  at  lower  pe  than  that  of  nozzle  12  and  their  growth 
curves  cross  over  at  a  pressure  above  onset.  When  replotted  as 
a  function  of  pQ  DQ  the  results  from  nozzles  12  and  13  (which 
have  the  same  nozzle  contour  but  different  DQ)  fall  nicely  on 
top  of  one  another  whereas  the  results  from  nozzle  11  move 
slightly  further  away  (see  Fig.  7).  Thus  the  correlation  with 
hinary  or  trimolecular  collision  frequency  is  not  perfect,  but 
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then  the  various  nozzles  have  flow  histories  which  are  not  u.-- 
universai  i«  t!ios>:  of  the  free  jet  expansions. 

Electron  Diffraction 

Electron  diffraction  patterns  have  been  taken  with  our 
single  channel  detection  system  for  all  three  condensable 
noble  gases  under  a  variety  of  conditions  and  using  several 
nozzles.  An  average  cluster  j ize  can  be  estimated  from  the 
measured  broadening  of  the  diffraction  rings.  Knowing  the 
density  of  the  condensed  phase,  the  number  of  atoms  per 
cluster  may  then  be  calculated.  The  results  for  Xe ,  v  =  0.06 
are  shown  in  Fig.  8.  The  sizes  are  in  the  g  =  100-1000  atoms 
per  cluster  range.  The  sizes  of  Fig.  8  correlate  with  the 
relative  intensities  of  Fig.  7  in  that  a  given  and  g  occur 
for  progressively  higher  pD  in  the  order  nozzle  11,  12,  and  13. 
recall  that  g  conies  from  a  diffraction  peak  which  is  a  sum  of 
the  diffraction  from  a  large  number  of  clusters  that  cross  the 
electron  beam.  There  is  always  a  distribution  of  sizes  in 

the  Seam,  and  since  the  scattering  intensity  at  a  diffraction 
peak  coos  -’s  the  square  of  the  number  of  atoms  in  a  cluster, 
the  larger  clusters  scatter  more  electrons.  Thus  g  is  weighted 
toward  the  high  end  of  the  size  distributions.  Theoretical 
mode!  for  the  noble  gas  clusters  in  this  size  regime  predict 
the  possibility  of  icosahedral  packing  in  contrast  to  the  bulk 
FCC  structure.'’^  The  comparisons  we  have  made  between  datJ 
and  cluster  models  have  so  far  shown  better  agreement  using 
icosahedral  rather  than  the  bulk  structure. 

Conclusions 

We  have  arrived  at  the  following  conclusions  with  regard 
to  snail  Laval  nozzle,  noble  gas  clustering: 

II  Noble  gases,  even  when  a  small  mole  fraction  expan  led 
In  helium,  readily  nucleate. 

2)  Because  of  the  high  kinematic  viscosity  of  helium  the 
nozzle  contour  must  be  carefullv  designed  to  balance  oft  the 
effects  of  increased  boundary  layer  growth  and  the  requirement 
that  throat  diameter  be  small  to  avoid  pumping  problems. 

31  There  is  tantalizing  evidence  that  a  parameter  p  P 
$/k  or  py  D  */k  '-ill  provide  correlation  or  scaling  wi’h 
regard  to  the  onset  o:  noble  gas  nucleatlon  in  La^al  nozzle 
iouro  . 

.)  Cluster  beam  intensities  (actually  cluster  d.nsiLios) 
hi  h  nov/h  to  conduct  "lectrou  diffraction  can.  be  attained  in 
a  size  regime  of  great  i merest  for  study  of  c! 


e  structure 


the  condensed  phase  and  possible  deviations  from  bulk  matter. 
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Table  1  Nozzle  Geometry 


a) 


Nozzle 


Number 

4 

11 

12 

13 

Vcn 

0.0063 

0.00S9 

0.0089 

0.0056 

0e-cm 

0.256 

0.47 

0.47 

r  .47 

Lr  cm 

.  3.15 

1.176 

t  .  ..  1.356  .  . 

1.233  . 

Jo 

0.00635 

0.0088 

0.009 

0.00 >5 

al 

-0.0248 

-0.0150 

0.015C 

-.1402 

<10 

0 • 1488 

4.6518 

-1.62C8 

2.473 

a  3 

-0.1156 

-43.562 

23.052 

-13.2382 

1  - 

0.0450 

+130.334 

-120.408 

39.233 

a5 

-0.006J7 

-361.34 

337.06 

-63.9J6 

.1 

0 

+378.83 

-561.0 

59- o2> 

■1 

0 

-201.50 

579.44 

-29.72',- 

ab 

0 

43.006 

-366. 

6.0812 

a  0 

0 

0 

129.404 

0 

■Ho 

0 

0 

-19.6 

0 

0 


A  best  fit  polynomial  equation  of  the  form, 
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D(X)  =  ai  X 

f.  lifted  to  measured  contours  with  D  and  ^  in  cm. 
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A 


B 


Fig.  1  The  experimental  arrangement  shows  the  nozzle  and 
stagnation  chamber  movable  in  xyz  space  with  either  a  skimmer 
S  or  a  pitot  probe  and  transducer  located  downstream  of  the 
nozzle  exit.  With  the  skimmer  in  place  and  collimator  C  a 
molecular  beam  is  crossed  with  a  40  keV  electron  beam  which  i 
trapped  in  a  beam  trap  BT.  Diffraction  patterns  arc  taken  wi 
a  detector  Det.  The  nozzle  contours  are  shown  in  lb  with 
additional  details  in  Table  I. 


NOZZLE 
12 
1 1 


6  *  in  He 


Fig.  2  Typical  nozzle  exit  Vich  r.^.rDers  obtained  using 
are  shown  for  nozzle  12  ns  the  solid  lines  fo*“  Ar,  Kr, 
and  as  dashed  lines  for  Xe  in  nozzles  11  and  13.  All 
expansions  begin  from  room  temperature  T0  =  293  K  and 
in  a  He  carrier  gas. 
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Xs/De 


Fig.  3  Molecular  beam  intensities  1^  are  shown  as  a  function  of 
distance  from  the  nozzle  exit  to  the  skiusrer.  The  location  of 
maximum  beam  intensity  is  indicated  with  an  arrow. 


Nozzle  12 
6  X  in  He 


13,6/ k  (Bar-K) 


Fig.  5  Bean  intensities  as  a  function  of  pQ  are  shown  as  the 
solid  lines  for  nozzle  12  with  *  0.06  and  a  He  carrier  gas 
Vith  a  scale  change  to  pQ  pe/k,  normalized  to  Xe,  the  curves 
for  Ar  and  Kr  move  to  the  dashed  lines  as  indicated  by  the 


arrows . 


2  3  4  5  6 


■P0  ;  Bar. 


Jig.  8  The  average  number  of  atoms  per  cluster  g  is  si 
function  of  pQ  for  Xe ,  x0  3  in  a  He  carrier  '.as 

nozzles  11-13. 
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